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PREFACE 


This  report  contains  papers  on  various  wind  tunnel  correction  methods  used  in  high 
angles  of  attack  tests.  The  papers  were  solicited  from  the  various  NATO  countries  and 
presented  in  a  round  table  discussion  following  the  AGARD  Fluid  Dynamics  Panel 
Symposium  in  Neubiberg,  Germany  in  May  1980.  Papers  given  and  published  here  are  from 
Canada,  France,  Germany,  Netherlands,  Sweden,  United  Kingdom  and  the  United  States. 

Several  wind  tunnel  wall  correction  methods  in  use  or  under  study  are  presented  here 
for  closed,  open  and  ventilated  wall  wind  tunnels.  The  Mach  number  range  is  generally 
limited  up  to  high  subsonic  speeds  with  some  techniques  only  useful  for  incompressible 
flow.  Wall  correction  techniques  discussed  along  with  their  attributes  and  disadvantages 
include  vortex  lattice,  panel,  system  of  images,  wall  pressure  and  adaptive  walls.  The 
adaptive  wall  technique  is  a  method  to  actively  reduce  or  eliminate  the  need  for  wall 
correction  and  is  becoming  more  favorable  as  d  welopment  problems  are  solved. 

Current  testing  in  wind  tunnels  over  this  speed  range  should  include  documentation  of 
pressures  and/or  flow  angles  in  the  inviscid  flow  near  the  tunnel  walls.  Such  measurements 
facilitate  (1)  determination  of  magnitude  of  w.U  interference,  (2)  determination  of 
equivalent  free-flight  test  conditions,  (3)  correction  for  wall  interference  and 
(4)  calculations  for  adaptive  wall  adjustments  in  t  innels  which  have  variable  geometry. 
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CANADIAN  STUDIES  OP  WIND  TUNNEL  CORRECTIONS 
FOR  HIGH  ANGLE  OF  ATTACK  MODELS 
M.  Mokry 

National  Research  Council  Canada 
National  Aeronautical  Establishment 
Ottawa,  Ontario 


This  paper  is  a  brief  account  of  wind  tunnel  interference  studies  relating  to  testing 
of  high  angle  of  attack  models,  carried  out  in  Canada  during  the  last  decade.  If  any 
relevart  Canadian  work  is  omitted  here,  it  is  only  because  the  author  was  unaware  of  its 
existence.  The  work  of  other  AGARD  countries  is  referenced  here  only  when  needed  to  indi¬ 
cate  the  origin  of  an  idea  or  parallel  development. 

In  the  early  seventies,  the  problem  of  wind  tunnel  interference  on  high  angle  of 
attack  models  was  dealt  with  at  Canadair,  in  connection  with  the  research  program  aimed 
at  the  development  of  high  lift  technology  for  advanced  aircraft.  A  vorticity  distribu¬ 
tion  method1  for  the  calculation  of  potential  flow  past  sulti-component  airfoils  in  free 
air  was  extended2  to  flow  past  a  multi-component  airfoil  in  an  infinite  cascade  and  used 
to  evaluate  wall  interference  effects  due  to  solid  wind  tunnel  walls.  Figure  1,  taken 
from  Ref.  2,  shows  schematically  the  generation  of  an  infinite  cascade  used  to  model  the 
flow  past  a  four-component  airfoil  between  two  wind  tunnel  walls.  The  computation,  utiliz¬ 
ing  vorticity  panels  in  the  infinite  cascade  mode,  is  performed  on  an  eight-component 
configuration,  consisting  of  the  original  airfoil  and  its  mirror  image.  Unlike  the  usual 
approximations  of  wall  effects  by  a  few  neighbouring  images,  this  solution  is  exact  (in 
the  limit  of  decreasing  panel  lengths) . 

The  computed  pressure  distributions  for  a  40  deg  flap  deflection  are  shown  in  Fig.  2; 
the  sc  lid  line  is  the  wind  tunnel  computation  and  the  dashed  line  is  the  free  air  computa¬ 
tion.  Also  shown  for  comparison  is  the  experimental  pressure  distribution  on  a  2  ft  chord 
model  in  the  NAE  6  ft  x  9  ft  low  speed  wind  tunnel.  The  theoretical  curves  indicat  t  that 
there  should  be  a  negligible  wall  effect  on  lift  at  this  flap  deflection.  The  discrepancy 
betWvien  experiment  and  theory  is  apparently  caused  by  viscous  and  three-dimensional  effects 
which  had  not  been  accounted  for  in  the  method. 

Figure  3  shows  the  result  for  the  same  airfoil  with  a  70  deg  flap  deflection.  The 
potential  flow  theory  shows  a  4%  loss  in  lift  due  to  wind  tunnel  interference.  The 
difference  between  experimental  and  calculated  pressures  is  more  pronounced  than  in  the 
previous  case,  since  the  configuration -is'  already  in  stall  (beyond  CLmax) • 

Using  the  methods  based  on  viscous  -  inviscid  interactions,  e.g.  Ref.  3,  it  should  be 
possible  to  evaluate  the  effects  of  airfoil  boundary  layers  on  the  pressure  distributions 
and  establish  a  better  agreement  between  theory  and  experiment.  However,  the  press:  ce  of 
boundarv  layers  on  the  test  section  walls  is  an  obstacle,  since  then  the  walls  do  noi 
behave  like  perfect  reflection  planes  and  the  concept  of  images  would  have  to  be  moditied, 
if  not  abandoned  altogether. 

At  the  National  Aeronautical  Establishment  (NAE) ,  the  method  of  singularity  distribu¬ 
tions  was  later  extended  to  compute  wall  interference  effects  in  perforated  wall  test 
sections,  assuming  a  linear  relationship  between  the  pressure  coefficient  and  the  normal 
component  of  velocity  at  the  wall.  The  solution  of  the  problem  was  separated  into  two 
parts:  first,  the  construction  of  modified  singularities  (influence  or  Green's  functions) 

satisfying  the  wall  boundary  conditions  and,  second,  the  solution  for  contour  distributions 
of  these  singularities  to  satisfy  the  airfoil  boundary  conditions,  in  a  procedure  equiva¬ 
lent  to  that  used  in  free  air.  The  computation  by  a  panel  method  is  thus  performed  only 
on  the  original  airfoil;  no  panels  are  placed  on  the  walls.  As  the  previously  discussed 
cascade  approach,  this  method  is  again  exact.  The  case  of  solid  walls  is  obtained  as  a 
special  case  of  zero  porosity.  The  source  distribution  version  of  the  method  is  described 
in  Refs.  4  and  5;  the  vortex  distribution  version  in  Ref.  6. 

Unfortunately,  experience  with  perforated  walls  indicates  that  the  normal  (cross-flow) 
component  of  velocity  is  a  nonlinear  function  of  the  local  pressure  coefficient  and,  there¬ 
fore,  the  applicability  of  wall  interference  corrections  evaluated  from  linear  porous 
wall  theory  is  only  limited.  Particularly  on  the  wall  facing  the  suction  side  of  the 
airfoil  the  normal  velocity  is  strongly  modulated  by  the  rate  of  growth  of  the  boundary 
layer,  as  discovered  in  Ref.  7.  Figure  4  gives  an  example  for  a  20%  perforated  wall  of 
the  15  in  x  60  in  test  section  of  the  NAE  high  speed  wind  tunnel.  The  flow  disturbances 
at  the  walls  were  generated  by  a  10  in  chord  BGK  1  airfoil,  located  midway  between  the 
walls.  Investigated  is  the  boundary  layer  on  the  wall  facing  the  suction  side  of  the  air¬ 
foil;  on  the  other  wall  (pressure  side)  the  boundary  layer  was  almost  completely  bled 
into  the  plenum.  The  wall  pressure  coefficient  Cp  ,  measured  by  a  static  pressure  pipe, 
is  plotted  against  the  ratio  v0/ue  of  the  normal  and  tangent  velocities  on  the  outer  edge 
of  the  boundary  layer.  The  normal  velocity  is  deduced  from  the  boundary  layer  develop¬ 
ment,  computed  according  to  the  method,  of  Ref.  8  which  can  handle  suction  and  injection, 
matching  the  velocity  profiles  measured  by  boundary  rakes  at  several  stations  along  the 
wall.  From  Fig.  4  it  is  sufficiently  clear  that  a  linear  relation,  used  in  porous  wall 
theory,  does  not  adequately  represent  the  basic  nonlinear  characteristics  of  the  cross-flow. 
The -curve  shapes  and  the  dependence  on  Mach  number  M  and  geometry  of  the  tested  airfoil 
(here  represented  by  the  angle  of  attack  a)  are  in  general  agreement  with  observations  in 
Ref.  9  (U.S.A.),  based  on  direct  measurements  by  a  laser  velocimeter. 
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Besides  including  the  computation  of  wall  boundary  layers  as  an  integral  part  of  wall 
correction  procedure,  there  are  several  ways  of  avoiding  the  difficulties  with  complicated 
boundary  conditions  of  ventilated  wind  tunnel  walls.  Two  of  them  have  been  extensively 
studied  in  Canada  and  are  discussed  here: 

(1)  development  of  wall  correction  procedure  which  uses  boundary  measurements,  but 
does  not  require  the  knowledge  of  cross-flow  properties,  and 

(2)  development  of  ventilated  walls  which  behave  in  a  well  predictable  manner. 

The  first  approach,  based  on  the  measurement  of  static  pressures  on  the  test  section 
boundary  was  first  successfully  demonstrated1 11  at  ONERA,  France.  Assuming  that  the  flow 
near  test  section  walls  is  subsonic  linearized,  it  is  possible  to  estimate  the  far  field 
of  the  modal  and  separate  from  the  wall  pressure  disturbance  the  part  which  is  due  to  wall 
interference.  Since  the  linearized  pressure  coefficient  is  proportional  to  the  axial 
(streamwise)  component  of  disturbance  velocity,  the  velocity  correction  at  the  model  posi¬ 
tion  can  be  obtained  as  a  solution  of  the  Dirichlet  problem  in  the  test  section  interior. 

The  cross-flow  properties  of  the  walls  thus  do  not  enter  the  problem. 

The  component  of  wall  interference  velocity  normal  to  the  tunnel  axis  is  obtained  as 
the  harmonic  conjugate  of  the  axial  component.  Since  the  integration  of  Cauchy-Riemann 
equations  involves  an  arbitrary  constant,  the  determination  of  the  angle  of  attack  correc¬ 
tion  requires  specification  of  the  flow  angle  at  an  arbitrary  distant  point  within  the 
test  section.  A  simple  approach  is  to  assume  that,  the  flow  enters  the  test  section  in  the 
direction  parallel  to  the  wind  tunnel  axis,  but  for  a  more  reliable  evaluation  of  the  angle 
of  attack  correction  an  actual  flow  angle  measurement  should  be  carried  out. 

The  difference  between  the  ONERA  method  and  the  method  developed  at  NAE,  Ref.  11, 
is  in  the  selection  of  test  section  boundaries  ar d  the  construction  of  the  solution.  The 
ONERA  method  utilizes  the  concept  of  an  infinitely  long  test  section,  Fig.  5a.  However, 
since  the  pressures  can  be  measured  only  over  a  finite  length  of  the  test  section,  the 
boundary  values  have  to  be  extrapolated  to  infinitely  distant  upstream  and  downstream 
points.  In  Ref.  10,  both  components  of  wall  interference  velocity  are  obtained  in  one 

step,  using  the  conformal  mapping  of  the  infinite  strip  onto  the  upper  half-plane,  wher.- 

the  solution  is  given  by  Cauchy's  integral.  The  angle  of  attack  correction  at  the  position 
of  the  model  is  obtained  assuming  the  flow  parallel  to  the  wind  tunnel  axis  at  the  infinitely 
distant  upstream  point. 

The  NAE  method1 1  is  based  on  the  finite  test  section  concept  and  utilizes  the  Fourier 
solution  of  the  Dirichlet  problem  in  a  rectangle.  Fig.  5b.  This,  in  principle,  requires 

that  the  static  pressure  be  measured  along  the  ups  cream  and  downstream  sides  of  the  rec¬ 

tangle  in  addition  to  the  floor  and  ceiling  pressures.  However,  if  the  latter  are 
measured  sufficiently  far  upstream  and  downstream  of  the  model  location,  a  simple  inter¬ 
polation  between  the  floor  and  ceiling  pressures  will  suffice.  The  evaluation  of  the 
coefficients  of  the  infinite  series  solution  by  the  fast  Fourier  transform  makes  the  method 
very  efficient  and  suitable  for  routine  correcting  of  wind  tunnel  test  data. 

Examples  of  the  evaluated  wall  interference  corrections  from  measured  wall  pressures 
and  forces  on  a  high  angle  of  attack  model  are  given  in  Figs.  6-8,  using  unpublished  test 
data  by  F.C.  Tang  of  NAE.  The  static  pressures  along  the  20%  perforated  walls  of  the 
15  in  x  60  in  test  section  of  the  NAE  high  speed  wind  tunnel  were  measured  by  1  in  diameter 
static  pressure  pipes,  attached  to  the  centerlines  of  the  walls.  The  force  and  pitching 
moment  coefficients  were  measured  by  a  three-component  balance  and  compared  with  the  values 
obtained  by  surface  pressure  integration.  In  addition,  the  wake  drag  was  measured  by  the 
pitot-traverse  method.  The  tested  12  in  chord  airfoil  was  equipped  with  a  4.5  in  flap  at 
the  deflection  angle  of  35  deg,  zero  overlap  and  gap  of  3%  of  the  chord  length. 

In  the  coordinate  system  of  Figs.  6-8,  the  quarter-chord  of  the  airfoil  is  at  x/c  =  0 
and  CP  is  the  wall  pressure  coefficient.  The  column  of  measured  quantities  contains:  Rc 
the  Reynolds  number  based  on  chord  length,  M  the  tunnel  Mach  number,  a  the  geometrical 
angle  of  attack,  CN,  CX,  and  CM  the  normal  force,  axial  force,  and  pitching  moment  coeffi¬ 
cients  obtained  by  balance  measurements,  A  the  cross-sectional  area  of  the  airfoil  normal¬ 
ized  by  c* ,  and  c/h  the  chord  to  test  section  height  ratio.  In  the  column  of  computed 
quantities  are:  AM  the  Mach  number  correction,  Aa  the  angle  of  attack  correction,  c*3AM/3x 
and  c*3Aa/3x  the  variations  of  AM  and  Aaover  the  chord  length,  and  Mco_  and  acor  are  the 
corrected  values  of  M  and  a.  The  values  CLcor,  CDcor,  and  CMcor  are  the  corrected  lift, 
drag,  and  pitching  moment  coefficients,  obtained  from  the  above  corrections  to  stream 
quantities  by  standard  techniques.  The  chordwise  variations  of  AM  and  Act  are  small  in  all 
three  cases,  indicating  that  the  airfoil  was  tested  in  a  reasonably  uniform  flow  environ¬ 
ment,  corresponding  to  free  air  conditions  at  _  and  a  _ . 

cor  cor 

The  wall  pressures  in  Fig.  6  were  obtained  with  the  airfoil  at  low  incidence,  in  Fig.  7 
at  high  incidence  near  and  in  Fig.  8  in  deep  stall.  The  lift  coefficients  in  Figs.  6 

and  8  are  nearly  the  same  and  so  is  the  magnitude  of  wall  pressure  coefficients.  However, 
contrary  to  the  expectations  from  classical  wall  interference  theory,  the  values  of  Aa 
differ  markedly. 

Figs.  9,  10,  and  11  show  the  corrected  and  uncorrected  CL  vs.  a  curves,  the  drag  polars, 
and  the  pitching  moment  polars  respectively.  The  agreement  between  the  lift  coefficients 
obtained  by  the -balance  measurement  and  pressure  integration  in  Fig.  9  is  excellent. 

However,  the  correlation  of  drag  polars  in  Fig.  10  is  poor.  It  has  been  suggested  in  Ref. 

12,  that,  unlike  the  wake  drag  which  reflects  only  the  viscous  effects,  the  balance  and 
pressure  drags  are  subjects  to  inviscid  forces  which  are  due  to  the  momentum  transfer 
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through  the  test  section  boundaries. 


Besides  the  analysis  of  drag  forces  on  high  angle  of  attack  airfoils,  the  work  at  NAE 
is  also  continuing  on  extensions  of  the  method  of  Ref.  11  to  three-dimensional  flows. 
Considered  is  the  evaluation  of  wall  interference  corrections  using  the  static  pressure 
pipes  attached  to  the  test  section  floor  and  ceiling  (two-rail  arrangement)  and  also  the 
additional  pipes  attached  to  the  sidewalls  (four-rail  arrangement) . 

The  task  of  developing  a  new  type  of  r.  ventilated  test  section  was  undertaken  in  the 
Aerodynamics  Laboratory  of  the  Department  ci  Mechanical  Engineering  at  the  University  of 
British  Columbia  (UBC) .  A  paper1 *  at  the  seventh  ICAS  Congress  described  the  unsatisfactory 
performance  of  the  longitudinally  slotted  walls  for  high  lift  airfoil  testing.  Similarly 
to  the  case  of  perforated  walls,  one  of  the  most  disturbing  features  of  slotted  walls  is 
the  dependence  of  the  theory  on  empirical  factor  which  varies  not  only  with  the  wall  con¬ 
figuration,  but  also  with  the  test  model.  Recently,  more  detailed  studies  of  flow  through 
slotted  walls  were  published  in  Ref.  14  (Sweden)  and  Ref.  15  (U.S.A.). 

The  unsatisfactory  empiricism  of  the  longitudinally  slotted  and  porous  walls  is  a 
result  of  the  partly  separated  flow  through  the  walls.  Therefore,  it  was  proposed16'17 
to  design  a  ventilated  wall  consisting  of  transverse  slats,  shaped  as  symmetrical  airfoil 
sections,  Fig.‘  12,  which  would  operate  within  their  unstalled  incidence  range.  A  singu¬ 
larity  distribution  method,  such  as  the  one  described  in  Ref.  18,  is  a  natural  tool  for 
studying  the  effects  of  different  wall  configurations.  The  boundary  conditions  to  be 
satisfied  are  tangent  flow  on  the  solid  wind  tunnel  wall,  the  airfoil  wall  slats,  and  the 
test  airfoil,  Kutta  conditions  at  the  trailing  edges  of  all  airfoil  elements,  and  a  free 
streamline  pressure  condition  on  the  shear  la/er  generated  in  the  plenum.  However,  tue 
potential  flow  calculations  serve  merely  as  a  guide  to  a  suitable  wall  configuration, 
which  then  has  to  be  verified  experimentally. 

The  reasons  for  using  the  transverse  slotted  wall  on  the  suction  side  of  the  airfoil 
only  are  following: 

(n)  The  pressure  distribution  on  the  pressure  side  of  a  high  lift  airfoil  is  aominated 
by  stagnation  pressure,  see  Figs.  2  and  3,  and  should  undergo  negligible  change 
from  free-air  to  tunnel  conditions  when  the  opposing  wall  is  solid. 

(b)  On  the  suction  side,  the  shear  layer  produced  at  the  interface  between  the  tunnel 
airflow  and  the  stagnant  plenum  air  enter  the  test  section  downstream  of  th': 
model,  see  Fig.  13.  If  the  other  wall  were  slotted  as  well,  the  stagnant  plenum 
air  would  be  brought  into  the  test  section  upstream  of  the  model,  and  this  air 
plus  its  bounding  shear  layer  would  degrade  the  test  section  flow  near  the 
airfoil. 

(c)  By  analogy  to  the  theoretical  test  section  having  one  wall  solid  and  the  other 
formed  by  an  open  jet  boundary,  th.'  angle  of  attack  correction  is  expected  to  be 
minimized. 


The  last  point  brings  us  to  the  most  important  objective  of  this  work,  namely  the 
development  of  a  test  section  producing  adequately  low  corrections  to  test  data  for  a 
wide  range  of  sizes,  shapes  and  angles  of  attack  of  test  airfoil.  For  many  purposes 
"adequately  low"  would  mean  corrections  of  the  order  of  1%  to  pressure,  force  and  moment 
coefficients,  and  it  should  be  possible  to  maintain  such  low  corrections  for  relatively 
large  models  at  high  angles  of  attack. 


An  experimental  verification  of  the  above  test  section  concept  was  demonstrated  in 
the  388  mm  x  915  mm  insert  of  the  UBC  low-speed  closed-circuit  wind  tunnel.  One  wall  was 
surrounded  by  plenum  and  could  be  fitted  with  airfoil  shaped  slats  of  NACA  0015  section 
and  chords  of  46  mm  (small  slats)  or  92  mm  (large  slats) ,  at  zero  incidence.  A  full  range 
of  wall  open  area  ratios  (OAR)  could  be  tested.  The  test  airfoils  were  mounted  on  the 
turntable  of  a  six-component  balance.  Figure  14,  adapted  from  Ref.  16,  shows  experimental 
values  of  the  lift  curve  slope,  m,  (per  degree)  of  NACA  0015  in  the  presence  of  walls 
having  OAR  =0,  60,  70,  and  80%.  The  results  show  a  convergence  toward  a  free  air  (c/h=0) 
lift  curve  slope  value  of  0.093/deg  and  indicate  small  corrections  for  all  airfoil  sizes 
in  the  range  of  OAR  between  60%  and  70%,  in  accordance  with  potential  flow  predictions. 

This  encouraging  finding  together  with  other  experimental  verifications  in  Refs.  16 
and  17  is  an . indication  that  the  described  concept  of  the  low-correction  or  "tolerant" 
test  section  could  be  a  viable  alternative  to  the  far  more  complex  concept  of  a  "self- 
correcting"  wind  tunnel. 
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(Courtesy  of  F.Mavriplis , Ref . 2) 


Fig.l  Eight-element  cascade  system  simulating 
a  high  lift  airfoil  between  wind  tunnel  walls 
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(Courtesy  of  F.Mavriplis, Ref . 2) 


Fig.2  Computed  W.T,  wall  effect  on  pressure  distribution 
and  comparison  with  experimental  data  at  Sf  =  40° 


(Courtesy  of  F.Mavriplis,Re£ .2} 


Fig.3  Computed  W.T.  wall  effect  on  pressure  distribution 
and  comparison  with  experimental  data  at  6f  =  70° 


(courtesy  of  V •  Y •  Chan » Ref* 7 ) 


Fig.4  Normal  velocities  at  the  edge  of  the  boundary  layer 
as  a  function  of  pressure  coefficient 
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Fig.5  Determination  of  boundary  pressures 
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Fig.6  Measured  wall  pressures  for  a  two-component  airfoil  at  low  lift. 
Wall  corrections  evaluated  according  to  Ref.l  1 
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Fig.7  Measured  wall  pressures  for  a  two-component  airfoil  near  maximum  lift. 
Wall  corrections  evaluated  acci  rding  to  Ref.l  1 
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Fig.8  Measured  wall  pressures  for  a  two-component  airfoil  in  stall. 
Wall  corrections  evaluated  according  to  Ref.l  l 
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Fig.9  Corrected  and  uncorrected  lift  curves  for  a  two-component  airfoil 
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Fig.10  Corrected  and  uncorrected  drag  polars  for  a  two-component  airfoil 
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Recent  U.  S,  developments  concerning  correction  technlqies  for  high  angle-of-attack  testing  are  reviewed 
In  general  terms  and  the  results  are  presented  of  a  letter  st rvey  on  the  methods  now  In  use.  The  majority  of 
the  paper  concerns  the  determination  of  corrections  from  wall  pressure  measurements.  The  application  of  the 
"wall  pressure  signature"  technique  Is  demonstrated  in  experiments  on  several  types  of  models.  The  method  is 
shown  to  provide  good  estimates  of  tunnel  blockage  effects  and  extension  to  lift  interference  Is  discussed. 

It  appears  that  correctab? 11 ty  will  be  limited  more  by  the  problem  of  determining  the  effects  of  tunnel- 
Induced  velocity  gradients  than  by  our  ability  to  determine  the  flow  field.  It  is  suggested  that  passive, 
boundary  measurement  technology  will  diffuse  first  Into  high  angle-of-attack  production  testing,  possibly 
followed  by  partially  adaptive  tunnel  techniques. 


1 , 0  Introduction 
1 . 1  Overview 

The  present  meeting  concerns  methods  of  correcting  for  wind  tunnel  constraint  effects  on  models  at 
high  angle  of-attack  throughout  the  subsonic  and  transonic  flight  regime,  including  the  effects  of  powered 
lift  and  directed  thrust.  These  requirements  are  a  substantial  challenge  to  present-day  and  to  evolving 
technologies  so  a  review  of  all  the  relevant  topics  woqjd  be  very  comprehensive.  An  attempt  will  be  made 
In  this  overview  to  at  least  mention  most  of  the  relevant  t »pics,  but  the  scope  of  the  remainder  of  the 
paper  will  be  more  limited. 

New  approaches  to  correcting  tunnel  data  for  bcuidary  Interference  have  recently  become  practical 
as  a  result  of  continuing  improvements  in  data  acquisition  i  nd  in  both  on-  and  off-line  computing.  Flow 
measurements  at  or  near  the  tunnel  boundaries  can  now  be  included  as  a  starting  point  for  tunnel  correction 
procedures,  rather  than  solely  model -geomet ry  as  in  the  past.  This  avoids  two  problems:  the  need  to  esti¬ 
mate  explicitly  the  effective  aerodynamic  shapes  of  models  w  th  separated  flows  and,  in  high  speed  tunnels, 
the  theoret’cal  difficulties  which  arise  for  ventilated  test  section  walls  which  have  differing  inflow  and 
outflow  characteristics. 

Correction  for  tunnel  effects  involves  two  stages  which  are  not  always  distinguishable  in 
traditional  approaches.  These  are  determination  of  tunnel  -  induced  changes  In  flow  field  velocities  and 
determinatirn  of  consequent  changes  In  surface  pressure  or  force.  Neither  .problem  is  easy  but  the  first 
can  at  l  ;a  t  be  "self-contained"  in  the  sense  that  interference  velocities  can  be  found  independently  of 
model  mearu.-ements  if  conditions  at  the  tunnel  boundaries  are  determined. 

Interference  Flows  4 

In  solld-walled  tunnels,  pressure  measurements  at  the  tunnel  surfaces  may  be  used  to  determine 
singularity  strengths  which  define  the  model  envelope,  Including  both  regions  of  open  or  closed  separation 
and  the  far-wake  displacement  surface.  The. degree  of  detail  available  in  the  model  description  is  a  function 
of  the  number  and  disposition  of  the  tunnel  surface  measurements  employed  as  boundary  conditions  to  the 
problem.  It  will  be  seen  in  a  later  section  of  this  paper  that  the  number  of  measurements  actually  required 
is  practical  for  modern  data  acquisition  and  reduction  systems.  Having  determined  the  singularity  strengths 
which  represent  the  model  envelope,  image  methods  or  their  equivalent  may  be  used  to  compute  the  interference 
flow  field. 


In  porpus  or  slotted-wal led  tunnels,  the  procedure  Is  similar  In  principle  but  less  straight¬ 
forward  in  application.  Normal  velocities  are  now  present  at  the  tunnel  boundaries  and  these  must  be 
measured.  To  avoid  local  effects  near  slots  or  holes  In  the  tunnel  surface,  a  control  surface  Is  usually 
defined  a  short  distance  inside  the  physical  boundary.  The  experimental  techniques  are  more  difficult  than 
before  because  the  normal  velocities  are  small.  Pressure  probes  used  to  treasure  pitch  angle  are  quite 
inconvenient  and  may  lack  the  desired  sensitivity.  The  use  of  the  laser  velocimeter  has  been  suggested  for 
flow  measurements  near  the  boundary. 

The  above  methods  are  essentially  extensions  of  traditional  techniques.  The  important 
difference  is  that,  rather  than  relying  on  known  model  geometry  and  "guesstimated"  wake  profiles,  new  more 
reliable  Information,  i.e.  measured  boundary  conditions,  Is  employed  to  size  and  locate  the  singularities 
which  describe  the  effective  model/wake  surface. 


Effect  at  the  Model 

The  fact  that  the  Interference  flow  field  is  calculable  does  not  guarantee  that  Its  effects 
on  the  model  will  be.  For  simple  attached  flow  cases  surface  pressure  correction  may  be  straightforward. 
Measured  forces  or  moments  may  be  also  correctable,  in  moderate  velocity  gradients,  though  with  somewhat 
less  confidence.  However,  if  tunnel -induced  velocity  gradients  are  large  enough  to  modify  the  character  of 
the  model's  viscous  flows  substantially,  proper  correction  may  be  difficult  or  Impossible.  Little  guidance 
can  be  derived  from  the  literature  concerning  what  constitutes  an  excessive  interference  gradient  in  this 
sense. 
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Even  if  the  nature  of  the  viscous  effects  Is  proper,  moment  correction  (in  particular)  In 
high  angle-of-attack  tests  may  not  be  trivial.  As  pointed  out  by  Heyson,  each  of  the  three  moments  involved 
In  typical  angle-of-attack  pius  yaw  cases  could  be  affected  by  one  or  more  of  the  three  spatial  gradient!  of 
each  velocity  component.  Though  not  all  of  the  nine  gradient  terms  may  affect  a  particular  moment,  it  Is 
evident  that  careful  thought  and  work  is  needed  to  produce  an  effective  and  reliable  correction  schefne. 

The  traditional  solution  to  the  above  problem  Is  to  use  models  which  are  small  in  relation 
to  tunnel  cross  section  and  accept  the  resulting  small  Reynolds  number.  This  is  the  approach  now  being 
taken  n  testing  models  of  fighters  at  high  angles-of-attack.  However,  power  simulation  requirements  fre¬ 
quently  dictate  the  use  of  relatively  large  models,  particularly  in  V/STOL  tests.  The  continuing  use  of 
large,  existing  facilities  for  these  tests  indicates  that  there  will  be  a  continuing  need  for  correction 
procedu  es  for  tests  with  high  levels  of  interference  - 

Adaptive  and  Partially  Adaptive  Wind  Tunnels 

Adaptive  wind  tunnel  technology  has  been  developing  steadily  during  the  past  decade.  Here, 
the  fir3.t,  wall-sensing  stage  Is  the  same  as  described  above  but  the  measured  wall  data  are  used  to  determine 
surface  contour  modifications  or  wall  porosity  changes  which  eliminate  interference  velocities  within  the 
test  section.  The  technology  is  currently  well  advanced  for  two  dimensional  testing  and  the  principles  for 
extension  to  three  dimensions  are  understood.  However,  the  practical  difficulties  for  three  dimensional 
applicat'on  are  quite  impressive.  Even  when  these  are  overcome,  there  is  a  substantial  body  of  opinion  that 
no  more  than  one  or  two  such  tunnels  will  be  commissioned  in  the  U.  S.  within  the  next  decade.  Nonetheless, 
such  tunnels  will  be  very  valuable  for  providing  "benchmark"  test  results  for  checking  high  angle-of-attack 
results  from  conventional  tunnels. 

A  number  of  compromise  solutions  have  been  suggested  which  lie  between  the  entirely  passive 
and  fully  adaptive  approaches.  Some  of  these  recognize  that  the  principle  problem  concerns  velocity  gradi¬ 
ents,  since  even  gross  tunnel -induced  changes  In  mean  angle-of-attack  or  speed  cause  little  difficulty  in 
applying  corrections.  Kemp  has  suggested  the  interesting  concept  of  the  "correctable"  wind  tunnel  which  is 
adaptive  only  to  the. degree  that  velocity  gradients  are  removed.  Once  the  mean  interference  velocities  are 
determined  in  such  a  tunnel,  completion  of  the  correct!  «r»  process  is  likely  to  be  trivial. 

Another  compromise  solution,  suggested  Y/  Joppa,  is  to  employ  a  three-dimensional  system 
which  is  adaptive  only  in  the  lifting  direction.  While  this  can  remove  longitudinal  variations  in  angle* 
of-attack  and  blockage,  the  impact  upon  transverse  gradients,  which  may  also  be  significant,  is  unknown 
One  of  the  difficulties  with  such  "mixed"  systems  Is  de  ermining  the  residual  corrections.  The  fact  i  h.  t 
the  effective  floor  and  roof  contours  are  irregular  in  shape  and  change  from  data-polnt  to  data-point  na  •. 
limit  the  usefulness  of  this  particular  compromise. 

Research  along  the  traditional  line  of  tailoring  porous  or  slotted  walls  for  minimum  inter 
ference  has  also  continued  in  recent  years.  Variation  o'  porous  wall  open  area  and  the  use  of  "tailored" 
slots  designed  via  recently  developed  methods  are  being  pursued  for  tests  on  limited  classes  of  models. 

This  work,  also,  is  relevant  to  high  angle-of-attack  testing,  particularly  if  it  can  be  combined  with 
passive  wall  sensing  schemes,  since  it  is  applicable  to  the  present-day  tunnels  which  will  be  the  mainstay 
for  production  testing  for  a  considerable  time  to  come. 

1.2  Literature  and  Letter  Surveys 


In  preparation  for  the  present  paper  literature  and  letter  surveys  were  made  on  the  subject  of 
high  angle-of-attack  testing. 

The  literature  survey  involved  both  standard  computer  search  techniques  and  personal  follow-up. 

To  keep  the  references  within  bounds  the  survey  was  restricted  to  U.  S.  work  published  after  1973:  prior 
work  may  be  found  via  AGARD  CP  174  (1975)  or  AGAROograph  109  (1956).  After  eliminating  papers  not  relevant 
to  high  angle-of-attack  testing,  about  fifty  remained  which  are  categorized  in  the  Reference  section  of  this 
paper.  Many  of  the  References  will  not  be  cited  directly;  however,  all  were  reviewed  while  preparing  the 
present  paper. 

A  survey'  letter  was  sent  to  appropriate  personnel  in  nine  aircraft  companies,  seven  government 
agencies  and  four  other  organizations.  Questions  related  to  currently  used  test  and  correction  methods  and 
their  shortcomings  and  to  future  plans.  Opinions  were  solicited  regarding  computer  applications  and 
adaptive-wall  tunnels.  The  replies  showed  an  almost  universal  conversatlsm  in  the  production  testing  com¬ 
munity;  models  are  made  sufficiently  small  that  corrections  may  be  neglected  or  can  be  calculated  conven¬ 
tionally.  In  most  cases  the  techniques  In  use  are  backed  by  large-  versus  small-tunnel  test  comparisons. 
The  exceptions  to  the  "small  model"  philosophy  concerned  fighter  testing  (Reference  22)  and  automobile 
testing  (Reference  8). 


Peitzman,  in  Reference  22,  describes  an  empirical  approach  which  combines  the  standard  "Pope"- 
type  solid-plus-wake  blockage -correction  method  at  low  angles  of  attack  with  a  "continuity"  method  at  high 
angles.  In  the  "continuity"  method,  a  q-value  is  used  which  corresponds  to  uniform  flow  through  a  tunnel 
cross  section  reduced  by  the  projected  frontal  area  of  the  model.  A  linear  combination  of  the  two  methods 
Is  used  between  40-  and  90-degrees.  The  method  was  applied  to  tests  in  a  7*  xlO'  tunnel  involving  two 
models  of  the  YF1 7  having  1.05  ft.  and  2.80  ft.  span  (.03-  and  ,08-scale),  wing  planform  areas  of  ,45%  and 
3.2%,  respectively  of  the  tunnel  cross  sectional  area  and  total  planform  areas  of  0.9?  and  6.4%.  A  blockage 
correction,  to  q,  of  approximately  14%  was  obtained  at  90-degrees  angle-of-attack  for  the  larger  model. 
Though  the  Reference  22  method  worked  well,  it  Is  empirical  and  the  "continuity"  method,  in  particular,  Is 
not  well  founded.  Its  use  would  be  unwise  without  further  correlation  except  for  fighter  configurations 
similar  to  the  YF17  at  small  relative  model  size. 
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In  Reference  8,  Wllsden  describes  the  application  of  the  ''wall  pressure  signature"  method  In 
automobile  testing.  This  method  is  a  developed  version  of  that  given  in  AGARD  CP  174.  The  method  Is 
applied  in  tests  upon  production  cars  and  vans  In  the  Lockheed  l6ix23£‘  wind  tunnel  in  which  model  * 
cross-sections  of  5  to  10%  of  tunnel  area  are  Involved.  Recent  correlations  and  newer  theoretical 
developments  will  be  described  In  Section  2. 


2 . 0  Recent  Developments  In  the  Lockheed  "Wall  Pressure  Signature"  Method 


2.  i  I ntroduct Ion 


Since  Its  first  appearance  In  Reference  24,  the  wall  pressure  signature  method  for  blockage 
estimation  has  been  refined  In  certain  aspects,  various  new  correlations  have  been  performed  and  some 
progress  has  been  made  in  the  development  of  pressure  and  force  correction  procedures  for  nonuniform  Inter* 
ference  flows.  These  developments  are  described  In  detal  in  Reference  26  and  will  be  reviewed  In  sub¬ 
section  2.2.  Recently,  a  new,  non-1 terat 1 ve  approach  to  wall  signature  analysis  has  been  developed  which 
Is  more  rapid  and  more  flexible  than  the  original  method.  The  same  principles  have  been  applied  in  pilot 
studies  regarding  the  determination  of  tunnel-induced  upwash  for  lifting  models.  These  new  developments 
will  be  described  in  subsection  2.3.  As  a  point  of  departure,  the  remainder  of  this  subsection  will 
describe  the  general  approach  and  Its  relationship  to  other  methods. 


Outline,  of  Blockage  Estimation  Method 


In  a  typical  experiment  (Figure  l)  wall  pressures  are  measured  along  the  centerlines  of  the 
tunnel  walls  or  roof.  If  lift  Is  present,  the  mean  of  coi responding  floor  and  roof  pressures  can  be  taken, 
or  sldewsll  pressures  are  used  to  determine  souce  and  sink  strengths,  spans,  and  locations  on  the  tunnel 
centerline  (Figure  2)  which  define  an  effective  body  which  corresponds  to  the  displacement  surface  of  the 
test  modal  and  Its  wake. 


FIGURE  1  DETERMINATION  OF  WALL  PRESSURE  SIGNATURE  FIGURE  2  EFFECTS,  AT  A  WIND  TUNNEL  WALL, 

IN  A  TYPICAL  CASE.  OF  SOLID/BUBBLE  AND  OF 

VISCOUS  WAKE  BLOCKAGE. 


A  typical  wall  supervelocity  signature  for  a  separated  flow  case  (Figure  2,  lower)  Includes 
a  velocity  peak,  just  aft  of  the  model,  which  reflects,  the  presence  of  the  expanding  separation  bubble 
which  closes  further  downstream  leaving  a  viscous  wake.  The  asymptote  downstream  Is  a  result  of  the 
displacement  thickness  of  this  wake. 

For  analysis,  the  measured  profile  is  resolved  Into  symmetric  and  antl-symmetr Ic  parts,  as 
shown  In  the  lower  part  of  Figure  2.  This  simplifies  the  subsequent  steps  since  the  symmetric  part  corre¬ 
sponds  to  solid/bubble  blockage  and  the  ant  1 -symmetric  part  reflects  wake  blockage.  Resolution  of  these 
components  from  the  measured  signature  Is  achieved  via  an  Iterative  procedure. 

The  wake  blockage  parameters,  (^v  and  X2,  emerge  directly  from  the  wake  signature  analysis. 
However,  to  obtain  the  solid  blockage  parameters  Qs  and  c^,  a  chart  look-up  technique  Is  employed  as 
Illustrated  In  Figure  3*  Preparation  of  the  working  charts  is  discussed  In  References  24  and  25.  Inputs 
to  the  charts  are  now  the  symmetric  part  peak  height  (Au/U)max  and  half  width  at  half  peak  height,  Ax/B, 
from  the  resolved  signature,  together  with  the  source  and  sink  span  (bs/B)  which  Is  estimated  from  model 
geometry. 
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FIGURE  3  DETERMINATION  OF  SOLI D/3UBBLE  AND  WAKE  BLOCKAGE  PARAMETERS 
FROM  A  WALL  PRESSURE  SIGNATURE 


In  essence,  the  above  analysis  is  an  engineering  solution  of  a  set  of  algebraic  equations 
involving  five  nonlinear  terms  (source  and  sink  spans  and  locations)  and  two  linear  terms  (source  and  sink 
strengths).  Multiple  solutions  are  possible  and  it  is  necessary  to  demonstrate  that  appropriate  roots  are 
obtained  by  the  above  procedure.  Reference  25  describes  a  Lockheed  study  using  a  nonlinear  solver  algorithm 
which  accepts  sets  of  seven  wall  pressures  as  inputs' and  returns  as  roots  the  seven  aerodynamic  varfrbles. 
The  program  was  provided  with  an  initial  guess  of  tie.  roots  and  on  varying  the  guess  converged  to  differing 
roots.  It  was  established  that  the  source/sink  strengths  were  '‘hard"  variables  (i.e.,  lead  to  closety 
spaced  roots)  and  their  streamwise  locations  were  "soft".  The  table  to  the  right  of  Figure  *4  shows  a  number 
of  sets  of  roots  obtained  for  a  normal-plate  exoenment.  The  left  part  of  the  figure  shows  that  the 
corresponding  distribution  of  interference  velocity  at  the  tunnel  centerline  was  affected  very  little.  The 
chart  method  produces  solutions  very  close  to  thost  from  the  nonlinear  analysis. 

UTfEKFE  HENCE 
VELOCITY  AT 
TJNNEL  [ 
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FIGURE  *4  THE  NATURE  AND  CONSEQUENCES  OF  MULTIPLE  SOLUTIONS  TO  THE 
WALL  PRESSURE  INVERSE  PROBLEM  (NORMAL  PLATE  S/C  = 9.5$) . 

Calculated  Flows 

Interference  effects  may  be  calculated  anywhere  in  the  test  section  by  the  wall  pressure 
signature  method.  Figure  5  shows  Incremental  pressure  coefficients,  at  four  tunnel  cross  sections,  induced 
by  the  image  system  for  a  large  sphere  and  its  wake  at  a  subcrltical  Reynolds  number.  The  corresponding 
supervelocities  on  the  tunnel  centerline  increase  from  7*5%  at  the  front  of  the  sphere  to  9*8%  at  its 
maximum  diameter  and  rise  to  11.3%  and  11.2%  at  planes  C  and  D  due  to  the  expanding  separation  bubble. 
Interference  velocity  decreases  downstream  of  C  as  the  separation  bubble  closes. 

Though  axial  variation  Is  very  marked,  the  lower  part  of  Figure  5  shows  that  the  percentage 
variations  at  any  one  cross  section  are  not  large  —  particularly  in  the  region  occupied  by  the  model,  it 
may  be  Inferred  from  continuity  considerations  that  the  blockage-induced  lateral  velocities  are  also  small. 
This  Is  confirmed  by  calculation. 

It  is  sometimes  instructive  to  calculate  streamlines  for  the  effective  body  corresponding 
to  the  source,  source,  sink  solution.  Figure  6  shows  some  body  shapes  for  a  15“inch  sphere  and  a  9“Inch 
normal  disc  tested  in  a  30-  by  A3**inch  wind  tunnel.  Considering  that  only  three  point  singularities  are 
used,  the  effective  shapes  blend  with  the  model  outlines  very  well.  As  Is  evident  from  Figures  5(a)  and 
5(b),  the  change  In  wake  character  through  sphere  transition  Reynolds  number  is  well  depicted.  It  is 
Interesting  to  note  that  the  subcrltical  wake  occupies  2*4.9%  of  the  tunnel  cross-section,  at  maximum 
diameter,  compared  with  13.7%  for  the  model  Itself.  Generally  similar  results  are  obtained  for  a  circular 
disc  mounted  normally  to  the  flow  [Figure  5(c)].  The  calculated  body  diameter  at  the  model  station 
approximates  the  disc  diameter  very  closely. 


.OCATlON  or  MAXIMUM 
I  WALL  DISTURBANCE 


FIGURE  5  DISTRIBUTION  OF  BLOCKAGE- I NOUCED  PRESSURE  COEFFICIENTS,  ACROSS  THE 
TUNNEL  CROSS-SECTION,  FOR  A  LARGE  SUBCRITICAL  SPHERE  (S/C  =13-7%). 


In  the  subcritical  case  (Figure  5)  •  the  calculated  wake  occupies  a  substantial  proportion 
of  the  tunnel  area  and  It  appears  likely  that  Interference  will  distort  it.  To  estimate  the  magnitude  of 
the  change,  body  shapes  were  calculated  with  and  without  iiterference  velocity  Increments.  It  was  found 
that  the  change  in  shape,  due  to  the  tunnel,  was  almost  too  small  to  draw  separately  (see  Reference  26). 
As  the  model  Is  very  large  compared  with  those  used  for  routine  tests,  it  appears  that  blockage- Induced 
distort! jn  of  separation  bubbles  may  not  be  a  limiting  factor. 


<b)  SPHERE,  £-650,000,  S/C -0.137 


FIGURE  6  EFFECTIVE  BODY  SHAPES  IN  30"  xW"  TUNNEL. 

Comparisons  With  Other  Methods 

Blockage  correction  methods  fall  broadly  into  two  classes,  requiring  geometric  or  geometric 
plus  measured  tunnel  data,  respectively.  In  the  absence  of  pressure  data  from  the  tunnel,  aerodynamic 
estimates  have  to  b^  made  which  are  often  conjectural.  Figure  7  shows  a  comparison  between  results  from 
five  blockage  estimation  methods  applied  to  bluff  models,  using  data  taken  largely  from  Reference  8.  The 
"free  air"  data,  in  the  last  column,  is  derived  either  from  sub-scale  models  (for  the  Idealized  sedan)  or 
correspond  to  accepted  values  for  the  Idealized  shapes  (sphere  and  normal  plate). 


AREA  RATIO 
HETH00 

THOM 

HENSEL 

HASKELL 

WALL 

PRESSURE 

SIGNATURE 

FREE 

AIR 

0ATA 

IDEALIZED  SEDAN 

S/C  -  5.65% 

-2.8% 

-7.4% 

-7.2% 

-9.6% 

-8.9% 

-9.0% 

VAN 

S/C  -10.16% 

-5.1% 

-22.1% 

-20.3% 

-- 

-22.9% 

-- 

SUBCRITICAL  SPHERE 

S/C  -13.7% 

-7.0% 

-17.9% 

-27+  -32% 

-29-3% 

-30.0% 

NORMAL  FLAT  PLATE 

S/C  -16.7% 

-9.0% 

-19.0% 

-27.0% 

-52.0% 

-47-0% 

-*5.5» 

FIGURE  7  BLOCKAGE  CORRECTIONS  TO  0RAG  COEFFICIENT  BY  VARIOUS  METHODS. 
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It  is  evident  that  the  two  geometric  methods  give  lower  blockage  estimates  than  two  of  the 
three  methods  employing  tunnel  measurements.  Hensel's  method,  published  In  1951  (Reference  23)  relies  upon 
the  wall  pressure  measurement  as  does  the  present  wal 1 -pressure  signature  method.  However,  it  has  been 
found  that  use  of  pressures  "opposite  to  the  model  location"  as  recommended  by  Hensel  leads  to  uncertainty, 
since  strong  pressure  gradients  are  usually  present  at  this  location  (see  Figure  A).  Hensel's  method  has 
been  revised  recently,  in  a  pilot  study  at  Lockheed,  to  use  the  entire  pressure  signature.  Another  feature 
of  the  Hensel  method  Is  that  It  uses  doublets  to  model  the  flow.  This  renders  the  method  less  flexible 
than  the  wall  pressure  signature  method,  which  permits  variable  source-sink  spacing.  This  can  be  Important 
for  long  models  or  models  with  thick  wakes.  Corrections  via  the  Hensel  method  are  generally  lower  than  by 
wall  pressure  signature  method  -  noticeably  so  in  the  case  of  the  sphere  in  Figure  7. 

Haskell's  method  determines  dynamic  pressure  correction  from  an  empirical  relationship 
involving  measured  drag  or,  in  a  variant  of  the  method,' measured  drag  In  combination  with  base  pressure. 
Haskell  shows  that  the  drag  coefficient  is  proportional  to  base  pressure,  for  square  plates,  and  notes 
that  non-square  plates  exhibit  varying  base  pressure  across  their  span,  in  the  equation*  Aq/q=qCg  S/C, 
n  is  assigned  tabulated  values  which  vary  from  2.77  to  2.13  for  plate  aspect  ratios  In  the  range  1  to  10. 
Haskell's  method  Is  fundamentally  based  and  gives  good  results  for  normal  flat  plates.  However,  in  other 
applications,  determination  of  base  pressure  for  direct  use  In  the  method  may  be  questionable  not  only 
because  of  aspect  ratio  effects,  but  also  if  vortex  separation  exist  In  the  vicinity  of  the  separation 
bubble.  If  the  tabulated  values  of  n  are  used,  there  may  be  practical  difficulties  in  distinguishing 
bubble  drag  from  other  drag  components  or  in  estimating  the  aspect  ratio  of  the  separation  bubble.  Though 
Haskell  asserts  that  an  n  value  of  2.5  40.25  is  reasonably  accurate  In  situations  where  q-correctlons  are 
of  order  10%,  the  spread  in  estimated  corrections  becomes  unacceptable  for  large  models  (see  Figure  7). 

it  has  been  noted  already  that  the  wall  pressure  signature  method  may  be  regarded  as  an 
extension  or  generalization  of  the  Hensel  method,  it  may  also  be  shown  that  the  present  method,  when  imple¬ 
mented  in  its  original  form  (Ref.  2A)  exactly  parallels  Haskell's  approach  while  avoiding  both  Base  pressure 
measurements  and  estimation  of  q.  The  relationshi p  between  the  two  methods  is  examined  further  In  Reference 
26.  • 


2.2  Correction  of  Pressures  and  Forces 

Having  determined  the  spattal  distribution  of  the  three  components  of  Interference  velocity,  as 
described  above,  the  correction  process  is  completed  by  determining  their  effects  upon  the  test  model.  This 
presents  as  great  a  challenge  as  determining  the  interference  velocities  themselves,  and  our  ability  to  ob¬ 
tain  reliable,  accurate  corrections  depends  upon  the  amount  of  detailed  aerodynamic  data  available  from  the 
model . 


Surface  Pressures 


It  Is  demonstrated  in  Reference  26  that  surface  pressure  measurements  may  be  corrected 
successfully  for  blockage  by  applying  x-dependent  Interference  velocity  ratios  to  local  surface  velocities 
and  then  recalculating  Cp.  it  is  shown  that  it  Is  unnecessary  to  consider  transverse  velocities,  even  for 
large  models,  but  it  is  essential  to  ratio  the  surface  velocities  since  simple  superposition  impl les 
velocity  vectors  which  violate  local  tangency  conditions  at  the  model  surface.  The  corresponding  surface 
pressure  correction  equation  is 


CPCW 


cPuM  ■  1 

(i  + 


Figure  8  shows  pressures  measured  on  the  surface  of  the  15-inch  sphere  in  a  30"xii3"  wind 
tunnel,  designated  Cp..  and  in  a  164'  x234'  wind  tunnel,  designated  Cp^.  Blockage  increases  both  the  peak 
and  the  wake  suctions  signif icantly  in  the  small  tunnel.  Correction  “via  the  above  equation  leads  to  good 
correlations  except  that  large  tunnei  peak  suctions  are  slightly  higher. 


FIGURE  8  APPLICATION  OF  BLOCKAGE  CORRECTIONS  TO  SPHERE  SURFACE  PRESSURE  MEASUREMENTS. 
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To  avoid  confusion  with  blockage  velocity  ratio ,  Maskell's  e  has  been  replaced  by  n. 


Similar  tests  encompassing  the  transl t io;ial  Reynolds  number,  reported  In  Reference  2b,  Lead 
to  the  conclusion  that  use  of  the  above  correction  equation  gives  corrected  small  tunnel  data  which  are 
comparable  with  and  no  less  credible  than  results  from  similar  tests  In  the  large  tunnel.  In  correcting 
Reynolds  number,  the  Interference  velocity  ratio  calculated  at  the  peak  suction  location  was  employed.  Modal 
center  v? lues  were  a  reasonable  approximation  to  this. 

Force  Corrections 

it  is  possible  In  principle  to  determine  force  corrections  from  the  pressure  corrections 
Indicated  above.  However,  this  Is  frequently  Impractical  and  alternative  means  are  needed  which  are  less 
detailed  but  which  reflect  interference  gradient  effects  properly.  To  correct  spatially-distributed  loads 
for  spatially-distributed  Interference  effects,  two  general  approaches  may  be  considered,  based  respectively 
upon  test  model  properties  or  measurements  and  upon  fundai  entai  results  for  a  theoretical  flow  model  derived 
from  wall  pressures  and  other  data. 

If  a  good  analytic  counterpart  to  the  test  configuration  exists,  this  may  be  used  In  dividing 
measured  total  loads  between  the  model's  components,  prior  to  correcting  them.  Measured  pressures  may  also 
be  used  In  this  procedure.  In  the  absence  of  these  options,  It  may  be  necessary  to  instrument  a  model  to 
measure  those  component  loads  which  cannot  be  estimated  confidently,  (in  multi-tunnel  tests  on  a  NASA  fan- 
In-wlng  model,  for  example,  a  metric  tailplane  was  used  and  a  particular  division  of  loads  between  the  wing 
surface  and  the  embedded  fans  was  assumed.)  Next,  the  sensitivity  of  each  component  to  changes  in  dynamic 
pressure  or  angle-of-attack  must  be  estimated.  This  will  generally  be  straightforward  for  dynamic  pressure 
changes  but  will  Involve  uncertainties  with  regard  to  angle  of  attack.  Finally,  the  various  effects  are 
summed  and  used  to  correct  the  measured  loads. 

The  source,  source,  sink  blockage  flow  model,  or  Its  lifting  counterpart,  may  be  used  in 
principle  not  only  to  estimate  the  interference  flow  field,  but  also  for  estimating  the  corresponding  forci 
and  moment  changes.  So  far,  this  general  approach  has  b> en  applied  only  to  the  evaluation  of  “horizontal 
buoyancy11  drag.  This  analysis  will  be  used  to  Illustrate  the  principles  Involved. 


FIGURE  9  FLOW  MDDEL  FOR  BLDCKAGE  ESTIMATION 


As' was  seen  previously,  each  source,  source,  sink  combination  determined  from  wall  pressures 
may  be  used  *to  generate  an  effective  body  with  a  defined  area  distribution  and  with  defined  pressures  on  its 
surface.  In  principle,  the  body  shape,  Its  free-air  surface  pressures  and  the  tunnel -Induced  Increments 
could  be  calculated,  the  pressure  correction  could  be  applied,  and  appropriate  Integrations  could  be  per¬ 
formed  to  determine  the  corresponding  drag  Increments.  However,  the  same  result  is  obtained  much  more 
easily  by  considering  the  drag  of  the  sources  and  sinks  themselves.  Using  the  notation  shown  in  Figure  9, 
this  gives,  for  the  forward  body 

adMDDEL  “  p(~Aus  "  Auw  Qw  +  Aus  Qs^ 


Recognizing  that  the  whole-body,  potent ial -flow  drag  on  the  body  sketched  in  Figure  9  is  zero,  we  may  write 
alternatively 

adM0DEL  “  ‘P  uw  0*1 


Substituting  u“  *  i  (l^/BH)  and  nondl mens  Iona] [zing  yields 


ACp 


BH  r  Ow  l2 

S  Lu„  BHj 


where  BH  represents  tunnel  area, 
pressure  correction  to  give 


The  above  horizontal  buoyancy  correction  is  combined  with  the  dynamic 

CDU  + ACp 

Cn_  c  - 

0  +  e^)2 


where  is  the  value  of  Au/Uw  at  the  model  location. 
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Application  of  the  above  equation  to  the  flat  plate  data  analyzed  In  Reference  2k  gave  very 
similar  but  slightly  smaller  corrections  compared  with  those  calculated  by  the  older  method.  The  normal 
flat  plate  drag  coefficient  which  resulted  was  1.20,  slightly  greater  than  Hoerner's  quoted  value.  Figure 
10  shows  new  data,  for  a  family  of  floor-mounted  rectangular  flat  plates.  For  \ik,  the  largest  plate,  the 
unccrrected  drag  coefficient  is  almost  double  the  corrected  value  and  the  deviation  of  this  point  from  the 
earlier  trend  is  probably  significant.  However,  it  ippears  that  the  above  method  can  be  applied  success¬ 
fully  to  models  having  Cq0  S/C  values  of  at  least  0.12  and  highly  separated  flows. 


;/c 

FIGURE  10  ORAG  OF  RECTANGULAR  PLATES  IN  THE  30"  xk3"  TUNNEL. 


The  same  rectangular  models  were  tested  as  flat-plate  wings  which,  being  mounted  from  the 
floor  as  shown  In  Figure  10,  had  a  whole  wing  aspect  ratio  of  3.  Angl e-of-attack  corrections,  determined 
via  standard  methods,  were  generally  less  than  one  degree.  The  results  are  shown  In  Figure  11. 

Oespite  the  very  large  size  of  wing  Wk  (circles)  and  despite  the  presence  of  an  uncontrolled 
boundary  layer  on  the  tunnel  floor,  the  I  i  ft  .'correlation  between  the  four  models  is  excellent  over  the 
lineal  -  a  range.  There  was  a  tendency  for 'smaller  models  to  stall  earlier,  which  was  probably  due  t">  the 
greater  adverse  influence  of  the  floor  boundary  layer.  The  lift  curve  slope  is  about  )0%  greater  than  that 
given  by  simple  wing  theory;  probably  on  account  of  lift  from  edge  vortices  at  the  wing  tip.  The  I ,r  est 
wing  gave  the  best  overall  result  and  was  the  least  affected  by  the  uncontrolled  floor  boundary  laye:*  .'hich 
degraded  the  drag  correlation  somewhat  in  the  linear-CL  range  (see  Reference  26). 

Figure  12  extends  the  lift  and  drag  data  to  90-degrees  angle  of  attack,  i.e.  the  normal  flat 
plate  case.  Uncorrected  data  are  also  shown  for  the  largest  plate.  The  correlation  between  models  remains 
good  except  in  the  immediate  post-stall  region.  Here,  it  appears  that  the  floor  boundary  layer  again  had  an 
adverse  effect  on  the  smaller  models,  this  time  In  delaying  the  approach  to  the  second  lift  peak. 


SYMBOL  WIMG  5/C 

O  W*  16.7* 

+  Vi  9.8* 

X  V2  6.3* 


FIGURE  II  LIFT  CORRELATION  FOR  FOUR  ,  FIGURE  12  DRAG  POLAR  FOR  AR-3  FLAT  PLATE  WINGS 

AR-3  FLAT  PLATE  WINGS.  AT  HIGH  ANGLES  OF  ATTACK. 
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The  above  results  show  that  the  combined  blockage  and  conventional  angle-of-attack  corrections 
produce  good  results  for  large  wings.  Including  one  for  which  Cl  $/C  was  0.13***  In  the  presence  of  a  floor 
boundary  layer,  the  largest  semi-span  model  gave  the  best  results.  Similar  force  studies  on  spheres  yielded 
less  conclusive  results  due  to  hysteresis  and  other  test  difficulties  described  In  Reference  26.  However,  It 
was  concluded  that  the  corrected  small  tunnel  drag  data  were  no  less  credible  than  those  measured  In  the 
large  tunnel. 

2.3  Reformulation  and  Hatrlx  Solution 

The  wall  pressure  signature  method,  described  above,  Is  now  being  used  at  the  Lockheed-Georgia 
l6£x23f-foot  wind  tunnel,  and  Is  also  being  Implement -d  at  other  tunnels  In  the  U.S.  In  the  period  from 
1975  to  1978  the  method  was  refined  as  described  above  and  the  run  time  per  data  point  was  reduced  from 
about  30  seconds  to  3  seconds  on  a  late  '60's  vintage  computer  (CDC  1700),  Nonetheless  the  computation  time 
was  considered  excessive  for  on-line  applications  in  a  ‘production1  wind  tunnel. 

In  setting  up  the  source,  source,  sink  blockage  model  for  the  present  wall  pressure  signature 
program,  a  decision  was  made  to  accept  the  difficulty  of  determining  unknown  source  positions  in  exchange 
for  better  physical  insight  and  the  ability  to  implement  the  procedure  in  chart  form.  However  the  need  for 
faster  operation  and  for  an  additional,  angle-of-attack  correction  capability  has  lead  to  a  study  of  a  dif¬ 
ferent  method  which  employs  multiple  singularities  at  fixed  positions.  This  approach  Is  inherently  faster 
because  wall  influence  matrices,  their  inverses  and  centeiline  interference  matrices  can  be  precalculated 
and  used  non-i teratl vely.  In  view  of  previous  experience  the  revised  techniques  have  been  applied  first  to 
the  blockage  problem  (see  below).  Significant  improvements  In  run  time  have  been  achieved,  but  the  new 
method  Is  neither  well  enough  understood,  nor  correlated  sufficiently  to  be  considered  yet  as  a  substitute 
for  the  present  approach. 

A  welcome  consequence  of  the  matrix  approach,  when  applied  to  angle-of-attack  correction  Is  that 
an  estimate  Is  obtained  of  the  axial  distribution  of  loads'.  This  may  provide  at  least  a  partial  solution 
to  the  problem  mentioned  In  Subsection  2.2,  of  loads  assignment  during  blockage  correction.  It  remains  to 
be  seen  'iow  this  will  work  out  In  practice. 


Influence  Matrices  for  Source  and  Vortex  Arrays 

Figure  13  shows  example  Influence  coefficient  matrices  for  line  source  and  horseshoe  vortex 
arrays  for  a  rectangular  tunnel  having  fl  width-to-he I ght  ratio.  All  of  the  coefficients  correspond  to 
axial  velocity  increments  except  the  influence  of  the  horseshoe  vortex  system  at  the  sidewall  surface, 
which  Is  of  course  a  vertical  velocity,  as  indicated  In  the  table.  Source  coefficients  are  normalized  upon 
mainstream  velocity  and  tunnel  area;  vortex  coefficients  are  normalized  upon  mainstream  velocity  and  the 
square  root  of  tunnel  area.  At  a  given  point,  the  roof  and  floor  coefficients  have  like  signs  for  sources 
but  opposite  signs  for  the  vortex  array:  angle-of-attack  and  blockage  effects  can  therefore  be  separated 
readily  by  considering  the  sum  and  difference  velocities  from  corresponding  locations  on  the  tunnel  roof 
and  floor. 

There  Is  no  theoretical  need  to  use  sidewall  pressures.  In  fact,  we  shall  see  that  there  are 
reasons  not  to.  In  practice  and  particularly  in  large  tunnels,  the  use  of  floor  Instrumentation  Is  unde¬ 
sirable  on  account  of  susceptibility  to  damage  and  complications  related  to  turntables  and  model  supports. 
From  this  viewpoint  the  combined  use  of  sidewall  and  roof  static  pressures  Is  preferable.  Measurement  of 
sidewall  vertical  velocities  is  undesirable  because  of  the  potential  for  damage  to  pitch  probes,  and  because 
of  possible  errors  arising  from  trailing  vortex  drift. 
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FIGURE  14  TUNNEL  SURFACE  PRESSURES  FOR  A  HORSESHOE  VORTEX,  (Cl/AR  =  1 . 0) 

Figure  14  shows  an  example  of  the  distribution  of  surface  pressures  Induced  by  a  horseshoe 
vortex.  This  vortex  spans  50%  of  the  tunnel  width  and  las  strength  corresponding  to  unit  (Ci/At)  •  The 
corresponding  Cl^l  value  of  /2  Is  well  within  the  limit  for  tunnel  flow  breakdown.  At  mid-sidewall  the 
vertical  velocity  opposite  to  the  model  Is  approximately  12%  of  mainstream  (see  inset).  This  causes  a 
sidewall  suction  Cp  of  about  0.015*  However,  far  downstream  the  suction  Is  four  times  this.  This  would 
have  .he  effect  of  reducing  slightly  the  vortex  strength  returned  by  a  Au-dependent  procedure  and  increas¬ 
ing  the  apparent  wake  source  strength.  Because  of  the  nature  of  the  Cp  equation,  correction  for  this 
effect  must  be  retrospective  and  a  second  pass  through  the  analysis  may  be  needed  when  wall  signatures  a'e 
used  In  high  1  if t  tests. 

Application  of  the  Matrix  Method 

In  essence,  the  source  Influence  matrice.  (Figure  13)  comprise  a  0.5000  term  at  every 
location,  associated  with  a  downstream  sink,  superposed  upon  an  antisymmetric  array  representing  the  effects 
of  the  local  sources.  Both  constituent  matrices  are  singular  but  their  sum  is  not.  The  left  side  of  Figure 
15  shows,  a  tunnel  roof  Influence  matrix  (upper  array),  its  inverse  (middle)  and  the  inverse  times  a  center- 
tunnel  influence  matrix  (lower  array).  The  middle  array  is  used  to  determine  source  strengths  from  an 
experimental ,  roof  Au,  blockage  signature.  The  combined  array  in  the  lower  block  is  used  to  determine 
centerline  interference  directly  from  the  roof  signature. 

It  Is  evident  from  Figure  15  that  the  source  Inverse  matrix  has  alternating  signs,  Is  poorly 
dampen  and  has  large  elements  below  the  leading  diagonal.  In  application  to  a  ‘noisy1  experimental  signa¬ 
ture,  v1  strongly-varying,  alternating-sign  source  array  can  be  anticipated.  The  upper  part  of  Figure  16 
confirms  this.  For  this  example  three  wall  signatures  v/ere  computed  for  source,  source,  sink  arrays  at 
0.1B  spacing,  using  source  and  sink  spans  b/B  of  0.30,  0.50  and  0.70.  When  these  were  input  to  an  Influence 
matrix  corresponding  to  a  0.50-span  case  the  source  arrays  depicted  in  the  upper  part  of  Figure  16  were  ob¬ 
tained.  Checks  showed  that,  in  every  case,  the  input  wajl  signature  was  properly  matched  by  the  solution. 
For  the  0.5_span  case  the  original  source  strengths  were  properly  returned. 
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FIGURE  15  APPLICATION  OF  INFLUENCE  MATRICES. 
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The  lower  part  of  Figure  16  shows  that,  despite  the  large  errors  In  source  and  sink  strengths, 
the  distributions  at  the  centerline  deviated  surprisingly  little  from  the  correct,  0.5_span,  signature.  *.hls 
Is  encouraging  and  Indicates  that  even  quite  major  errors  in  source  span  estimation  may  not  have  serious 
repercussions.  Figure  17  shows  the  results  of  a  similar  Exercise,  with  the  source  array  mlsphased  relative 
to  the  source  which  generated  the  "exper Imenta 1  "  signature.  Mlsphaslng  had  almost  no  effect  on  calculated 
lnterfp*  ence. 


To  explore  the  ability  of  the  new  method  to  analyze  complex  signatures,  an  experimental 
example  laving  a  double-peaked  wall  pressure  signature  was  selected  (see  Figure  18).  This  was  generated 
experimentally  by  a  highly  lifting,  Jet  flap  model  which  i  aused  tunnel  flow  breakdown  in  the  case  shown. 

The  first  peak  in  the  wall  signature  corresponds  to  model  blockage  and  the  second  Is  related  to  flow 
separation  from  the  tunnel  floor  aft  of  the  model.  Though  the  second  peak  should  be  represented  by  sources 
or  sinks  located  at  floor  level,  these  were  placed  at  the  tunnel  centerline  for  the  present  1 1 lustratlon- 
The  number  of  centerline  singularities  equalled  the  numbe**  of  experimental  points  (circles,  upper  plot) 
and  the  matrix  solution  passed  through  every  poi.nt.  The  corresponding  centerline  Interference  signature  is 
illustrated  in  the  lower  plot  (circles).  To  determine  the  impact  at  the  model  position  of  the  second  peak, 
a  further  case  was  run  with  It  removed  (triangles).  The  effect  at  the  model  location  (x/B=0)  Is 
surprisingly  smal 1 . 


As  mentioned  previously,  the  matrix  method  is  sensitive  to  data  scatter. 

In  the  Figure  18  case  produced  no  problems  It  has  become  apparent  that  some  smoothing  Is 
has  been  achieved  within  the  matrix  solution  Itself.  An  influence  matrix  Is  set  up  with 
point  but  with  a  reduced  number  of  columns  (i.e.  singularities).  The  matrix  Is  made  squa 
squares  analysis  in  which  all  data  Is  used.  Figure  19  shows  a  least  squares  solution  to 
trated  In  Figure  18.  Nine  singularities,  rather  than  the  previous  nineteen  were  used  and 
approximates  the  wall  pressure  signature  quite' well.  There  is  some  change  to  the  center  I 
signature  In  Figure  19,  as  compared  with  Figure  18;  however,  we  have  no  standard  by  which 
"correct". 


Though  the  scalar 
desirable  and  this 
a  row  for  each  data 
re  by  using  a  least 
the  probl em  i  1  lus- 
the  smoothed  curve 
ine  interference 
to  Judge  which  Is 


OKATCN  WITN  D0U6LE-PEAKE0 
PRESSURE  SIGNATURE 
A  HATCH  WITN  SIHPLIflEO 
PRESSURE  SIGNATURE 


CORRESPONDING  CENTERLINE  SIGNATURES 


CENTERLINE  *  1 
SIGNATURE  .  A 
.A***6 


O  EXPERIMENTAL  OATA 
— 6-H.EAST-S0UARES  OATA 


-o.a  -0.1*  0  0.*  O.fl  1.2  1.6  . 

X/B 

FIGURE  18  APPLICATION  OF  MATRIX  METHOD  TO  A 
COMPLEX  WALL  PRESSURE  SIGNATURE. 


FIGURE  19  APPLICATION  OF  LEAST-SQUARES  METHOD. 
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Application  to  Angle-of- Attack  Correction 

The  tunnel  roof  Influence  matrix  for  the  horseshoe  vortex  array  (see  Figure  13)  Is  symmetric 
about  the  leading  diagonal,  rather  than  anti -symmetric  as  for  sources.  The  Inverse  retains  the  alternating 
sign  property,  but  the  leading  diagonal  Is  now  dominant.  In  view  of  this,  the  solution  may  be  better 
behaved  than  for  sources. 

Pilot,  angle-of-attack  correction  runs  were  made  following  the  same  general  approach  as  for 
the  source  study  of  Figure  16.  The  overall  results  gave  angle-of-attack  corrections  which  agreed  well  with 
conventional  estimates  despite  oscillating  vortex  strength  solutions.  However,  the  overall  sensitivity  to 
span  effects  was  greater  for  the  vortex  model  than  for  sources.  Figure  20,  upper,  shows  three  tunnel  roof 
pressure  signatures  for  Idealized  models  having  the  s'me  lift  (C^  (S/C)=0.20),  but  differing  spans.  It 
was  assumed  that.  In  the  hypothetical  experimental  case,  the  span  was  unknown  and  a  0.508  span  was  assumed 
throughout  the  analysis.  The  vortex  strengths  returned  >y  applying  the  0.5B  matrix  to  each  signature 
deviated  from  each  other  much  less  than  for  the  corresponding  source-source-sink  study  (Figure  16). 

Despite  this,  the  percentage  errors  In  centerline  Interference  were  somewhat  greater. 


FIGURE  20  ANGLE-OF-ATTACK  SOLUTION:  SENSITIVITY  TO  SPAN  EFFECTS 
[CONSTANT  MATRIX  FOR  b=  0.5B.  CL  (S/C)  =0.20] 

Bearing  In  mind  that  the  angle-of-attac'  analysis  of  Figure  20  and  the  blockage  analysis  of 
Figure  16  were  both  made  at  the  same  force  level,  it  ma/  be  shown  that  the  correct  selection  of  span  's 
consideiably  more  Important  for  vortices  than  it  Is  for  sources  at  likely  values  of  model  L/D.  However,  • 
the  sensitivity  Is  such  that  span  selection  within  ±0.10'1  of  the  actual  value  will  hold  errors  to  an 
acceptable  level.  The  possibility  also  exists  of  determining  vortex  span  from  tunnel  surface  mea su remen tr . 
Figure  21  (a)  shows  that  there  Is  little  prospect  for  discriminating  between  vortex  strength  and  vortex  sp  in 
on  the  basis  of  the  shapes  of  streamwise  ceiling  pressure  signatures.  However,  it  Is  evident  from  Figure 
21(b)  that  there  are  good  prospects  for  deriving  vortex  ^pan  from  transverse  pressure  profiles  at  the  model 
locatlor. . 
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FIGURE  21  EFFECT  OF  B0UN0  VORTEX  SPAN  ON  ROOF  SUPERVELOCITIES 
NEAR  M00EL  LOCATION. 


3.0  Discuss  ion 

The  demands  of  high  angle-of-attack  testing  tax  existing  wind  tunnel  corrections  technology  across  a 
broad  front.  However,  after  many  years  of  relatively  slow  progress,  new  advances  are  being  made  which 
promise  to  bring  production  testing  out  of  the  "horse-and-buggy"  stage  and  change  the  present,  small -model, 
small  corrections  philosophy.  This  has  especial  relevance  today  because  of  the  need  for  energy  conserva¬ 
tion  and  because  of  the  increasing  cost  of  wind  tunnel  construction. 

The  new  advances  are  made  possible  predominantly  by  recent  increases  in  computer  power  and  reductions 
In  computation  costs.  The  advances  have  also  fostered  the  growing  science  of  computer  fluid  dynamics  and 
have  raised  the  feasib 1 1 1 ty  of  flow  recomputation  as  a  "correction"  method.  Though  this  is  becoming  a 
reality  for  airfoil  section  design  and  testing  it  will  be  a  considerable  time  before  a  similar  attack  will 
be  practical  for  highly  three-dimensional  flows  such  as  those  over  a i rcraft,  model s  at  high  angles-of-attack. 
Correction  of  test  data,  i.e.  the  semi -empl r leal  removal  of  unwanted  tunnel  effects,  will  therefore  continue 
to  hold  prime  Importance. 

Increased  model  size  brings  with  It  the  ability  to  Improve  Internal  Instrumentation,  to  Install  power 
simulators  more  readily  and  to  improve  reading  accuracy  as  well  as  to  increase  Reynolds  number.  Before  any 
of  these  desirable  features  can  be  exploited,  however,  the  test  engineer  must  ask  "Will  the  data  be  correc¬ 
tive?"  At  all  levels  of  testing,  the  question  resolves  Into  two  aspects,  concerning  the  flow  physics  and 
the  data  needed  to  compute  the  correction. 
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It  was  shown  In  Section  2  that  the  prospects  are  good  for  defining  the  tunnel** Induced  flow  field  from  r 
a  practical  number  of  measurements  at  or  near  the  test-section  boundaries.  However,  an  exception  to  this 
occurs  when  a  sonic  region  occupies  a  substantial  proportion  of  the  tunnel  cross  section  or  approaches  a 
tunnel  surface.  The  prospects  are  also  good  for  deriving  corresponding  corrections  successfully  In  pre- 
and  post-stall  situations,  provided  the  loads  on  major  components  can  be  measured  or  estimated.  However, 
Incipient  stall  and  immediate  post-stall  flows  are  likely  to  be  gradient-sensitive  and  this  may  limit  model 
size  In  f l,  ed-geometry  tunnels.  There  Is  a  clear  need  for  sensitivity  studies  which  determine  the  permis¬ 
sible  Imposed  axial  gradients  for  typical  airfoil  sections  or  wings. 

Tunnel-related  wake  distortion  may  limit  correctabl 1 i ty  In  V/STOL  tests,  especially  If  powered  flows 
Impinge  upon  tunnel  surfaces,  This  Is  a  pressing  problem  in  large  tunnels  where  full-scale  vehicles  are 
tested.  Despite  a  number  of  attempts  to  use  adaptive  tcc  irilques  In  this  situation,  no  practical  solution 
has  emerged.  It  is  not  yet  clear  to  what  extent  the  difficulty  can  be  relieved  using  the  new  correction 
techniques  described  In  Section  2.3  (see  Figure  l 8) . 

The  second  aspect  of  correctabl 1 1 ty  -  concerning  data  needed  to  compute  corrections  -  divides  naturally 
Into  tunnel  surface-  and  model -rel ated  parts.  Tunnel  surface  instrumentation  for  solid-wailed  tunnels 
represents  a  fairly  modest,  one-time  cost  Involving  probably  two  scanivalves.  The  requirement  Is  somewhat 
greater  foi  porous-walled  tunnels  and  better  methods  are  needed  for  measuring  normal -to-wal 1  velocities. 

Part  of  the  data  needed  to  calculate  the  correction  to  a  particular  model  measurement  Is  often  the  measure¬ 
ment  Itself.  This  Is  all  that  Is  required  from  the  model  in  pressure  distribution  tests  and  In  some  force 
tests.  However  extra  data  are  needed  to  correct  six-component  balance  measurements  and  means  must  be 
available  for  distributing  loads  between  components  in  order  to  compute  gradient-dependent  corrections, 
particularly  to  moments. 

It  Is  possible  to  eliminate  one  of  the  weaker  links  In  the  above  procedure  -  that  of  assigning  loads 
to  individual  model  components-  by  using  adaptive  walls.  It  is  of  Interest  to  consider  a  minimally  adap¬ 
tive  tunne*  In  which  changes  are  made  only  to  "upper"  and  "lower"  tunnel  surfaces  and  only  as  needed  to 
remove  axial  gradients.  Such  an  approach  would  reduce  the  recurring  cost  of  model  Instrumentation.  How¬ 
ever,  the  magnitude  of  the  task  of  providing  partially,  as  opposed  to  fully  adaptive  walls,  has  not  been 
defined  ani  cost  effectiveness  studies  are  needed  which  consider  such  trades.  Possible  difficulties  In 
defining  the  mean  flow  and  determining  residual  corrections  for  partially  adaptive  schemes  should  also  be 
reviewed. 

A  significant  factor  In  promulgating  emerging  techno’ogies  concerns  acceptance  by  the  test  engineer. 

When  the  pressure-signature,  blockage-correction  method  wts  first  developed  (see  Section  2),  this  was  a 
major  reason  for  adopting  the  "chart"  method  (Section  2.2)  rather  than  the  more  direct  but  less  Intelllg'b  e 
matrix  approach  (Section  2.3).  Conversations  with  the  testing  community  concerning  fully  adaptive  tunnels 
have  revealed  deeply  Ingrained  reservations  concerning  the  "calibration"  of  these  tunnels.  A  partially 
adaptive  tunnel  used  with  a  baseline  calibration  would  be  more  acceptable  to  the  test  engineer. 

4.0  Conclusions 

Recent  U.  S.  developments  wind  tunnel  test  and  corrections  technology  have  been  reviewed  which  are 
applicable  to  high  angle-of-attack  testing.  The  needs  are  f jeh  that  the  applicable  techniques  encompass  a 
large  proportion  of  current  research  on  tunnel  test  technology. 

An  overview  of  the  relevant  factors  Is  given  in  Section  1,  together  with  the  results  of  a  limited 
literature  review  and  of  a  letter  survey  involving  members  of  the  tunnel  testing  community  In  the  United 
States.  Developments  In  the  use  of  measurements  at  tunnel  surfaces  for  determining  corrections  are 
reviewed.  These  Include  both  improvements  to  the  iterative  approach  published  by  Hackett  and  Wilsden  in 
1975  (see  Sections  2.1  and  2.2)  and  some  early  results  from  a  new,  non-iterative  technique  (Section  2.3). 

The  definition  of  the  tunnel  -  induced  flow  field,  using  wall  measurements,  has  been  accomplished  for 
low  speed  flows  and  extension  to  high  subsonic  speeds  appears  feasible  provided  that  supercritical  flow  is 
not  extensive.  Experiments  have  shown  that  surface  pressures  may  be  corrected  successfully  even  on  quite 
a  large  model . 

The  most  difficult  aspect  of  the  correction  process  is  to  determine  the  consequences  of  tunnel  -  induced 
velocity  gradients.  The  effects  on  moments,  in  particular,  may  require  either  experimental  or  theoretical 
input  in  order  to  assign  loads  to  particular  model  components.  However  new  work  is  needed  to  determine  the 
sensitivity  of  viscous  effects  -  particularly  the  stall-  to  tunnel -induced  gradients. 

The  above  topics  are  reviewed  In  Section  3,  which  also  includes  a  discussion  of  adaptive-wall  techniques. 
A  possible  compromise  Is  to  adapt  the  walls  only  to  a  zero  gradient  condition,  rather  than  to  a  configuration 
for  zero  Interference.  This  might  involve  only  upper  and  lower  tunnel  surfaces  and  might  prove  more  accep¬ 
table  to  the  practicing  test  engineer  than  a  fully  adaptive  scheme. 

It  appears  that  passive,  boundary  measurement  technology  may  diffuse  first  Into  high  angle-of-attack 
production  testing,  possibly  followed  by  some  form  of  partially  adaptive  tunnel  technology.  However,  much 
remains  to  be  done  concerning  the  application  of  adaptive  techniques  to  configurations  with  highly  three- 
dimensional  flows. 
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AMELIORATIONS  ENVISAGEES  POUR  RESOUDRE  LES  PROBLEHES  RENCONTRES 
AU  COURS  D'ESSAIS  A  GRANDE  INCIDENCE  DE  MAQUETTES  EN  SOUFFLERIE 
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RESUME 

Parmi  les  problemes  rencontres  lors  de  1' execution  d'essais  en  soufflerie  de  ma- 
quettes  5  grande  incidence,  trois  sont  selectionnes  ici  :  effets  de  parois^ interaction 
de  dards  supports  et  vibrations  au-dela  du  decrochage, 

L'&tat  de  l'art  des  corrections  des  effets  de  parois,  sys tenia tiquement  appliquees 
en  soufflerie. a  caractere  industriel,  est  expose  avec  de  nombreuses  comparaisons  entre 
aouffleries  ou  entre  vol  et  soufflerie.  Les  modelisations  necessaires  des  maquettes 
an  atmosphere  illimitee  font  apparaitre  des  lac  mes  dans  le  cas  de  tourbillons  d'apex 
ftt  de  jets  actifs  des  maquettes. 

Simultanement,  d'autres  methodes  de  corrections  sont  developpees  en  FRANCE  pour 
pallier  aux  hypotheses  simplif icatrices  et  aux  limitations  de  la  methode  classique. 
Ainsi  AMD-BA  a  developpe  la  methode  des  petits  paves  et  ONERA,  la  methode  des  signa¬ 
tures  maquettes  sur  les  parois  des  veines  d’essais.  Cette  derniere  methode,  dej& 
utilisee  en  essais  bidimensionnels industriels  pour  des  veines  ventilees  cylindriques, 
est  etendu  au  tridimensionnel.  Son  extension  au  cas  des  veines  deformees  constituera 
une  etape  intermgdiaire  de  1' etude  des  parois  autoadaptables. 

II  est  insist^,  en  conclusion,  sur  les  methodes  de  controle  de  la  validity  des 
methodes  de  correction  des  effets  de  parois. 


EXPECTED  IMPROVEMENTS  ON  HIGH  ANGLE  OF  ATTACK  MODEL  TESTING 

Office  National  d'Etudes  et ,de .Reche *ches  Aerospatiales  (ONERA) 
92320  Chatillon  (France) 


SUMMARY 

Several  problems  were  encountered  during  tests  at  high  angle  of  attack  in  wind 
tunnels.  Three  are  selected  here  :  wall  interference,  sting  interference  and  vibra¬ 
tions  beyond  the  stall. 

The  state  of  the  art  on  wall  interference,  systematically  applied  to  the  develop¬ 
ment  tests,  is  showed  with  several  comparisons  between  wind  tunnels  or  between  flight 
and  tunnels  tests.  The  models  used  iri’ unconfined  flow  point  out  some  deficiencies  as 
regards  apex  vortex  and  active  jets. 

Simultaneously,  other  correction  methods  are  developed  in  FRANCE  to  palliate  the 
assumptions  and  limits  of  the  classical  method.  Thus  AMD-BA  works  on  the  vortex  lat¬ 
tice  method  and  ONERA  on  the  model  signatures  on  the  test  section  walls.  This  method, 
already  used  in  2D  flow  for  development  tests  with  ventilated  cylindrical  test  sec¬ 
tions,  is  extended  at  the  3D  case,  as  a  first  stage  for  distorted  but  non  stream¬ 
lined  walls. 

As  a  conclusion,  the  control  of  the  validity  of  the  wall  interference  correction 
method  is  analysed. 
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INTRODUCTION 

Les  provisions  de  performances  d' avion  sont  Atsblies  A  partir  de  norabreuses  donnOes  - 
parmi  lesquelles  les  essais  en  soufflerie  tiennent  un.rdle  majeur.  En  m&me  temps  que  la  fourniturc 
des  rOsultats  de  mesures,  le  spOcialiste  de  soufflerie  doit  effectuer  les  corrections  nOcessaires 
pour  que  le  con3tructeur  puisse  disposer  des  meilleurs  OlOments  d 'apprOciation  de  la  validity  et 
de  la  prOcision  des  rOsultats.  En  dehors  des  comparaisons  entre  rOsultats  d'essais  en  voL  et 
en  soufflerie  qui  nOcessitent  de  recourir  A  des  mAthodes  d 'extrapolation,  les  essais  de  fidAlitA 
et  les  comparaisons  de  rAsultats  obtenus  dsns  diffOrentes  souffleries  ou  diffArentes  configurations 
de  ,7eines  d'esaais  constituent  des  guides  prAcieux  pour  juger  les  arnAliorations  A  apporter  aux 
mAthodes  de  corrections. 

Des  programmes  de  comparsison  de  rAsultats  d'esssis  ont  ainsi  AtA  rAalisAs  en  supersoniqu^ 
sur  maquettes  e talons  AGARD  Cl)  puis  plus  rAcemment  en  transsonique  sur  maquettes  Atalons  0NERA(2) . 
De  telles  confrontations  sont  bien  entendu  spAcifiquss  des  maquettes  utilises  et  en  particulier 
devrsient  Stre  developpAes  dans  le  domsine  des  grsndes  incidences  sur  avion  civil  ou  militaire. 

•  Lors  des  comparaisons  des  rAsultats  provenant  de  divers  essais  effectues  dans  diverses 
veines  d'essais,  il  y  a  lieu  de  tenir  compte  des  qualitAs  de  l’Acoulement  aussi  bien  stationnaire 
(ascendance,  gradient  longitudinal  de  pression)  qu 'imtstionnaire  (turbulence,  bruit),  des 
effete  de  parois  et  des  champs  des  dards  supports  et  -des  conditions  de  traitement  des  mesures 
surtout  lorsque  des  vibrations  interviennent . 

Si  certaines  corrections  peuvent  Atre  effectuAes  ssns  ambiguity  telles  que  celles 
d'sscendance  et  de  poussAe  d'ArchimAde,  d’autres  corrections  telles  que  celles  des  effets  de  parois 
peuvent  poser  des  problAmes  si  les  effets  des  autres  paramAtres  modifient  les  caractAristiques 
sArodynamiques  des  maquettes.  II  y  s  done  lieu  de  s  assurer  de  l'effet  de  ces  paramAtres  tels 
que  turbulence  ou  bruit,  avant  d'aborder  les  corrections  des  effets  de  parois,  sous  peine  de  mAlan- 
ger  diffArentes  influences  ce  qui  conduirait  slors  1  doutcr  des  mAthodes  de  corrections.  Ceci 
risque  d'etre  particuliArement  dAlicat  dans  les  domaines  oA  Is  susceptibility'  des  rAsultats  aux 
par&iAtres  indiquAs  devient  AlevAe  et  lorsque  l'itfiuence  des  divers  paramAtres  ne  peut  etre 
Atudiee  isolement  car  il  n'est  pas  encore  possiblt  \'e- les  fsire  varier. 

L'CNERA  s  dAJ  A  publiA  des  travaux  (3,4)  nur  des  comparaisons  de  rAsultats  dont  cer- 
tsins  concemaient  le  domaine  des  incidences  AlevAes  Les  sujets  qui  seront  ici  abordAs 
conce^tient  particuliArement  le  domaine  des  incidences  AlevAes  ou  plus  prAcisAment  le  domaine 
d 'incidence  pour  lequel  des  dAcollements  importsnta  existent.  Le  sujet  principal  concernera  les 
effets  de  parois  en  presence  de  dAcollements  et  les  f'iverses  mAthodes  de  corrections  en  service 
pour  essais  industriels  ou  envisagyes  dans  le  futur  (5)  .  L'influence  des  dards  supports  sera 

rapidement  abordye  ainsi  qu'une  tentative  de  modyiisation  du  comportement  vibrstoire  des  msquettes 
en  dad  autour  et  au  delA  du  dycrochage. 


I  -  CORRECTIONS  DES  EFFETS  DE  PAROIS 

1. 1  -  Mythode  analytique 

1.1.1  -  Corrections_en  1 'absence_de_dAcollements 

Il  est  seulement  rappeiy  ici  que  Is  mAthode  analytique  de  corrections  nycessite  la 
formulation  d 'hypotheses  de  conditions  aux  limites  sur  les  psrois  de  la  veine  d'essais  et  la 
modyiisation  des  potentiels  maquettes  en  atmosphere  illimitye.  Ces  modyiisations  peuvent  Stre 
plus  ou  moins  yiaboryes  selon  les  tallies  des  maquettes  vis  A  vis  des  dimensions  des  veines 
d'essais  (4)  .  On  reviendra  sur  ce  point  su  parsgrsphe  1.1.5. 

La  mythode  analytique  est  systymatiquement  sppliquye  pour  corriger  les  rysultats 
d'essais  dsns  les  souffleries  A  carsctAre  Indus triel  de  l'ONERA. 

1.1.2  -  Corrections  en  prAsence  de  dAcollements 

Pour  tenir  compte  de  Is  prysence  de  dycollements,  une  correction  suppiymentaire  est 
ajoutAe  aux  corrections  de  blocage  de  volume  et  de  sillsge,  A  1* inverse  de  ces  derniAres, Atabliea 
A  partir  d'une  modyiisation  de  maquette.  Is  correction  de  dycollement  est  dyrivye  de  la  thAorie 
semi-empirique  de  MASKELL  (6,7)  fondee  sur  Is  notion  d'esu  morte,  Les  deux  hypotheses  fon- 
damentales  de  cette  thAorie  sont  rappelAes  seulement  ici. 

La  premiAre  hypothAse  consiste  A  sssirailer  1 'ycouleraent  dycoliy  dans  le  cas  d'une  aile 
A  grande  incidence  au  bulbe  d'eau  morte  derriAre  une  plsque  perpendiculaire  au  vent.  Si  ce  bulbe 
est  bien  dAcelable  dans  les  visualisations  au  tunnel  hydrodynamique  derriAre  une  aile  en  inci¬ 
dence  ou  avec  des  voleta  braquAs,  VAYSSAIRE  (9)  note’ que  dans  le  css  de  tourbillons  d'apex  bien 
organisAs  et  portants,  la  notion  d'eau  morte  n'est  justifiAe  qufen  aval  du  point  d'Aclatement  du 
tourbillon  et  ce  d'autant  plus  que  le  tour billon  est  diffus  A  grande  incidence  et  avec  un  Aclatement 
proche  du  bord  d'attaque. 
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La  deuxidme  hypothdse  de  MASKELL  consiste  en  une  addition  pure  et  simple  du  coefficient 
6t>  aux  deux  autres  coefficients  de  blocage  de5' volume  £ «,  et  de  sillage  £w  (figure  1)  dans 
la  correction  du  nombre  de  Mach  M  ou  de  la  pression  cindtique  cf  .  Cette  superposition 
lindaire  correspondrait  k  une  apparition  de  ddcollement  au  bord  de  fuite  de  voilure  et  k  son 
extension  progressive  k  toute  l'envergure  de  l'aile, 

Le  coefficient  de  blocage  de  ddcollement  (figure  1)  est  proportionnel  k 

1* increment  de  trainee  engendrd  par  les ddcollements .  Ce  supplement  de  trainee  est  deduit 

de  1'dcart  entre  la  polaire  experimentale  et  la  mise  en  equation  polynomiale  de  la  polaire 
da;\s  le  domaine  exempt  de  decollements .  Le  coefficient  de  proportionnalifcd  Q  intervenant  dans 
est  deduit  de  tarage  des  veines  d' essais  k  J'aide  de  plaques  perpendiculaires  au  vent  de 
divers  allongements  Jf8  )  ^q.  .  La  fonction  O  ( Xa.)  indiquee  figure  1  a  ete  obtenue  en  veine 
guidee.  D'aprds  d' autres  tarages  effectues  par  VAYSSAIRE  T9)  en  veine  semi  guidee  de  la 
soufflerie  BREGUET,  une  valeur  voisine  de  2  serait  k  considerer  pour  G  .-En  outre,  d'aprds 
des  tarages  en  bidiraensionnel,  il  n'est  pas  exclu  que  la  valeur  de  O  depende  de  1' incidence. 


1.1.3  -  Applications  preiiminaires 

En' vue  de  valider  la  thdorie  de  MASKELL,  des  comparaisons  de  rdaultats  corrigds 
dana  diverses  souffleries  ont  dte  prdsentdes  en  1972  par  VAYSSAIRE  C9l  Un  bon  accord  dtait 
obttnu  entre  les  rdsultats  corrigds  d'une  maquette  d' avion  de  transport  k  fldche  de  voilure 
de  25  degrds  dans  deux  configurations  de  la  veine  de  la  soufflerie  274  de  St  Cyr  :  guidee 
et  k  parois  horizontales  perfordes.  Pour  un  rapport  entre  l'aire  de  la  section  veine  C  et 
la  surface  de  voilure  S  de  24  les  result at a  coriigds  dtaient  confondus  jusqu'A  Mach  0,85. 


Dans  le  cas  d'une  maquette  de  Mirage  F  k  fldche  de  voilure  de  48  degrds,  en  configura¬ 
tion  hypersustentde,  l'accord  des  courbes  de  por' ance  corrigdes  dana  trois  aouffleries  diffdrente^ 
S5  Toulouse,  S2  Chalais  et  Si  St  Cyr,  pour  des  rapports  C/S  respectivement  de  19.3,  9.6  et  5.1 
(figure  2)  dtait  assez  bon  malgrd  1 'importance  d rs  corrections  pour  le  cas  C/S  =  5.1. 


La  confrontation  de  ces  rdsultats  corrigds  avec  lea  rdsultats  des  essais  en  vol,  en 
configuration  trains  rentrds  ou  sortis,  (figure  3)  dtait  considdrde  comme  largement  satisfai- 
aante  par  le  constructeur  m3me  dans  des  cas  ou  1 ’estimation  du  terrae  £  £>  dtait  grande  par 
rapport  aux  autres  termes  €s  et  S  .. 


1.1.4  -  Comparaisons  de  rdsultats  dans  les  souffleries  SIMA  et  Fl 

Depuis  1 'entree  en  service  en  1977  de  1 p  uouvelle  soufflerie  basse  vitesse  pressurisde 
Fl  de  l'ONERA,  de  nombreuses' comparaisons  entre  l^j  rdsultats  obtenus  dans  cette  soufflerie  et 
dans  Id  soufflerie  Si  de  Modane  ont  dtd  effectuee.  k  m&ne  Reynolds  sur  mdmes  maquettes.  La 
soufflerie  Fl  possdde  une  veine  rectangulaire  guidte  de  15,75  m2  d'aire  ;  la  nouvelle  veine  k 
fentes  de  la  soufflerie  SIMA  k  une  aire  de  40,44  m2  (figure  6).  Les  maquettes  passdes  k  Fl  et 
Si  prdsenteront  done  des  corrections  d'effets  de  parois  notables  k  Fl  et  ndgligeables  k  Si,  de 
telle  sorte  que  les  rdsultats  obtenus  k  Si  peuvent  servir  de  rdfdrence  aux  corrections  appliqudes 
k  Fl. 

.  En  outre,  il  a  dtd  vdrif 14j sur  une  maquette  de  petites  dimensions  du  Mirage  G8,  (rapports 
C/S  de  28,4  et  73,0  k  Fl  et  SI)  que  les  qualitds  des  dcoulements  des  deux  veines  d'essais  condui- 
saient  k  des  rdsultats  parfaitement  identiques  (figure  4).  Un  ddcalage  constant  de  la  stability 
rdsulte  des  differences  de  montage,  en  dard  k  Si  et  sur  mAt  &  Fl.  Dfes  lors,  il  dtait  possible  de 
comparer  des  rdsultats  sur  maquettes  de  tailles  plus  grandes  et  d'dtudier  les  effets  de  parois 
sans  risque  d' influence,  d' autres  paramdtres. 


La  comparaison  des  rdsultats  corrigds  d'une  maquette  au  1/10  du  MERCURE  100  dans  les  deux 
souffleries  montre  une  bonne  concordance  (figure  5)  en  configuration  croiaifere  bien  que  l'envergure 
relative  de  la  maquette  soit  de  0,68  k  Fl.  Les  polaires  dquilibrdes  ob tenues  k  Fl  k  3,9  bars  de 
pression  gdndratrice  sont  compardes  aux  rdsultats  donnds  par  la  moddlisation  faite  par  les  AMD. BA 
et  correspondant  aux  essais  en  vol  (figure  6).  La  concordance  est  trfes  bonne  en  configuration 
croisifere  et  ddcollage.  Les  rdsultats  de  vol  k  1 'atterrissage  correspondent  k  la  position  train 
sorti  alors  que  les  essais  en  soufflerie  ont  dtd  effectuds  train  rentrd  :  on  vdrifie  cependant 
la  bonne  cohdrence  vol-soufflerie  gr3.ee  k  1 'existence  de  rdsultats  de  vol  au  ddcollage  donnant 
l1 Influence  de  la  sortie  du  train  d 'atterrissaga. 


Dana  tous  les  exemples  qui  viennent  d'etre  donnds  jusqu'ici,  les  flfeches  de  bord  d'atta- 
que  sont  telles  que  les  voilures  sont  exemptes  d'dcoulement  tourbillonnaire  d'apex.  La  figure  7 
illustre  les  difficultds  qui  apparaissent  dans  le  cas  de  fortes  flfeches  avec  des  portances 
tourbillonnaires  marqudes.  Les  courbes  correspondant  aux  valeurs  non  corrigdes  et  corrigdes  k 
Fl  permettent  d'apprdcier  1 'importance  des  corrections  :  la  concordance  entre  les  valeurs  corrigdes 
A  Fl  et  SIMA  n'est  pas  enti&rement  satisfaisantc  dans  les  zones  k  forts  ddcollements  au  delA  de 
l'infldchissement  de  la  courbe  de  portance.  Dans  le  cas  de  la  figure  7,  la  correction  de  blocage 
de  ddcollement  est  appliqude  k  partir  de  1* incidence  d 'infldchis semen t  de  la  courbe  C  t_  (  CX  )  ; 

51  (figure  8)  :  les  courbes  corrigdes  sont  intituldes"corrected  1"  sur  planches  7  et  8.  Si  la 
correction  de  blocage  de  ddcollement  avait  dtd  appliqude  k  1* incidence  d'apparition  du  tourbillon, 

52  de  figure  8,  les  courbes  corrigdes  seraient  cellea  intituldea  "corrected  2"  sur  les  planches 
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8  et  9.  La  figure  9  montre  l1  importance  relative  des  coefficients  de  trainee  de  ddcolleraent  Cj>^ 
et  Cj) S2,  tandis  que  la  figure  8  donne  les  corrections  de  nombre  de  Mach  (  -  )  dans  les 

2  cas  de  bornes  de  ddcolleraent.  La  figure  8  montre  bien  qu'aux  incidences  dlevdes  la  correction 
la  pluj  importante  est  celle  de  blocage  sur  la  pressioi  cindtique  (trajets  BC1  ou  BC2)  alors  que 
la  correction  d' incidence  (trajets'AB)  est  moddrde.  Or.  voit  done  apparaitre  dans  le  cas  de  portances 
tourbillonnaires  la  difficult^  de  mise  en  oeuvre  de  la  thdorie  de  MASKELL. 

D'une  part, on  a  vu  que  le  tourbillon  non  dclatd  n'dtait  pas  assimilable  k  une  zone  d'eau 
morte  done  que  la  borne  infdrieure  d1 incidence  de  ddcollement  doit  dtre  plutBt  &  1 'inf ldchis semen t 
de  Ct  (points  S2)  et  pourrait  dtre,  en  basse  vitesse,  bien  ddfinie  par  le  minimum  de  la  courbe 
C A fx)  (figure  9). 

D' autre  part,  la  mise  en  Equation  polynomiale  de  la  polaire  conduisant  k  la  valeur  de  la 
trainde  de  ddcollement  devrait  dtre  issue  de  la  polaire  induite  en  Cl^/mAo. 

distincte  de  la  polaire  expdrimentale  dds  l1 apparition  du  tourbillon  qui  donne  un  supplement  de  Cj>  . 

II  y  a  done  arabiguitd  dans  1 'application  de  la  thdorie  semi-empirique  de  MASKELL  par  suite 
d'une  absence  de  moddlisation  progressive  de  la  trainee  induite  par  le  tourbillon  lorsque  1' inci¬ 
dence  crolt  done  lorsqu'il  se  ddplace  vers  le  plan  de  symetrie  avion  et  delate  en  un  point  progres- 
sant  vers  le  bord  d'attaque. 

Dans  le  cas  d'une  portance  tourbillonnaire  moins  accusde  du  fait  d'un  braquage  de  bees 
par  exemple,  (figure  10)  la  confrontation  est  meilleure  en  Toutefois  la  comparaison  entre 

rdsultats  corrigds  k  FI  et  SIMA  est  moins  bonne  en  ce  qii  concerne  la  stability  et  la  trainde. 

1.2  -  Schdmatisation  des  maquettes 

La  mdthode  analytique,  comme  la  radthode  des  signatures  qui  sera  dvoqude  au  paragraphe 
1.3,  ndeessite  une  schdraatisation  de  l'dcoulement  aucour  de  la  maquette  en  atmosphdre  illimitde. 
S'agirsant  de  calculs  correctifs,  la  moddlisation,  fr'.dde  sur  un  support  theorique  admissible, 
doit  S'rre  raisonnable  et  rester  assez  simple  pour  re  pas  entrainer  de  trop  longues  durdes  de 
calcul.  II  serait  en  effet  inutile  d'apporter  des  .‘a/finements  supplementaires  si  ceux-ci  ne 
conduisrdent  qu'A  de  faibles  variations  des  correction^  .  La  difficultd  principale  reside  dans  la 
ddfinition  d'un  niveau  suffisant  d' dlaboration  des  moddlisations  ;  la  meilleure  mdthode  de  contr81e 
de  ce  niveau  reste  encore  la  vdrification  des  effets  sreondaires  entraines  par  la  prise  en  consi- 
ddratic  i  de  descriptifs  plus  dlabords  de  la  moddlisatiox. 

II  ne  faut  pas  oublier  que  le  support  thdorique  de  la  sch&natisation  de  la  maquette  est 
corapldtde  de  valeurs  expdrimentales ,  globales  telles  qua  C L  t  ou  locales  telles  que  des 

rdpartitions  de  pressions  ou  des  mesures  clinometriques .  En  outre,  des  coefficients  correcteurs 
sont  as.°ocids  k  certains  terraes  de  corrections  pour  tenir  compte  d'expdrimentations  de  base. 

1.2.1  -  Termes  de  blocage  ^ 

La  plupart  du  temps,  les  moddlisations  sont  trdfc  simples,  doublet  et  source,  pour  les 
termes  de  blocage  de  volume  et  de  sillage.  Les  intensitds  de  doublets  sont  ddduites  du  volume 
maquette  et  de  coefficients  correcteurs  tenant  compte  de  l'dlancement  du  fuselage  et  de  l'dpaisseur 
relative  de  voilure.  L'intensitd  des  sources  est  directement  liee  k  la  mesure  de  la  trainde  C. jjj 
maquette  ou  de  la  trainde  de  sillage  de  profil  £j>  . 


L'emploi  de  multisingularitds ,  rdparties  sur  la  maquette,  permet  de  respecter  mieux  la 
gdomdtrie  tie  celle-ci  mais  au  prix  d'un  alourdissement  des  calculs. 


Pour  Jes  termes  de  poussde  d'Archiradde,  induits  par  les  gradients  longitudinaux  de  blocage, 
la  loi  des  aires  maquette  est  incluse  pour  aradliorer  la  correction  rdsultante  de  trainde. 


Les  ddcolleraents  entrainent  une  correction  suppldmentaire  calculde  &  partir  de  la  thdorie 
serai-empirique  de  MASKELL.  Les  difficultds  d'emploi  en  prdsence  de  tourbillons  d'apex,  progressive- 
ment  dclatds,  ont  dtd  indiqudes  au  paragraphe  1.1.  Des  ameliorations  devraient  dtre  apportees 
dans  ce  cas  comme  en  prdsence  de  spoilers  pour  s'affranchir  des  tarages  prelirainaires  des  veines 
d'essais  &  l'aide  de  plaques  planes. 

Dans  le  cas  de  maquettes  actives,  pour  lesqueHes  les  jets  sont  simulds ,  des  dtudes  seraient 
ndeessaires  pour  effectuer  les  corrections  utiles  tenant  compte  de  la  forme  et  de  1 'emplacement 
des  jets  et  calculer  les  modifications  des  termes  de  volume  et  de  sillage  en  blocage. 

1.2.2  -  Termes  de  portance 

Les  corrections  d' incidence  et  de  c our bur e  de  champ  font  intervenir  les  valeurs  raesurdes 
Ci_  f  Cm  ,  et  x/c  et  des  coefficients  de  correction  prdcalculds  en  fonction  de 

divers  paramdtres.  Actuellement  ces  paramdtres  peuvent  dtre  : 
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-  forne  de  veine  et  porosity  des  parois 

-  envergure  relative  de  maquette  u-  2.S  /JS> 

-  flAche  de  voilure  -A. 

-  type  de  repartition  de  circulation  en  envergure 

-  exccntrement  vertical  de  maquette  ^/d. 

En  outre,  la  position  de  la  nappe  tourbillonnaire  par  rapport  A  la  maquette  peut  Stre 
imposAe,  par  exeraple  au  point  neutre  obtenu  en  esasis. 

Les  influences  de  ces  divers  paramAtrea  sont  illus trees  par  quelques  exemples  donnAs 
figures  11  et  12.  En  ce  qui  concerne  la  repartition  de  circulation  en  envergure  ,  elle  est 

suppcsee  constante  lorsque  1* incidence  varie  dans  les  rAsultats  corriges  donnes  au  paragraphe 
1.1.  Le  calcul  des  coefficients  de  corrections  d'.'ncidence  et  de  courbure  de  champ,  locaux  ou 
integres  en  envergure  <$""  ,  <T-/  ,  peut  etre  effectue  pour  un  cas  quelconque  de  repartition 

et  il  serait  possible  de  tenir  compte  de  1' evolution  de  celle-ci  en  incidence  ai  l'Aquipement 
maquette  permet  d'accAder  A  /YcjJ  . 

Les  exemples  qui  vont  Atre  donnes  ci-aprAs,  issus  d'essais  sur  une  aile  A  flAche  varia¬ 
ble  fioj  montrent  qu'il  y  aurait  lieu  de  tenir  compte  de  Involution  avec  c<  de  la  distri¬ 
bution  /Vy)  • 

La  figure  13  donne  la  courbe  de  portance  globale  de  l'aile  mise  en  flAche  A  60  degrAs 
en  incompressible  en  fonction  de  ctf  .  Cette  courbe,  ainsi  quc  la  polaire,  montre  nettement 

le  point  A  d' apparition  du  tourbillon  d'apex.  De  A  a  B  le  tourbillon  se  developpe  et  commence 
A  dclater  sur  la  voilure  en  B  avec  perte  progressive  de  portance  ;  en  C  1'Aclatement  est  intense 
et  le  trainee  croit  rapideracnt .  La  figure  14  donne  le  rAseau  des  forces  normales  C/y  ,  intAgrAcs 
4  pr rtir  des  pressions  relevAes  sur  chaque  corde,  et>  fonction  de  la  position  Q  de  la  corde  et  .  , 
de  1 'incidence  C<  .  La  perte  de  portance  locale  apparait  A  une  incidence  variant  de  20  A  32 
defies  pour  une  corde  se  ddplaqant  de  1 ‘extrAmitA  a  l'emplanture  ou  la  perte  de  portance  est  plus 
ECiusde.  De  ce  fait,  les  repartitions  plates  A  faible  incidence  deviennent  plutQt  trian- 

gclv ires  dAs  20  degrAs  d'incidence.  Les  figures  15  ft  16  montrent  le  mAme  phAnomAne  A  Mach  0,92 
avec  en  outre  des  discontinuitAs  nettes  des  courbes  C-fy  ,  ey  en  extrAmitA  de  voilure. 

L'influence  de  la  deflexion  de  la  nappe  tourbillonnaire  en  incidence  n'a  paa  encore  AtA 
AtudiAe  car  elle  nAcessite  la  miae  en  oeuvre  de  Bondages  clinoraAtriques  en  aval  de  la  voilure 
raren.ent  ex A cut As. 

Habituellement  les  corrections  d'incidence  et  de  courbure  de  champ  sont  calculAes  en 
valeur  moyenne  par  integration  sur  l'envergure  de  la  maquette.  En  fait,  il  est  possible  d'accAder 
aux  corrections  sur  chaque  corde  de  voilure  et  de  determiner  le  vrillage  de  l'aile  equivalente 
corrigAe  en  atmosphere  illimitAe.  La  figure  17  montre  que  la  correction  locale  d'incidence, 
proportionnelle  A  6o  ,  Avoluera  dif fAremment  en  envergure  selon  le  type  de  repartition  de  cir¬ 
culation  utilise  pour  une  mAme  envergure  relative  T  de  maquette  :  les  variations  do  (y  )  sont 
dans  le  cas  de  la  figure  17,  de  signes  opposes  sur  la  voilure,  sauf  en  son  extreraite.  De  mAme 
lors  d'une  comparaison  de  rAsultats  sur  merae  maquette  dans  des  veines  de  dimensions  diffArentes 
il  y  a  lieu  dc  considArer  que  le  vrillagei  dG  aux  corrections  peut  Stre  different.  Ainsi  pour 
une  mAme  repartition  elliptique  /VyJ  ,  la  figure  17  montre  que  pour  des  envergures  relatives  tr 
de  0,4  et  0,8  le  vrillage  sera  de  signe  oppose  done  que  les  voilures  Aquivalentes  corrigAes  seront 
diffArentes. 

Dans  tout  ce  qui  prAcAde,  la  nappe  tourbillonnaire  est  supposAe  horizontale.  Il  y  aurait 
lieu,  A  grande  incidence  et  dans  les  cas  hypersustentAs ,  de  tenir  compte  de  la  deflexion  die  la 
nappe  tourbillonnaire  A  l'aide  de  mAthodes  telles  que  celles  dAcrites  par  HEYSON  fill  .  Ici  encore 
un  support  experimental  A  partir  de  visualisations  serait  nAcessaire  pour  tenir  compte  d'une 
repartition  en  envergure  de  la  deflexion  et  son  influence  sur  les  corrections  devrait  Atre 
examinAe  avant  d'appliquer  des  raffinements  oe  schAmatisation  maquette. 

1.3  -  Autres  mAthodes  de  corrections  des  effets  de  parois 

1.3.1  -  HAthode_des_pe tits  pavAs 

La  mAthode  des  petits  pavAs,  initiAe  par  JOPPA,  a  AtA  dAveloppAe  par  AMD. BA  fl2*]  .  Cette 

mAthode  permet  de  tenir  compte  des  dimensions  de  la  veine  et  de  la  nature  des  parois  decouples  en 
AlAmenta  rectangulaires  d'intensitA  tourbillonnaire  inconnue.  Il  est  ainsi  possible  de  l'appliquer 
A  des  gAomAtries  quelconques  de  veine,  A  des  longueura  quelconques  de  parois  ventilAes  et  A  des 
emplacements  quelconques  de  maquettes  en  veine.  Bien  entendu  plus  le  descriptif  de  veine  sera  ddtail- 
1A  ,  plus  les  durAes  de  calcul  augmenteront .  Cette  mAthode  requiert,  comme  la  mAthode  analytique,  une 
modAlisation  des  maquettes. 

1.3.2  -  MAth6dc_des_ signatures  des  msquettes 

Devnnt  les  incertitudes  existant  sur  les  conditions  aux  liraites  imposAes  sur  les  parois 
ventilAes  (non  linAaritAs,  amplification  par  les  couches  limites) ,  et  les  limites  d'application  des 
mAthodes  classiques  dans  les  domaines  de  nombre  de  Mach  AlevA  ou  par  suite  de  la  prAsence  de 
dAcollements  plus  ou  moins  organisAs,  ONERA  a  dAveloppA  une  autre  mAthode  £5)  en  relation  avec  le 
concept  dea  parois  auto-adaptable  ^13}  •  Cette  mAthode  est  basAe  aur  la  mesure  des  signatures  de 

la  maquette  sur  les  parois  de  la  veine  d'essais. 
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Le  cslcul  des  ’corrections  a  tout  d'abord  EtE  basE  sur  Is  mesure  de  la  seule  composante 
longitudinale  de  la  vitesse  de  perturbation.  La  maquette  est  encore  reprEsentEe  par  un  nombre 
restreint  de  singularity  proportionnEes  aux  forces  mesurEes  dans  1'hypothAse  d'un  Ecoulement 
compressible  linEsrisE.  De  nombreuses  applications  en  courant  plan  ont  conduit  A  son  utilisation 
en  essais  industriels.  La  figure  18  donne  les  rEsultats  bruts  et  corrigEs  obtenus  dans  la 
aoufflerie  S3  de  Modane  sur  un  profil  dEcrochE .  Le  regroupement  des  rEsultats  est  correct  bien  que 
ler  corrections  en  parois  perforees  soient  ElevEes.  Les  coefficients  de  trainee  Cj>  p  ,  dEduits 
det  integrations  de  pressions  sur  le  profil  s'accordent  bien  aprAs  correction  alors  qu'en  parois 
per;:orEes  les  Cj>p  bruts  dEpassaient  les  trainees  de  sillage.  La  figure  19  donne  les  signatures 
net  :ement  differentes  selon  la  nature  des  parois  ayant  conduit  aux  mSmes  valeurs  de  o(  >  et 

au  voisinage  du  decrochage.  - 

Les  signatures  d'une  maquette  tridimensicnnelle  A  aile  delta  ont  ete  relevees  sur  les 
parois  de  la  veine,  rendue  guidee  ,  de  la  soufflerie  S3  Chalais  .  A  Mach  0,7  et  pour  deux  incidences 
de  10  et  32  degres,  la  demiAre  en  presence  d'un  dEcollement  prononce,  les  courbes  presentees  sort 
issues  des  signatures  :  demi  sommes  (figure  20)  et  differences  (figure  21)  directeraent  liees  respec- 
tivement  aux  tenues  de  blocsge  et  de  portance.  Le  caractAre  tridimensionnel  de  la  signature  en 
portance  (difference)  est  nettement  accuse  A  la  plus  forte  incidence. 

La  methode  des  signatures  est  en  voie  de  generalisation  en  faisant  appel  A  la  connaissance 
suppiementaire  de  Is  direction  de  l'ecoulement  sur  la  surface  de  contr&le.  Dans  le  domaine  compres¬ 
sible  lineaire,  le  champ  de  maquette,  y  compris  les  effets  des  supports  et  des  dEcollements ,  en  Er  ou- 
lement  illimite,  est  deduit  de  l'ensemble  de9  deux  composantes  de  perturbation  sur  la  surface  de 
contc-Qle  et  substituE  aux  champs  des  singularity  representatives  de  la  maquette  lorsque  celles-ci 
sont  incertaines. 

L'extension  au  trans sonique,  domaine  supcrsonique  lEger  inclus,  en  Ecoulement  bidimens ionnel, 
faissnt  appel  aux  mAmes  principes  et  aux  mEthodes  de  petites  perturbations  transsoniques  sera  dEvalop- 
pEe  prochsinement  par  ONERA.  L' spplicstion  locale,  et  non  plus  seulement  globale,  des  corrections, 
qui  apparaitra  nEcessaire  lors  des  domaines  linEsires,  a  4 tE  traitEe  sur  quelques  cas  subcritique*, 

1.3.3  -  Parois  sutoadaptables 

L'extension  de  la  mEthode  des  signatures  au  contrGle  de  la  composante  transversale  d>.  a 
vitesse  de  perturbation  sur  une  surfece  voisine  des  parois,  en  fonction  de  la  composante  longitud.'- 
nale,  permet  d'assurer  1 'Elimination  des  effets  de  parois.  L'application  de  ce  principe  a  EtE 
efiectuEe  dans  la  soufflerie  T2,  au  moyen  de  parois  pleines  dEformables  (l3l  .  L' adaptation  des 
parois  en  Ecoulement  bidimensionnel  et  le  calcul  Eventuel  de  corrections  rEsiduelles  par  la  mEth  >de 
des  signatures  ont  Ete  rEalisEs  sur  le  profil  NAl.  4  0012  A  T2  et  seront  completEs  par  des  esssis  sur 
des  profils  CAST  7. 

En  Ecouleraent  tridimensionnel,  l'application  des  parois  auto  adaptables  soulAve  des  diffi- 
cultEs  pratiques,  Une  adaptation  partielle  des  pa  ’ois  horizontales,  alliEe  A  des  corrections 
rEsiduelles  basEes  sur  les  signatures,  complAtere  les  essais  effectuEs  en  tridimensionnel  A  S3Ch 
avec  des  parois  non  dEformEes ,  indiquEes  au  para,  raphe  prEcEdent, 

II  -  INFLUENCE  DES  SUPPORTS  DE  MAQUETTE S 


Les  supports  de  maquettes  ne  sont  en  principe  pas  considErEs  dans  les  corrections  des  ef¬ 
fets  de  parois.  Le  principe  de  superposition  des  trois  Ecouleraents  suivants  est  admis  :  Ecoulement 
autour  de  la  maquette,  isolEe,  en  atmosphere  illimitEe,  Ecoulement  induit  par  les  parois  et  Ecoule- 
ment  perturbsteur  dG  A  la  prEsence  des  supports  de  msquettes  en  prEsence  des  parois. 

Des  calculs  thEoriques  du  champ  rAgnant  autour  des  supports  de  maquettes  sont  souvent 
effectuEs  A  psrtir  de  multisingularitEs  sur  la  pEriphErie  des  supports.  La  difficultE  dans  le  cas 
des  dards  placEs  A  grande  incidence  rEside  dans  le  calcul  du  champ  perturbateur  incluant  les 
sillsges  en  aval  des  dards.  En  l'absence  de  tels  calculs,  des  sondages  de  1 'Ecoulement  autour  des 
supports  A  1' emplacement  maquette  restent  encore  la  meilleure  mEthode  permettsnt  d'accEder  aux 
corrections induites  par  les  supports,  En  gEnEral  seules  des  mesures  de  rEpartitions  longitudinales 
de  pressions  sont  effectuEes  pour  accEder  aux  corrections  de  poussEe  d'ArchimAde  D  ^  nEcessai- 

res  pour  Etablir  le  bilan  de  trainEe  d'un  avion  civil.  Des  sondages  clinomAtriques  complEmentaires 
perraettraient  d'accEder  A  la  perturbation  dEflexion  sur  empennage  horizontal. 

Un  exemple  des  rEsultats  de  sondages  en  pressions  A  Mach  0,8  dans  la  soufflerie  S2  de 
Modane  est  donnE  figure  22.  Les  rEpartitions  longitudinales  de  pressions,  le  long  de  la  droite 
qu 'occuperait  l'axe  maquette  aux  diffErentes  incidences,  sont  mesurEes  A  l'aide  d'une  sonde 
multiprises  en  l'absence  et  en  prEsence  du  dard  support  de  maquette.  MSme  A  de  faibles  incidences, 
la  variation  du  gradient  de  pression,  entralne  des  variations  des  corrections  de  pou'ssEe  d'ArchimAde, 
sur  la  partie  arriAre  du  fuselage,  non  nEgligeables  avec  M  et  C*  . 


Dans  un  tout  autre  domaine  d' incidence,  une  Etude  de  1' influence  des  dards  supports  de 
missiles  a  EtE  effectuEe  en  subsonique  ElevE  dans  la  soufflerie  S3  de  Modane  £14}  .  Les  rEsultats^ 

bien  que  non  gEnEralisables  A  d'autres  configurations  de  msquette,  donnent  un  ordre  de  grandeur  des 
corrections  A  apporter  aux  rEsultats  de  mesure.  Les  rEductions  de  force  normale  provoquEes  par 

la  prEsence  du  support  de  dard  prolongeant  la  maquette  sont  infErieures  A  2  */•  A  Mach  0,7  et  attei- 
gnent  7  %  A  Mach  0,9  (figure  23).  Le  dsrd  de  balance,  prolongeant  le  missile,  mSrae  avec  un  diamAtre 
infErieur  A  celui  du  fuselage  (ici  0,65  D) ,  induit  un  lEger  supplEment  de  portsnee  par  attEnuation 
des  effets  d'extrEmitEs  sur  le  fuselage  en  incidence.  II  y  a  lieu  de  noter  (15}  qu'un  montage 
sur  mAt  dorsal  rEduirait  ,  A  1' inverse  du  dard  droit,  et  d'une  quantitE  importante  par  altEra- 
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tion  de  1 'organisation  tourbillonnaire  du  sillage  du  fuselage.  La  simulation  du  jet  du  propulseur 
(figure  23)  provoque,  k  partir  de  40  degrds  d'incidence,  un  accroissement  de  C*,  de  l'ordre  de 
10  /«  k  60  degrds,  done  un  effet  de  signe  opposd  k  ce lui  du  support  d'incidence. 

Le  champ  du  support  de  dard  a  pour  effet  d'augmenter  les  incidences  locales,  davantage  sur 
les  gouvernes  que  sur  les  ailes,  et  de  rdduire  les  vitesses  de  l'dcoulement.  Ce  ralentissement  est 
confirmd  d'une  part  par  1 'accroissement  de  la  press  ion  au  culot  de  la  maquette  et  par  la  signature 
de  1* ensemble  maquette-dard  mesuree  sur  le  plancher  de  la  veine  d'essais  (figure  23)  donnd  k  titre 
d'exemple  pour  OC  *  60°.  Cependant,  le  champ  des  sous-vitesses  locales  induit  par  le  support, 
qui  a  tendance  k  rdduire  les  charges  adrodynamiques  sur  la  maquette  est  preponderant  vis  k  vis  de 
l'jffet  d 'augmentation  de  1 'incidence  moyenne,  ces  deux  rdsultats  ayant  des  effets  contraires .  r 
L' influence  du  jet  du  propulseur  sur  la  signature  au  plancher  de  la  veine  est  ndgligeable  ce  qui 
setnbl.*  indiquer  que  le  jet  ne  erde  qu'une  modification  des  incidences  locales  sans  altdrer  les 
vitesses  de  1 'dcoulement. 


Ill  -  vibrations  des  maquette s  autour  du  decrogiage 

Les  maquettes,  montdes  en  dard,  sont  parfois  soumises  k  des  vibrations  de  fortes  amplitudes 
autour  et  au  delA  du  ddcrochage.  Outre  que  le  probldme  de  la  sdcuritd  des  maquettes  rdquiert  une  re¬ 
duction  des  amplitudes  de  vibrations,  des  filtrages  approprids  des  raeaures  s 'avdrent  ndeessaires 
pour  pouvoir  rdduire  la  dispersion  des  rdsultats  et  fournir  au  constructeur  les  meilleures  caractd- 
xistiques  adrodynamiques  moyennes  au  delk  du  ddcrochsge, 

Au  cours  d'essais  k  basse  vitesse,  dans  la  soufflerie  pressurisde  FI,  sur  une  maquette 
d' avion  civil  montde  en  dard,  des  vibrations  imports ites  sur  le  mode  de  flexion  fondamentale  du 
dard  ont  dtd  rencontrdes  pour  certains  cas  de  configuration  d 'atterrissage.  L' evolution  dans  le 
temps  des  amplitudes  conduiaait  alors  1 ' experiment at eur  k  limiter  le  domaine  des  incidences 
expJordes , 

Une  approche  de  calcul  de  prdvision  des  amplitudes  a  dtd  dlaborde  k  l'ONERA  :  elle  est 
bar  da  sur  la  seule  connaissance  des  coefficients  adrodynamiques  stationnaires .  L'dquation  gdndrale 
du  mouvement  est  dtablie,  sans  approximation,  pour  des  mouvements  d 'amplitude  moderde  autour  d'une 
incidence  moyenne  dlevde  variant  avec  le  temps  seJ >n  un  cycla  ddfini.  La  rdsolution  de  cette 
dqu^tion,  k  l'aide  des  radthodes  classiques  de  md unique  non  lindaire  conduit  aux  dvolutions  de  la 
frdquence  et  de  1 'amor tis semen t  avec  1 'amplitude  dw  raouvement,  Lorsque  des  instabilitds  sont 
rencx-ntrdes ,  provenant  des  seuls  termes  de  rigiditd*.  adrodynamiques  ndgatifs  au  del&  du  ddcrochsge, 

1 'effet  cumulatif  des  amortissements  ndgatifs  conduit  k  des  cas  d' amplification  des  amplitudes. 

La  figure  24  donne  ainsi  en  incompressible  pour  trois  dvolutions  types  des  forces  adrodyna- 
miquea,  autour  du  centre  instantand  de  la  ligne  de  dard,  avec  l1 incidence,  les  variations  de  l'amor- 
tissunent  et  de  1' amplitude  qui  seraient  obtenus  au  cours  d'un  cycle  de  montde  et  descente  en  inci 
dence.  Dans  ce  cas  1 'incidence  maximala  est  de  24  degrds  et  la  vitesse  de  balayage  d'incidence  de 
0,5  degrds  par  seconde.  Selon  l'intensitd  de  ddcrochage,  1 'amor tis semen t  devient  plus  ndgatif  done 
l'anriitude  croit.  Le  phdnomdne  se  compliqua  si  des  cycles  limites  peuvent  etre  prdsenta  et  en 
par  iculier  1 'amortissement  est  modifid  par  1' amplitude.  Dc  ce  fait  les  amplitudes  en  incidences 
croinsantes  et  ddcroisqantes  sont  diffdrentes.  En  particulier  1 'effet  cumulatif  peut  Stre  tel 
qu'une  divergence  dynamique  soit  atteinte,  f 


La  figure  25  montre  la  susceptibilitd  des  amplitudes  de  vibration  k  un  parara&tre  essentiel,  E 

la  vitesse  de  balayage  en  incidence,  lorsque  celle-c±  augmente  l'dchelle  des  temps  peut  Stre  telle 
que  les  instabilitds  n'ont  pas  le  temps  de  s'dtablir  et  les  amplitudes  sont,  dans  le  cas  illustrd,  | 

nettement  rdduites.  f. 


Selon  les  niveaux  d 'amplitudes ,  des  filtrages  plus  ou  moins  drastiques  peuvent  2tre  appro- 
prids  pour  accroitre  la  validitd  des  polaires  et  courbes  adrodynamiques  en  valeurs  moyennes. 


CONCLUSIONS  ;  Contrdle  de  la  validitd  des  corrections  des  effets  de  parois 

.  # 

Les  calculs  de  corrections  des  effets  de  parois  appliquds  A  des  maquettes  de  dimensions 
raisonnables  conduisent,dans  bien  des  cas,&  des  comparaisons  de  rdsultats,  obtenus  dans  diverses 
souffleries,  honorables.  Des  corrections  classiques  ddduites  de  la  mdthode  analytique,  ne  ndees- 
sitant  que  des  durdes  rdduites  de  calcul  sont  bien  adaptdes  aux  souffleries  k  caractfere  industriel 
et  effectudes  aystdmatiqueraent  par  ONERA.  Cette  mdthode  est  perfectible  par  amdlioratinn  des 
moddlisations  des  maquettes.  Ainsi  les  dvolutions  en  incidence  des  rdpartitions  de  portance  en 
envergure  et  des  ddflexions  des  nappes  tourbillonnaires  de  voilure  pourraient  Stre  considdrdes 
dfes  lors  qu'elles  seraient  connues .  A  forte  incidence,  dans  le  cas  de  voilure  k  forte  fldche, 
des  modfelea  rdalistes  dea  ddveloppements  du  tourbillon  d'apex  et  de  son  dclatement  restent  k 
ddvelopper.  II  en  va  de  mfcme  pour  la  prise  en  considdration  des  jets  actifa  des  maquettes. 


D'autres  mdthodes  de  corrections  des  effets  de  parois  sont  en  cours  de  ddveloppement .  La 
mdthode  des  petita  pavds,  ddvcloppde  par  AMD, BA  peut  s'avdrer  ndeessaire  pour  des  veines  et  des 
longueurs  de  parois  vcntildes  courtes  ainsi  que  dans  le  cas  dc  formes  compliqudes  de  veines.  La 
mdthode  des  signatures  en  pressions  des  maquettes  sur  les  parois  des  veines  d'essais  prdsente  un 
progrds  intdressant  par  l'affranchissement  des  conditions  limites  lindaires,  criticables  en  parois 
perfordes  et  des  dtalonnages  de  leur  porositd,  De  telles  mdthodes,  nettement  plus  lourdes  en  durdas 
d'essais  ou  de  calculs,  seraient  utilisdes  en  essais  industriels  dans  les  css  ou  la  mdthode  clas- 
sique  s'avdrerait  insuffisante. 


f 


Le  principal  problAme  restant,  est  la  validation  des  mdthodes  de  corrections  des  effets  de 
parols.  En  ce  qui  concerne  les  sch&natisations  des  maquettes,  il  est  toujours  loisiblede  verifier  que 
des  moddlisations  plus  yiabordes  n ' entrainent  que  des  modifications  mineures  des  corrections  done: 
qu'elles  sont  superflues.  Ls  comparaison  entre  des  rdsultats  provenant  de  diverses  origines  reste 
toujours  un  dldment  d£cisif  mais  n^cessitant  une  critique  approfondie  et  dylicate  pour  s'assurer 
que  les  y carts  observes  rdsultent  bien  des  corrections  des  effets  de  parois  et  non  de  1 'influence 
d'autres  paramAtres  parfois  prgponddrants ,  En  outre,  il  y  a  lieu  de  discriminer  la  part  due  aux 
dispersions  de  mesures  et  A  leurs  traitements .  (fig  ire  26) 

La  comparaison  vol  -  soufflerie,  dvidemment  la  plus  int£ressante  en  definitive,  fait  appel 
A  des  corrections  multiples  et  souvent  bien  supdrieures  aux  effets  de  parois  tellesque  les  influences 
des  supports  de  maquette  ["16/  ,  les  transpositions  en  Reynolds  et  la  prise  en  compte  de  1 'aeroy last!- 
cite  QL7j  .  -* 

A  partir  de  nombreuses  comparaisons  vol-soufflerie  effectudes  sur  le  Concorde,  AEROSPATIALE 
conclue  que  les  variations  et  les  grandeurs  des  car&ct£ristiques  adrodynamiques  A  haute  incidence 
ddtermindes  en  soufflerie  sont  retrouv^es  sur  1* avion  en  vol  Q.8J  . 

La  comparaison  entre  rdsultats  sur  mAme  maquette  dans  diverses  configurations  de  parois 
d'une  m£me  veine  d'essais  est  intdressante  par  1 'invariance  de  certains  paramAtres  tels  que  turbu¬ 
lence,  interaction  de  supports,  mais  insuffisante  esc  elle  peut  cacher  une  erreur  systdmatique  si 
une  correction  ddpend  peu  de  l'Atat  des  parois. 


La  comparaison  entre  rdsultats  sur  maquettes  homothdtiques  dans  diverses  souffleries  est 
de  loin  la  plus  dAlicate  car  elle  fait  intervenir  lei  qualitds  respectives  des  ycoulements  de 
diverses  souffleries  et  la  susceptibility  des  maquettes  A  ces  qualitAs.  D'autre  part,  il  y  a  lieu 
de  s'assurer  de  1 1 homo th^ tie  non  seulement  au  sens  macroscopique  mais  encore  au  sens  dtat  de  surf  ice 
et  8 irodlasticitd  des  maquettes.  Enfin  si  les  ytats  de  parois  des  veines  d'essais  comparAes  et  si 
les  tallies  relatives  de  maquettes  different,  les  ytats  de  msquettes  corrigyes  yquivalentes  peuve^c 
d'f.yrer  notablement  en  vrillage  et  carabrure. 


La  comparaison  entre  rysultats  sur  mAme  mavuette  et  mSme  support  dans  diffdrentes  souffle- 
ries  permet  d'yiiminer  des  causes  d'yearts  mais  non  d'obvier  aux  difficultys  rysultant  des  quality^ 
d'ycoilement  spycifiques  A  chaque  soufflerie  et  des  ytats  de  maquettes  corrigdes  yquivalentes. 

Une  perspective  d'avenir,  qui  relAve  indirecteraent  des  mythodes  de  corrections  des  effets 
de  parois,  ddveloppee  depuis  quelques  anndes,  est  le  concept  des  parois  sutoadaptables .Actuellement 
des  perspectives  encourageantes  s'ouvrent  pour  les  applications  en  Acoulement  bidimensionnel. 
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tig.  5  :  Comparaison  des  souffleries  FI  -  i IMA  MERCURE  100  (ech.  1/10  -  configuration  croisiere 
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Influence  de  la  correctiou  de  decollememt  Polaire  -  courbe  CA  (c<  ) 
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Fig.  i  0  s  Comparaison  des  Bouffleries  FI-  SIMA  Maquette  a  forte  fleche  -  configuration  avec 


Effect  of  the  model  schematisation 
on  4he  lift  correction  coefficient 


Fig.  11  :  Influence  de  la  schematisation  de 
voilure  sur  la  correction  d* inci¬ 
dence 


Fig.  12  :  Influence  d'un  excentrement  verti¬ 
cal  de  maquette 


EFFECT  OF  VERTICAL  MODEL  OFFSET 


LIFT  INTERFERENCE  PARAMETER 
ALONG  THE  SPAN 


HALF  MODEL  .CLOSED  WALL 


Fig.  17  :  Repartition  en  envergurp 
de  correction  df inci¬ 
dence 


Fig.  18  :  Ecoulement  2D  -  Corrections 


WALL  INTERFERENCE  FROM  WALL  PRESSURES 
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EFFECT  OF  THE  PARAMETER  A«/At 
ON  THE  VIBRATION  AMPLITUDE 


a" 


Pig,  25  :  Influence  de  la  vitesse  de  variation  d'incidence  sur  l'ampli-- 
tude  des  vibrations 


HOW  ARE  CHECKED  WALL  CORRECTIONS  VALIDITY  ? 

LEVEL  OF  MODEL  DESCRIPTION 
COMPARISONS 

W.T., FLIGHT  TESTS  :R  TRANSPOSITION  f  AEROELASTICITY 
W.T.  '  W.T.  :  SAME  REYNOLDS 

DIFFERENT  MODELS  OR  W.T.  OR  WALLS 

FLOW  QUALITY 

SIMILARITY 

CORRECTED  EQUIVALENT  MODEL 
STING  INTERFERENCE 

Fig,  26  :  Controle  de  la  validite  des  mgthodes  de  corrections 
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GERMAN  ACTIVITIES  ON  WIND  TUNNEL  CORRECTIONS 


H.  Holst 

Deutsche  Forschungs-  und  Versuchsanstalt  ffir  Luft-  und  Raumfahrt  e.V. 
Aerodynamische  Versuchsanstalt  Gottingen 
Bunsenstrasse  10,  D-3400  GSttingen,  W,  Germany 


Abstract 

To  achieve  better  results  from  wind  tunnel  measurements  it  is  necessary  to  carry  out 
interference!  free  measurements  or  to  improve  the  methods  to  correct  data  to  free  flight 
conditions.  This  paper  summarizes  the  main  German  aztivitities  in  this  field. 

o  Wind  Tunnel  interference  factors  have  been  calculated  for  open,  closed,  slotted 
and  perforated  walls  using  the  vortex  lattice  method  with  a  homogeneous  boundary 
condition. 

o  Taking  the  lift  dependent  relocation  of  the  trailing  vortices  into  account  gives  a 
more  realistic  pitching  moment  correction.  Sor.e  results  axe  presented. 

o  For  models  large  in  comparison  to  the  test  section  dimensions,  the  inhomogeneities 
of  lift  and  blockage  interference  parameters  throughout  the  test  section  have  been 
investigated. 

o  A  method  has  been  developed  using  measured  wall  pressures  for  the  correction  of 
drag  in  transonic  wind  tunnels. 

o  For  closed  test  sections,  the  image  method  and  a  modified  vortex  lattice  method 
have  been  used  to  evaluate  wall  pressure  sfcaals  for  correction  purposes. 

o  A  2-d  adaptive  wall  wind  tunnel  has  been  bu'.it  giving  encouraging  results.  Two 
different  concepts  of  a  3-d  adaptive  wall  wind  tunnel  are  being  pursued. 


I.  INTRODUCTION 

In  wind  tunnel  testing  of  models,  the  presence  of  the  boundaries  imposes  a  constraint 
on  the  flow  not  existing  in  free  flight.  To  correct  measured  data  to  free  flight  conditions 
different  Methods  have  been  developed  to  perform  corrections  of:  Angle  of  attack  (flow 
direction  \t  the  position  of  the  wing) ,  pitching  moment  (flow  direction  at  the  position 
of  the  tail)  ,  stagnation  pressure  (blockage  effect) ,  the  wall  induced  curvature  of  the 
flow  (induced  camber  of  the  wing) ,  and  the  inhomogeneities  in  spanwise  direction  (effective 
twist) . 

To  determine  these  corrections,  potential  theory  can  be  used  to  compute  interference 
factors  and  correct  wind  tunnel  data  to  free  flight  conditions  such  as 

o  the  image  method 
o  the  vortex  lattice  method 
o  the  wall  perturbation  method 

o  and  methods  using  measured  wall  pressure  distributions. 

In  most  cases  the  gifting  system  is  represented  by  simple  mathematical  models  such  as 
horseshoe  vortices,  which  are  horizontal  or  follow  a  given  trajectory.  However,  for  high 
lift  configurations  this  is  no  longer  adequate.  The  trailing  vortices  are  relocated  due 
to  the  presence  of  the  wind  tunnel  boundaries  which  means  that  the  trajectory  of  the 
trailing  vortices  in  the  wind  tunnel  is  different  from  the  free  flight  trajectory.  The 
pitching  moment  correction  is  different  from  that  using  horizontal  trailing  vortices. 

The  correction  now  depends  on  the  circulation  of  the  wing  which  is  not  the  case  for 
horizontal  trailing  vortices.  For  the  vortex  lattice  method  computer  programs  do  exist  to 
account  the  relocation  of  the  trailing  vortices  in  the  wind  tunnel  due  to  the  presence 
of  the  tunnel  boundaries  by  means  of  an  iterative  procedure. 

Increasing  the  model  size  and  the  angle  of  attack  (flow  separation,  wake  blockage)  not 
only  leads  to  an  increase  of  the  applied  corrections  but  also  to  an  increase  of  their 
inhomogeneities.  That  means,  e.g.,  that  it  is  probably  no  longer  permitted  to  apply  average 
corrections  for  the  whole  model.  One  has  to  be  aware  of  the  possibility  that  the  flow 
field  around  the  model  has  substantially  changed  due  to  the  presence  of  the  wind  tunnel 
boundaries.  Measured  data  can  then  not  be  corrected  to  free  flight  conditions.  Excluding 
these  cases  (when  the  flow  field  is  substantially  -changed)  it  is  most  important  to  know 
more  about  boundary-induced  effects  and  their  inhomogeneities.  This  knowledge  can  be 
used  to  estimate  the  maximum  tolerable  model  size  for  a  given  wind  tunnel. 

Generally  there  are  three  methods  to  avoid  large  corrections  and  high  levels  of  their 
inhomogeneities . 

o  reducing  model  size  or  very  large  test  sections 
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o  ventilated  test  sections  (longitudinal  slots  and/or  perforated  walls  to  reduce  the 
magnitude  of  the  corrections  and  their  inhomogeneities)  . 
o  adaptive  wall  concept  (streamline  matched  test  section)  to  cancel  the 
corrections  at  all. 

In  the  following  some  German  work  done  in  the  fields  mentioned  above  will  be  presented. 


II.  ANGLE  OF  ATTACK  AND  PITCHING  MOMENT  CORRECTIONS  IN  LOW  SPEED  WIND  TUNNELS  WITH 


ARBITRARY  WALL  CONFIGURATIONS 


a.  Open  and  closed  test  sections 


la  the  past,  corrections  were  computed  using  the  image  technique  representing  the  lift¬ 
ing  system  by  one  single  horseshoe  vortex.  Figure  1  shows  the  system  of  images.  -The  inter¬ 
ference  factor  is  defined  in  the  following  equation: 


boundary  induced  vertical  velocity 
spanwidth  of  the  model 
cross-sectional  area  of  tunnel 
wing  circulation. 


The  angle  of  attack  correction  is  then  given  by 


'  5o  a;  cl 


wing  area 
lift  coefficient 


5  au  well  as  its  derivative  in  x-direction  5.  =  d5  /dx  was  calculated  at  the  location  of 
the  ving.  The  interference  factor  at  the  tail  position  (pitching  moment  correction  could 
then  be  found  very  simply  by  extrapolation,  Figure  2  and  Ref.  1. 

•  '  *  2 

The  vortex  lattice  method  for  closed, test  sections,  reported  by  Joppa  in  1967  and  for 
open  test  sections  reported  by  Doberenz4  in  1974,  has  now  been  introduced  for  interference 
calculations  in  a  version  reported  by  Borovik4  in  1973.  This  method  can  compute  interference 
factors  for  open,  closed,  slotted  and  perforated  test  sections.  This  program  can  1  andle 
arbitrary  cross  sections;  lift  distributions  and  sweep-back  angles  of  the  wing  can  be 
prescribed.  The  representation  of  the  tunnel  boundaries  by  vortex  squares  or  doublet  sheets 
is  shown  in  Figure  3.  The  lifting  system  is  a  set  of  horseshoe  vortices  shown  in  Figure  4. 
Calculations  have  been  performed  to  check  the  validity  of  the  method.  The  results  have 
been  compared  (for  closed  test  sections)  with  results  of  the  Joppa^  computer  program 
modified  by  Holst®. 

Figure  5  and  Figure  6  show  the  variation  of  interference  factor  in  longitudinal  and 
spanwise  directions,  respectively;  parameter  is  the  sweep-back  angle  A.  In  1979  Schulz® 
presented  computational  results  of  interference  factors  for  closed  rectangular  test  sections 
using  the  image  method.  Figure  7  shows  a  typical  result  which  gives  an  impression  of  the 
inhomogeneity  of  the  wall  interference  in  the  main  cross  section  (model  position) . 

Knowing  the  wall  interference  throughout  the  test  section  (at  least  along  the  span  of  the 
wing  and  the  span  of  the  tail)  one  can  find  average  values  for  the  correction  factors  of 
wing  and  tail.  Corrections  of  angle  of  attack  and  pitching  moment  including  induced  camber 
and  effective  twist  can  then  be  performed. 


b .  Slotted  wall  test  sections 

As  mentioned  above  the  vortex  lattice  method  is  capable  of  calculating  interference  5 
factors  for  slotted  and/or  perforated  walls.  Calculations  were  performed  to  investigate 
the  influence  of  slot  geometry  parameter  K  and  porosity  parameter  R,  respectively.  Compa¬ 
risons  were  made  with  computational  results  of  an  NLR  program?  using  sources  and  sinks  to 
represent  the  wind  tunnel  boundaries.  Both  methods  use  a  homogeneous  boundary  condition 
which  is  shown  in  Figure  8  for  various  wall  configurations.  Figure  9  shows  results.  The 
agreement  of  the  two  methods  is  quite  good  for  small  values  of  Q  (Q  =  (1  +  1/Rf”1  )  . 

Figure  10  and  Figure  11  show  the  influence  of  slot  geometry  parameter  K  and  porosity 
parameter  Q,  respectively  (see  Appendix  A  for  definition  of  parameters  and  homogeneous 
boundary  condition) . 

Measurements  were  carried  out  in  a  wind  tunnel  with  a  square  cross  section  of 
1,3  m  X  1,3 -m  with  longitudinally  slotted  walls  (12%  open)  at  DFVLR-HNW  Braunschweig.  The 
experimental  results  indicated  that  the  effect  of  the  slots  was  smaller  than  predicted 
by  theory.  This  is  probably  caused  by  a  too  small  plenum  chamber,  and,  therefore,  experi¬ 
ments  have  begun  with  no  plenum  at  all  (i.e.,  very  large  plenum).  It  is  intended  to  carry 
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out  measurements  in  the  DNW  (German  Dutch  Wind  Tunnel)  with  an  8  m  x  6  m  rectangular  test 
section  equipped  with  longitudinally  slotted  walls  (0  -  12%  open)  to  check  the  existing 
theories . 

C.  Vortex  relocation  due  to  the  presence  of  the  tunnel  boundaries 

For  an  improvement  of  pitching  moment  corrections  for  high  lift  configurations  it  might 
be  necessary  to  take  into  account  the  vortex  relocation.  In  1973  Joppa®  presented  a  vortex 
lattice  method  using  a  single  horseshoe  vortex  as  lifting  system.  He  takes  into  account 
the  vortex  relocation  due  to  the  presence  of  the  wind  tunnel  boundaries  (closed  test 
section) .  Equilibrium  trajectories  for  the  trailing  vortices  are  found  for  the  wing  in 
tunnel  and  ir.  free  flight;  the  difference  of  the  induced  velocity  fields  is  due  to  the 
presence  of  the  wind  tunnel  walls.  The  horseshoe  vortex  model  used  for  the  iterative  Pro¬ 
cedure  is  shown  in  Figure  12.  This  method  was  extended  to  open  test  sections  by  Holst? 
using  doublet  sheets  instead  of  vortex  squares  to  represent  the  wind  tunnel  boundaries 
and  the  wing.  Figure  13  shows  a  sketch  of  the  expected  wake  trajectories  for  open  and 
closed  tunnels.  Figure  14  and  Figure  15  give  some  typical  results^0.  Figure  16  shows 
the  interference  factors  with  and  without  vortex  relocation  for  different  values  of  the 
downwash  factor  c  /AR  for  an  open  test  section  (c  =  lift  coefficient,  AR  =  aspect  ratio 
of  wing)  u  .  u 


III.  BLOCKAGE  CORRECTIONS 

11  12  13 

In  various  publications  G.  Schulz  '  '  reported  about  the  inhomogeneities  of  solid  . 

blockage  corrections  throughout  the  test  section.  Using  the  image  technique  and  superimpo¬ 
sing  point  singularities  representing  the  influence  of  model  blockage,  he  computed  t, 
defined  in  the  following  equation: 


— &Y 


K,1  model  shape  parameters 
V  volume  of  the  model 
S™  area  of  wing 

a”  cross-sectional  area  of  the  tunnel 
Av  boundary  induced  longitudinal  velocity 
V„  free  stream  velocity 

T  (=  ■j'j/ir' t)  tunnel  shape  parameter 

Figure  17  shows  an  iso-line  diagram  of  tin  the  main,  cross  section  (i.e.,  x  «  xmoiel'’ 
Schulz  investigated  the  influence  of  the  following  parameters: 

span  width/tunnel  width  b/W 

length  of  fuselage/tunnel  width  1,/W 

eccentricity  of  model  (model  position)  y,z 
tunnel  height/tunnel  width  H/W 

spanwise  distributions  of  volume  V(y) 

A  variety  of  computational  results  is  presented  in  the  references  mentioned  above. 
Figure  18  and  Figure  19  show  some  results  which  can  be  used  to  determine  the  maximum  to¬ 
lerable  model  size  for  a  given  wind  tunnel  or  to  estimate  the  achievable  accuracy  of 
measurements  (so  far  as  blockage  corrections  are  concerned) . 


IV.  CORRECTIONS  USING  MEASURED  WALL  PRESSURES 

a.  Simple  blockage  corrections  in  closed  subsonic  wind  tunnels 

1 9 

In  1978  Mercker  and  Fiedler  presented  a  simple  blockage  correction  method  for  flows 
around  bodies  of  arbitrary  shape.  The  correction  factor  n  is  obtained  experimentally  by 
measuring  the  wall  pressures  at  appropriate  locations.  The  geometrical  blockage  ratio 
S  /Afc,  the  thickness  of  the  boundary  layer  at  the  model,  the  thickness  of  the  boundary 
layer  at  the  tunnel  walls  and  the  model  shape  do  not  enter  into  the  correction  scheme 
explicitly.  They  are  taken  care  of  by  the  correction  factor  n.  The  method  can  be  applied 
for  total  ratios  of  blockage  up  to  25%,  including  boundary  layer  etc.  It  has  been  tested 
for  flows  about  two-dimensional  wedges  and  cylinders.  In  principle  the  method  is  also 
applicable  to  flows  around  three-dimensional  models.  For  more  details  see  Ref.  19. 

b.  Drag  measurements  in  transonic  wind  tunnels 


In  order  to  increase  the  accuracy  of  drag  measurements  of  slender  bodies  of  revolution 
at  zero  angle  of  attack  in  transonic  wind  tunnels  a  correction  method  using  measured  wall 
pressures  was  proposed  by  MBB  (Aulehla  and  Venghaus) 1 ^ < 1 ? ' 1 ®  .  This  method  can  correct 
pressure  distributions  measured  in  a  wind  tunnel  to  free  flight  conditions.  A  theoretical 
wall  pressure  distribution  (free  flight,  imaginary  wall)  is  computed  by  a  finite  element 
method  and  compared  with  the  wall  pressure  distribution  which  is  measured  together 
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with  that  of  the  model.  The  difference  of  measured  and  computed  wall  pressures  can  be  used 
to  find  a  perturbation  potential  thus  allowing  to  correct  the  model  pressure  distribution. 
Figure  20  and  Figure  21  show  the  principle  of  this  method.  Figure  22  shows  one  example  of 
a  corrected  pressure  distribution. 

.  Interference  free  data  are  necessary  to  prove  the  validity  of  this  correction  method 
and  possibly  extend  it  to  models  at  angle  of  attack. 

o .  Correction  method  using  measured  wall  pressures  - 

For  large  models,  high  angles  of  attack  and  generally  in  cases  with  flow  separation 
the  computational  correction  methods  based  on  potential  theory  are  somewhat  uncertain. 

As  a  contractor  of  DFVLR-HNW  G.  Schulz  develops  a  correction  method  using  measured  wall 
pressures . 

o  This  universal  method  of  correction  is  supposed  to  give  correction  factors 

for  angle  of  attack,  curvature  of  the  flow,  pitching  moment,  and  stagnation  pressure. 

o  The  method  should  provide  reliable  corrections  using  only  geometric  parameters 
of  model  and  wind  tunnel. 

The  problems  that  arise  are 

o  to  properly  position  the  static  pressure  holes  at  the  walls 
o  to  separate  the  different  corrections 

o  to  achieve  high  accuracy  of  computation  as  well  as  of  the  pressure  measurements. 

Schul’  uses  the  image  method  and  simple  mathematical  models  to  represent  the  lifting 
system  and  the  influence  of  blockage  (horseshoe  vortices,  doublets,  sources  and  sinks)  to 
find  influence  functions.  Figure  23  shows  a  vail  pressure  distribution  (influence  function) 
in  the  main  cross  section  for  a  single  horseshoe  vortex.  The  main  advantage  of  the  method 
proposed  is  that  it  takes  into  account  the  solid  blockage  and  wake  blockage  automatically. 
It  is  not  necessary  to  know  more  about  the  nndel  and  the  flow  field  in  its  vicinity  than 
just  a  few  geometric  parameters.  This  method  is  restricted  to  closed  test  sections. 

It  should  be  mentioned  that  the  image  method  is  adequate  for  evaluation  of  pressures 
at  the  tunnel  walls  because  the  boundary  condition  is  fulfilled  by  definition.  The  singu¬ 
larities  used  to  represent  the  model  and  its  Images  are  at  some  distance  from  the  walls. 
This  is  not  the  case  for  the  vortex  lattice  methods  where  the  boundary  condition  is  just 
fulfilled  at  a  set  of  control  points  and  the  singularities  representing  the  walls  arc 
embedded  in  the  wall.  . 


Schulz 
rimen  .s . 


presented  first  results;  the  validity  of  his  method  should  be  checked  by  expe¬ 


lt  was  mentioned  above  that  the  vortex  lattice  method  is  not  appropriate  for  wall  pressu¬ 
re  evaluations.  In  the  usual  vortex  lattice  method  the  singularities  and  the  control  points 
are  lying  in  the  plane  of  the  wall.  The  boundary  condition  is  fulfilled  at  the  control 
points.  At  points  in  the  plane  of  the  wall  which  are  too  close  to  the  singularities  the 
velocity  field  can  not  be  evaluated. 

On  the  other  side  the  image  method  is  not  appropriate  for  slotted  and  perforated  walls, 
it  is  restricted  to  open  and  closed  test  sections. 

It  was,  theiefore,  tried  to  modify  the  vortex  lattice  method.  This  was  done  by  shifting 
the  plane  containing  the  singularities  outwards  while  keeping  the  control  points  at  the 
same  locations  as  before.  As  far  as  the  angle  of  attack  correction  is  concernd,  this  new 
method  was  successfully  applied  to  closed  test  sections. 

For  the  calculation  of  blockage  corrections  in  closed  test  sections  the  vortex  squares 
were  substituted  by  5  doublets  per  panel,  also  using  this  modified  lattice  approach. 

More  details  of  these  methods  including  first  results  are  given  in  Appendix  B. 
Calculations  for  slotted  and/or  perforated  walls  will  be  carried  out  in  the  future  to 
check  the  validity  of  these  new  methods . 


V.  ADAPTIVE  WALL  CONCEPT 


a.  2-dimensional  adaptive  wall  wind  tunnel 


U.  Ganzer  reported  investigations  in  a  2-dimensional  adaptive  wall  wind  tunnel  of  TU 
Berlin.  Figure  24  shows  schematically  the  adaptive  wall  concept.  The  flow  field  inside  the 
wind  tunnel  is  supplemented  by  a  fictitious  flow  field  outside  the  tunnel.  Knowing  the 
inner  static  wall  pressure  distribution  and  the  shape  of  the  wall  a  pressure  distribution 
for  the  fictitious  flow  field  outside  the  tunnel  can  be  found.  The  pressure  differences 


P  ^inside  ^outside 

can  be  taken  to  define  a  new  pressure  distribution  for  the  fictitious  flow  outside  the 
tunnel . 


c  (n)  +  K  Ac  =  c  (n  +  1) 

•outside  P  ^outside 

with  0  <  X  <  1 . 

This  pressure  distribution  gives  a  new  contour  of  ,thei  wall  for  an  iterative  procedure. 

The  advantage  of  this  method  is,  that  measurements  us' well  as  computations  are  very  easy 
to  perform  due  to  relatively  small  curvature  of  the  wall.  This  method  works  without  know¬ 
ing  anything  about  the  geometry  of  the  model  and  the  flow  field  in  its  vicinity. 

Results  of  pressure  distribution  measurements  are  shown  in  Figure  25  and  Figure  26 
demonstrating  the  validity  of  the  adaptive  wall  concept  by  comparison  with  nearly  inter¬ 
ference  free  results.  It  should  be  noted  that  difficulties  arise  for  flows  with  high 
pressure  gradients  or  shock  waves. 

More  recent  pressure  distribution  measurements  v.ith  a  NACA  0012  model  of  higher 
accuracy  than  that  used  before  show  even  better  results  than  those  mentioned  above.  The 
iterative  procedure  of  wall  adaptation  is  fully  automated.  The  computer  program  for  the 
fictitious  outer  flow  field  has  been  optimized,  thus  reducing  the  computation  time  from 
initially  30  sec  for  one  iteration  to  now  8  sec  on  a  TR  440. 

Theoretical  work  for  2-dimensional  investigations  is  carried  out  with  a  subsonic-panel- 
program  (TU  Berlin)  and  the  transonic  Bauer-Garabe dien-Korn-program.  This  is  necessary  tg 
estimate  the  magnitude  of  Reynolds  number  effects  (measurements  at  Re  =  1.0  ./.  1.6  x  10, 
comparable  measurements  mostly  at  Re  =  6  x  10®)  ard  to  investigate  the  influences  of  model- 
accuracy  and  remaining  wall  interference  respectively. 


b.  3-dimens jonal  investigations 

Two  different  models  will  be  investigated:  the  body  of  revolution  ONERA  C5  and  the  wing 
fuselage  combination  ZKP-F4 .  Both  models  are  being  constructed  now.  Measurements  will  be 
performed  in  the  2-dimensional  tunnel  and  in  the  3-dimensional  octagonal  tunnel  with  8 
flexible  walls  shown  in  Figure  27.  It  is  expected  that  the  automatic  procedure  of  wall 
adaptation  for  the  2-d  test  section  will  be  applicable  to  the  3-d  case.  Theoretical  in¬ 
vestigations  are  carried  out  for  estimation  of  the  accuracy  of  computation  for  the  3-d 
outer  flow  field  and  for  simulation  of  the  iterative  procedure  of  wall  adaptation. 


c.  Another  3-d  concept 

At  DFVLR  Gottingen  a  feasibility  study  is  being  made  for  a  different  concept  of  a  3-d 
adaptive  wall  test  section.  This  test  section  is  planned  to  be  a  0.8  m  diameter  deformable 
rubber  tube.  An  existing  high  speed  intermittent  wind  tunnel  could  then  easily  be  equipped 
with  these  adaptive  walls.  It  is  planned  to  optimize  the  iterative  procedure  for  the 
acquisition  of  one  data  point  within  40  seconds. 


VI.  CONCLUSIONS 

The  magnitude  and  inhomogeneity  of  corrections  set  a  limit  to  the  model  size.  For  real 
time  application  during  wind  tunnel  tests  appropriate  average  values  should  be  taken  for 
the  correction  factors. 

A  better  pitching  moment  correction  could  be  achieved  by  taking  into  account  the  vortex 
relocation  because  the  interference  factors  are  lift  dependent.  Time  consuming  computa¬ 
tions  are  necessary,  so  that  this  correction  cannot  be  applied  in  real  time. 

Ventilated  test  sections  reduce  the  magnitude  of  interference  factors.  The  inhomogenei¬ 
ties  of  interference  will  be  investigated  in  the  future.  Measurements  to  check  theory 
have  begun. 

A  promising  approach  to  correct  pressure  distributions  with  the  aid  of  simultaneously 
measured  wall  pressures  has  been  shown. 

A  method  for  corrections  of  angle  of  attack,  pitching  moment  and  stagnation  pressure 
using  measured  wall  pressures  is  being  developed  for  closed  test  sections . 

The  methods  mentioned  above  find  their  application  in  the  range  of  "no  substantial 
change  of  the  flow  field  around  the  model"  compared  to  free  flight  conditions.  This  sets 
a  limit  to  or  is  a  criterion  for  estimating  maximum  tolerable  model  size.  In  some  cases 
the  adaptive  wall  wind  tunnel  will  be  the  only  way  to  get  corrigible  or  even  interference 
free  results. 
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Existing  2-d  adaptive  tunnel  results  are  encouraging.  A  3-d  octagonal  adaptive  tunnel 
is  being  built  in  Berlin.  The  GSttingen  high  speed  wind  tunnel  is  being  equipped  with  an 
adaptive-rubber-tube-test-section. 
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Appendix  A 


Homogeneous  boundary  conditions  (compare  Ref.  4) 

a.  Longitudinally  slotted  wall 

For  an  ideally  slotted  wall  the  pressure  at  the  Slots  is  assumed  constant  and  equal  to 
the  free  stream  pressure,  there  is  zero  normal  flow  at  the  strips. 

=  •  u  =0  at  the  slot  (open  boundary) 

3  x 

=  v  =  0  at  the  strips  (closed  boundary) 

3n  n 

4  velocity  potential 

u  velocity  component  parallel  to  free  stream 
•  vn  velocity  component  normal  to  the  wall. 

The  slot  geometry  parameter  K  is  defined  according  to 

„  ,  1  .  r  .  ,7ra,  ,  1  slot  separation 

*  *  "  i  tsln  2"l  *  a  slot  width. 

An  approximate  average  boundary  condition  for  a  slotted  wall  was  derived  by  Baldwin, 
Turner  and  Knee h tell 8 : 

14  +  K  lit.  =  0 
3X  Vn 

or  using  the  velocity  components 

3v„  ■  ■ 

u  +  K  — -it  «  0 
3x 

This  is  the  boundary  condition  for  an  ideally  slotted  wall,  disregarding  any  viscous  effe.  ts 
in  the  slots. 

b.  Perforated  walls 

The  linearized  approximation  to  a  viscous  flow  through  a  porous  medium  is: 


3*  R  . 

3n  '  vn  pV0  Ap 


p  density 

Ap  pressure  drop  across  the  wall 
R  =  ^  1 

r  porosity  parameters 


Introducing  the  compressiblity  correction  factor  6  =  Vi  -  M  gives  the  boundary  condition 
for  an  ideally  perforated  wall:  . 

u  +  |  vn  =  °. 

The  value  for  the  porosity  parameter  R  must  be  found  experimentally, 

c .  Homogeneous  boundary  condition 

A  homogeneous  boundary  condition  is  fulfilled  at  a  set  of  control  points.  It  can  be 
written  in  the  following  way  to  give  a  general  boundary  condition 

a  14  +  K  s2*  +  1  14  =  n 

3x  3x3n  R  3n 
or,  using  the  velocity  components 


"V^ 
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9v 

Au  +  K  IF1  +  I  vn  -  °- 

This  general  boundary  condition  covers  all  cases  of  test  section  boundaries  with  the  proper 
constants  for  A  (0.0  or  1.0),  K  and  R. 


APPENDIX  B 

Modified  vortex  or  doublet  lattice  approach 
a.  Angle  of  attac  correction 

For  the  evaluation  of  wall  pressures  the  vortex  lattice  method  has  been  modified. 
This  is  shown  in  the  following  sketch. 


The  planes  containing  the  vortex  squares  are  shifted  outwards;  the  distance  from  the 
tip  of  the  wing  to  the  wall  is  the  same  as  the  distance  from  the  wall  to  the  vortex  lattice. 
The  control  points  are,  of  course,  kept  at  the  same  locations  before,  i.e.,  in  the  plane 
of  the  tunnel  walls.  The  homogeneous  boundary  condition  for  slotted  or  perforated  walls 
can  easily  be  applied  at  the  control  points. 

For  closed  test  sections  this  method  has  been  proved  to  give  correct  results  as  is 
shown  in  the  following  plot  of  the  wall  pressure  influence  function  Sxw/6zo= 


wAfc 

2bT7 


vxAt 


6xw  ”  2bT 


w 


boundary  induced  vertical  velocity 
wing  circulation 

test  section  cross  sectional  area 


(at  the  position  of  the  model) 

(at  the  walls) 

V  boundary  and  wing  induced  velocities 
x  in  x-direction  at  the  walls 

b  spanwidth 


•  -  >/» 

nlllln  tftmUtt 


—V 


b.  Blockage  correction 

In  the  case  of  calculating  blockage  interference  parameters  the  lifting  system  (repre¬ 
sented  by  a  single  horseshoe  vortex  or  a  set  of  horseshoe  vortices)  is  substituted  by  a 
doublet  or  a  set  of  doublets  to  represent  model  blockage.  The  vortex  squares  are  substitu¬ 
ted  by  "doublet  squares"  with  5  doublets  each.  The  result  is  presented  in  the  following 
figure  as  an  isoline  diagram  of  t ,  A  single  doublet  represents  the  model. 

z  /H 


This  new  method  will  be  applied  to  slotted  and  perforated  walls  using  a  homogeneous 
boundary  condition. 


Fig.  It  System  of  Images  for  an  open  test  section 
h  (  x  ■  xw|ng  main  cross-section). 


W  width  of  test  section 
H  height  of  test  section 


\  spanwldth  of  model 

1  H  height  of  model 


f 
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Fig.  Mathematical  model  of  the  wing 


'V# 


0,30 


Fig.  6:  Interference  factors  along  span  for  a  ‘ 
square  closed  test  sectle*',  comparison 
of  constant  and  elliptic  lift  distribution. 
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b  spanwidth  of  horseshoe  vortex 

W  width  of  wind  tunnel 

H  height  of  wind  tunnel 

.A  0°  sweepback  angle  of  wing 

k  M  0,7  relative  spin  of  horseshoe  vortex 
W 


i 

l 

■i 


1.  closed  walls 


2.  ideally  perforatsd 


3.  ideally  slotted 


slotted  wall  with 
viscous  effects 


5.  open  boundaries 


u  + 


u  + 


ti 


VL 


9  Vrt 
3  x 


+  4-Vn  .  0*) 


u  induced  velocity  at  ths  wall 

(  x— direction,  parallel  to  V®  and  ths  wall) 

induced  velocity  at  the  wall  Cnormal) 


Flg»  Ss  Homogeneous  boundary  condition  at  ths 

wall  for  various  cases  C  compare  appendix  A  ) • 


In  this  case  the  porosity  parameter  R  takes  into 
account  viscous  cross  flow  In  the  slots, see  also 
appendix  A. 
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f*9«  9*  Comparison  of  vortex  lattice  method  (ref.  5) 
with  the  method  using  sources  and  sinks 
(ref.7)  and  Heyson  program,  ventilated 
walls,  square  cross  section. 


b  ^9  geometric  spanwldth  of  model 
W  width  of  wind  tunnel 

H  height  of  wind  tunnel 

-i.  *0,75 (rectangular  cross  section) 

-k_agp^  0,6  relative  span 

W 

.A  *0°  sweepback  angle 


elliptic  lift  distribution 


Fig.  12:  Representation  of  Lifting  system  (horseshoe  vortex)  for  iterative  procedure 
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Wind  tunnel  boundaries 

w 

Lifting  system  „ 

in  a  closed  test  section 
in  free  flight 
in  an  open  test  section 


pig.  Relocation  of  the  trailing  vortieea 

due  to  the  preeenee  of  wind  tunnel  boundaries 


-on 


Fig.  14:  Interference  factors  at  wing  and  tail .including 

relocation  effects  of  the  trailing  vortlcesj 
tail  on  tunnel  centerline. closed  test  section. 

fcy  _  X  i3eo  -  0,?y  s panwidth  of  horseshoe  vortex 
W  ■  i(5  width  of  wind  tunnel 
H  —  -ftO  height  of  wind  tunnel. 

Lift  coefficient 
/R  aspect  ratio 

c j_/AR  downwash  factor 
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W  tunnel  width 
H  tunnel  height 

b  spanwidth 
^  *  0,6  relative  span 


y  jo-4 


Fig.  17:  Open  test  section^ so-llnes  of  parameter'?' In  the  main  cross-section  (  x-x^^), 
constant  volume  distribution  along  the  span  represents  simple  airfoil  model. 


*  ST0'* 


"ground 

°'W  A„./A  t 


Fig.  18:  Inhomogeneity  of  stagnation  pressure  for 
bulky  models,  open  test  section,  with 
and  without  ground  simulation. 

Constant  volume  distribution  along  span 
represents  simple  airfoil  model . 

H  tunnel  height 
W  tunnel  width 

-tL  *  o,7 
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Fig.  19: 
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Inhomogeneity  of  stagnation  pressure  for 
free  jet  with  ground  simulation,  the  model 
is  a  body  of  revolution  with  a  thin  airfoil. 
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20:  Correction  method  using  measured  wall  pressures 
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1.  The  measured  wall  pressure  distribution  is  influenced 
by  the  presence  of  the  wind  tunnel  walls. 

2.  The  difference  between  measured  and  computed  pressure 
distributions  A epalong  the  tunnel  wall  or  a  corresponding 
imaginary  wall  is  a  criterion  for  the  magnitude  of  the  mistake 

3.  and  can  be  used  to  compute  the  correction  an  the  tunnel  centerline 

4.  The  correction  on  the  surface  of  the  body  corresponds 
to  the  correction  on  the  tunnel  centerline. 


Fig.  21:  Correction  method  using  measured  wall  pressures. 
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Fig.  23:  Wall  pressure  distribution  for  a  closed  square  tunnel,  stain  cross  section  (x  » 

6q  Interference  factor  (vertical  component  of  Induced  velocity  due  to  tunnel  walls  at  the  position 
of  the  model  -  not  averaged  along  span) 

^Xw  represents  Induced  velocities  In  x-directlon  of  model  and  tunnel  walls 


Fig.  24  The  principle  of  adaptive  walls. 

there  Is  no  pressure  jump  through  thewall(f1ct1tious  flow  field  outside  the  tunnel  with  computed  pressure  distribution 
-  real  flow  field  In  the  tunnel  with  measured  pressure  distribution)  If  the  wall  contour  Is  a  streamline.  The  shape 
of  the  wall  Is  changed  by  an  Iterative  procedure.  If  the  measured  pressure  agrees  with  the  computed  pressure  of  the 
fictitious  outer  flow  field,  the  Inner  flow  field  Is  Interference  free. 
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Pig.  25:  Pressure  distribution  f  a*  NACA  0012,  M  «  0,76,  ot  «  0° 
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Fig,  26s  Preaaure  distribution  for  NACA  0012,  ct  *  6°,  M  *  0 
Comparison  of  different  measurements. 
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SUMMARY 

A  survey  is  given  of  past,  current  and  planned  work  at  NLR  ii  the  field  of  vind-tunnel  vail 
interference  on  models  at  high  angles  of  attack  at  lov-subsonic  speeds .  Among  the  nev  long-term  research 
activities  those  concerning  slotted-vall  test  sections  play  a  dominant  part.  It  is  felt  that  an  approach 
which  makes  use  of  measured  distributions  of  flow  quantities  near  the  vails  is  the  most  promising  one.  It 
may  be  recommended  also  as  a  short-term  solution  for  some  special  vail  correction  problems  arising  from 
modern  lov-speed  vind-tunnel  testing  in  closed  test  sections. 

LIST  OF  SYMBOLS 

A'  =  effective  aspect  ratio  of  wing 

A  *  influence  coefficient  matrix 

b  =  tunnel  width 

=  vector  of  model-induced  influences  at  the  collocation  points 

C  ■  quantity  related  to  static-pressure  difference  between  plenum  chamber  and  upstream  infinity 

(see  Fig,  5) 

C  -  cross-sectional  area  of  test  section 

CD  =  drag  coeTicient  of  model 

=  lift  coefficient  of  model 
h  -  tunnel  height 

K  =  slot  parameter 

Kjj,  Ky  =  slot  parameter,  for  horizontal  and  vertical  tunnel  vails  respectively 
L  =  iength  of  ventilated  part  of  test  section 

L  *  model  lift 

n  =  normal  direction  (to  tunnel  vail) 

Pc  =  static  piessure  in  plenum  chamber 

Pro~  static  pressure  of  undisturbed  flow  far  upstream  « 

q  =  dynamic  pressure 

Aq  =  dynamic  pressure  correction 

R  =  radius  of  rotor  model 

Q,  R  =  porosity  parameter 

Qjj,  Rjj  -  porosity  parameter  for  horizontal  tunnel  walls 
Qy,  Ry.  =  porosity  parameter  for  vertical  tunnel  vails 
s  =  model  semi-span 

t 

5  -  model  reference  area  or  area  of  tunnel  wall 

V  -  velocity  in  test  section 

v  =  wall-induced  upwash  velocity 

x,y,z  =  co-ordinate  system  with  origin  in  test-section  centre 
o  -  angle  of  attack  of  model 

T  =  horseshoe  vortex  strength 

Sq,  6^  =  lift  interference  factor 

Act  =  angle  of  attack  correction 

6  =  blockage  factor 

a  23  vector  of  source  strengths 

<p  *=  perturbation  velocity  potential 

1 .  INTRODUCTION 

The  NLR  research  on  wind-tunnel  wall  interference  on  models  at  high  angles  of  attack  concentrates 
almost  entirely  on  problems  in  the  lov-subsonic  speed  regime.  Especially  in  this  field  radical  changes  in 
test  capabilities  and  measuring  techniques  are  taking  place  by  the  completion  of  the  DNW,  the  large  low- 


speed  Gern  an -Dutch  Wind  Tunnel-,  which  is  entering  the  operational  phase  now  *). 

These  charges  include  firstly  a  larger  variety  of  models  and  test  conditions  in  test  sections  of 
various  dimensions  and  wall  configurations  (closed,  slotted  and  open),  but  secondly  they  comprise 
higher  requirements  on  quality,  versatility  and  productivity.  An  enumeration  of  the  prevailing  design 
principles  and  of  the  expected  activities  is  given  in  figure  1.  The  msin  tunnel  design  data  and  test 
equipment  are  summarized  in  figure  2.  Further  details  car  be  found  in  previous  papers  (Ref.  1,2  and  3). 
Considering  both  the  several  ambitious  goals  and  the  special  features  of  this  facility,  it  may  be 
concluded  that  the  very  combination  of  these  development.*,  the  unusually  wide  range  of  test 
capabilities  on  the  one  hand  and  the  high  quality  of  the  -test  data  on  the  other,  requires  special 
attention  for  correction  procedures  in  general  and  for  wg 11  interference  in  particular. 

It  is  well  known  that  the  greater  part  of  the  low-speed  wall  correction  methods  has  been  developed 
many  decades  ago  and  is  based  on  elementary  model  representations  and  on  linear  theory.  In  several 
respects  tnese  methods  have  not  kept  pace  with  modern  trends  in  low-speed  wind  tunnel  testing,  like 
higher  angles  of  attack  (large  wakes  and  separated  flow)  ,  higher  lift  (vortex  wake  relocation),  and 
more  complete  engine  simulation  (larger  jets) .  Thus,  it  may  be  anticipated  that  real  problems  may  arise 
when  in  addition  also  relatively  large  model  dimensions  are  applied,  attended  with  requirements  on  data 
accuracy  (in  particular  with  respect  to  drag  measurement)  which  are  significantly  higher  than  in  former 
days. 

In  search  of  means  which  could  help  to  cope  with  this  new  situation  and  which  could  be  useful  to 
attain  that  relatively  large  models  can  be  applied  without  making  concessions  to  the  accuracy  of  the  data, 
it  was  decided  that  all  DNW  test  sections  would  be  provided  with  longitudinal  slots  of  a  continuously 
variable  •width  (see  Fig.  3).  Obviously,  the  adoption  of  slotted  walls  can  not  be  considered  as  a 
substitute  for  necessary  renewals  and  improvements  of  wall  correction  theory.  On  the  contrary,  even  new 
problems  wise,  since  it  is  not  to  be  expected  that  an  interference-free  test  section  could  be  created 
in  this  way.  Therefore,  one  is  still  forced  to  develop  a  reliable  method  for  the  correction  of  the 
remaining  (reduced  and/or  more  simple)  wall  effects. 

All  in  all,  these  developments  have  given  rise  to  a  thorough  reorientation  in  the  field  of  low- 
speed  wall-interference  correction  procedures  for  both  slotted  and  closed  test  sections.  In  the  present 
paper  a  i wiew  is  given  of  recent,  current  and  planned  work  (to  be)  executed  at  NLR.  Some  of  this 
research,  intended  both  to  improve  wall  correction  methods  and  to  develop  new  ones,  is  carried  out 
within  tie  framework  of  existing  DFVLR-NLR  co-operation  programs. 

A  survey  of  the  main  items  of  the  low-speed  wall-interference  studies  is  given  in  figure  1*.  This 

survey  is  not  complete;  further  items  will  be  mentioned  i: i  the  text  of  the  present  paper.  Not  mentioned 

at  all,  neither  in  figure  U  nor  in  the  text,  is  a  separate  research  program  set  up  some  years  ago  for 
the  high-s  ieed  regime  and  intended  for  closed  and  slotted-wall  test  sections.  In  this  frame  a  new  two- 
dimensional  wall  correction  method  is  being  developed  which  is  based  on  measured  wall  pressures  and 
which  uses  a  model  representation  by  point  singularities .This  approach,  started  in  behalf  of  the 
transsonic  wind  tunnels  at  NLR  and  intended  to  be  extended  also  to  three-dimensional  test  conditions  in 
course  of  time,  concentrates  on  compressible  flow’  rather  than  on  typical  features  like  high  lift,  high 

angles  of  attack,  etc.  For  this  reason  this  method  is  out  of  order  now  and  will  not  be  discussed  here 

further.  A  description  of  this  method  is  given  in  reference  19. 


2.  •  PAST  AND  CURRENT  RESEARCH  ON  SLOTTED-WALL  INTERFERENCE 
2.1  The  design  philosophy  behind  the  DNW  slotted  test  sections 

Formerly, the  application  of  longitudinal  slots  in  test-section  walls  of  low-speed  wind  tunnels  was 
not  very  usual,  at  least  up  to  several  years  ago  when,  particularly  in  the  U.S.A.  ,  special  V/STOL  wind 
tunnels  of  large  size  were  built.  It  was  argued  then,  that  typical  V/STOL  aircraft  models  in  closed  test 
sections  can  create  so  large  and  unpredictable  flow  distortions  that  it  would  be  no  longer  possible  to 

derive  reliable  wall  corrections.  Based  on  the  simple  reasoning  that  open  and  closed  test  sections 

create  wall  effects  of  different  sign,  it  seemed  logical  to  assume  that  ventilated  walls  could  largely 
reduce  or  even  eliminate  at  least  an  important  part  of  these  wall  effects.  In  fact,  the  problem  was 
postponed  rather  than  solved,  for  it  was  highly  uncertain  which  slot  configuration  should  be  chosen  and 
to  what  extent  the  wall  effects  could  be  reduced  actually  under  tbe  widely  varying  test  conditions  which 
may  occur  for  V/STOL  models.  However,  uncorrected  wall  effects  of  unknown  but  certainly  reduced  magnitude 
often  were  preferred  to  large  wall  effects  which  should  be  corrected  by  highly  faulty  or  even  elearly 
failing  methods. 

The  design  philosophy  adopted  for  the  DNW  test  sections  (except  for  the  9.5  m  x  9.5  test  section) 

is  of  a  different  nature,  however.  In  fact  the  present  approach  is  a  middle  course  between  two 

possibilities:  the  development  of  a  very  complicated  analytical  correction  procedure  for  closed  test 
sections  on  the  one  side  and  the  self- correcting  wind  tunnel  with  adaptive  walls  on  the  other. 

It  was  realized,  that  important  progress  has  been  and  may  still  be  made  in  the  field  of  computational 
aerodynamics,  so  that  advanced  correction  methods  might  be  developed  in  the  near  future.  As  an  early 
example  the  nonlinear  method  of  Joppa  (Ref. 5)  may  be  mentioned,  in  which  the  wall-induced  vortex  wake 
relocation  is  taken  into  account  by  applying  a  vortex-lattice  method  in  an  interative  way.  Such  methods, 
which  are  indispensable  in  several  respects,  e.g.  in  assessing  the  limits  of  validity  of  "classical" 
pitching-moment  corrections,  do  not  seem  to  yield  the  most  practical  solution,  however,  as  they  are  too 
cumbersome  to  be  used  as  part  of  a  standard  correction  procedure  in  routine  wind  tunnel  testing. 

*)  The  DNW  has  been  built  and  will" be  managed  and  operated  by  the  DNW  Foundation,  a  new  organisation 
established  as  a  subsidiary  institute  by  the  German  DFVLR  and  the  Dutch  NLR. 
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Since,  in  addition,  for  DNW  consequent  on-line  data  processing  is  pursued,  such  that  all  data  will 
"become  available  as  much  as  possible  in  corrected  form,  it  was  considered  to  be  more  worthwhile  to 
attempt  whether  slotted  walls  could  be  applied  in  combination  with  relatively  simple  correction 
procedures . 

The  function  of  the  slotted  walls  in  the  DNW,  therefore,  is  limited:  they  are  intended  principally 
for  executing  some  c >ntrol  over  the  most  troublesome  wall  effects,  i.e.  the  variations  of  the  wall- 
induced  velocities  over  the  model.  If  such  an  attempt  would  turn  out  to  be  successful,  i.e.  if  the  wall- 
interference  gradients  could  be  reduced  sufficiently,  the  role  of  the  theoretical  corrections  method 
also  could  he  limited  to  derive  the  remaining  corrections  for  a  homogeneous  wall- interference  field  in 
a  reliable  way.  Besides,  some  special  intentions  are  present  for  tue  slotted  floor  in  the  low-speed 
9.5  m  x  9-5  m  test  section  *)  viz. to  shift  the  low-speed  test  limit  .which  stems  from  flow  breakdown  or 
recirculation  phenomena  in  V/STOL  model  testing  in  closed  test  sections. 

2.2  Numerical  studies 

As  it  was  pointed  out  above,  the  adoption  of  slotted  walls  in  the  interchangeable  test  sections 
of  the  DNW  was  considered  neither  as  a  final  solution  of  all  wall  interference  problems  nor  as  a 
substitute  for  necessary  renewals  and  improvements  of  wall  correction  theory.  On  the  contrary,  even 
two  new  problems  were  introduced  really: 

(i)  to  find  out  how  effective  slotted  walls  would  be  in  suppressing  the  inhomogeneity  of  wall- interfer¬ 
ence  over  the  model  and 

(ii)  to  develop  a  reliable  method  for  the  correction  of  the  remaining  wall  effects. 

In  order  to  explore  the  first  question  theoretically,  a  computer  program  was  developed  for  calculat¬ 
ing  lift  interference  in  rectangular  test  sections  with  ventilated  (slotted  of  perforated)  walls  of 
finite  length  (Ref.  6'.  This  program  is  based  on  a  three-dimensiona.1  extension  of  the  first-order 
(linearized,  incompressible)  version  of  the  method  described  by  Slooff  and  Piers  (Ref. 7).  It  proceeds 
from  a  source-panel  singularity  distribution  as  a  representation  of  the  tunnel  walls  and  it  uses  a 
modified  form  of  the  ’'classical"  linear  homogeneous  boundary  condition  due  to  Baldwin  et  al.(Ref.  11). 

The  essential  features,  of  this  method  are  summarized  in  figure  5. 

The  modification  of  this  analytical  boundary  condition  occurs  is  a  consequence  of  the  finite 
length  of  the  ventilated  part  of  any  actual  test  section,  and  consist,  of  the  introduction  of  the 
quantity  C  on  the  right  hand  side  which  is  proportional  to  the  statin-pressure  difference  between  the 
plenum  chamber  (adjacent  to  the  ventilated  wall)  and  the  undisturbed  flow  far  upstream.  In  this  respect 
the  NLR  approach  differs  from  that  of  reference  8  for  instance,  where  the  finite-length  aspect  is 
considered  separately  from  the  pressure  difference  between  plenum  ard  upstream  flow.  Even  more 
important,  and  also  different  from  reference  8,  is  the  fact  that  in  vhe  present  approach  C  is  shown  to 
belong  to  the  unknowns  of  the  problem  and  should  be  determined  there 'rom  after  adopting  an  additional 
equation.  This  equation  expresses  the  plausible  physical  fact  of  zero  net  mass  flow  through  the 
ventilated  wall(s).  However,  since  this  modification  does  not  change  the  basic  form  (the  left  hand  side) 
of  this  boundary  condition,  some  serious  objections  as  to  its  insufficient  theoretical  basis  and  its 
problematical  use  in  practice  remain  the  same  as  in  its  original  ‘cm. 

The  principal  objection  is  the  rather  intuitive  introduction  of  the  porosity  or  viscosity  parameter 
1/R  or  Q  «  (l  +  1/R)”1,  which  turns  out  to  have  a  very  large  influence  on  the  calculated  wall- interfer¬ 
ence  while  its  magnitude  can  not  be  predicted  from  the  actual  slotted  wall  configuration.  Nevertheless, 
this  type  of  boundary  condition  was  used  as  a  basis  because  (i)  it  was  the  only  available  one  at  the  time 
and  (ii)  it  was  attractive  because  of  its  linearity.**)  It  is  also  obvious,  that  for  the  limited  purpose 
of  this  investigation,  i.e.  an  assessment  of  the  effectiveness  of  wall  slots  in  reducing  upwash 
gradients  over  the  model, ‘somewhat  lower  requirements  may  be  accepted  than  when  it  was  intended  to 
derive  definitive  wall  corrections.  In  addition,  from  the  very  outset  a  rather  extensive  experimental 
investigation  was  planned  for  verifying  the  computational  results. 

A  parameter  study  has  been  carried  out  using  the  computer  program  of  reference  6  and  a  set  of  input 
data  as  summarized  in  figure  6.  This  survey  represents  only  the  main  part  of  this  study;  actually  also  a 
square  test  section  and  both  axial  and  vertical  variations  of  model  location  have  been  considered. 

Further  details  and  the  full  results  have  been  reported  in  reference  9  and  10;  a  summary  of  the  latter 
report  has  been  included  in  reference  2. 

The  principal  findings  are  shown  in  figure  7  through  10  and  may  be  summarized  as  follows: 

-  when  not  only  minimization  of  streamwise  but  also  of  spanwise  variations  in  wall-induced  upwash  is 
pursued,  all  four  walls  should  be  ventilated, 

-  a  kind  of  optimum  has  been  found  for  a  rectangular  test  section  with  perforated  walls  (iCj  v  «  0)  at 

Q„  =  0.66  and  =  0.72  (R^  =  1.9^5  and  Ry  =  2.565) »  for  in  that  case  0  and  6g(y)  «  constant  along 
the  lifting  line  of  an  unswept  wing  (see  Fig. 7);  this  result  was  flattering,  however,  since  the  absence 
of  streamline  curvature  occurs  only  at  the  lifting  line. 


*)  Strictly  speaking,  the  design  considerations  pointed  out  here  do  not  apply  at  all  to  the  9.5  m  x 
9.5  m  test  section;  in  this  test  section  only  the  floor  has  been  provided  with  slots,  whereas  in 
the  other  test  section*  slots  are  present  in  all  four  walls . 

**TAs  has  been  indicated  in  figure  the  theoretical  work  at  NLR  directed  on  a  more  firmly  founded 
analytical  boundary  condition  is  still  continuing.  Presently  it  has  no  high  priority,  however, 
although  it  is  realized  that  it  is  of  interest  for  the  design  of  slotted  walls. 
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-  in  order  to  reduce  streamline  curvature  also  between  wing  and  tail,  and  thus  to  minimize  the  pitching- 
mameiit  correction  for  a  complete  model  with  horizontal  tail,  a  slotted  test  section  is  to  he  preferred 
to  a  perforated  one, 

-  most  favourable  with  respect  to  longitudinal  variations  in  lift  interference  is  found  to  he  a 

(rectangular)  test  section  with  almost  ideal  slots  1  or  1/Rw  0)  of  small  width  and/or  large 

spacing  (K  large),  (see  figure  6) , 

-  reductions  of  both  spanwise  and  streamwise  variations  may  be  achieved  also  for  smaller  and  possibly 
even  more  realistic  values  of  Q  and  K  (see  Fig. 9). 

As  a  tentative  conclusion  from  this  exploratory  investigation  it  may  be  stated  therefore,  that,  an 
actual  slotted  test  section  is  conceivable  which  induces  only  homogeneous,  and  thus  easily  correctable, 
lift  interference.  It  should  be  noted,  that  this  result  was  obtained  for  a  relatively  large  wing  model 
(span  «  0.6  tunnel  width)  as  well  as  for  a  rotor  model  (diameter  «  0.5  tunnel  width)  at  not  too  large 
lift  coefficients  (see  Fig. 19).  The  effect  of  slot  length  is  small  at  the  lifting  line,  and  remains 
within  bounds  elsewhere  when  this  length  is  not  small  el*  than  two  times  the  test  section  width,  at  least 
for  not  too  long  models. 

2.3  Experiments 

Especially  in  the  present  case,  where  the  available  computational  method  was  imperfect  because  of  a 
disputable  boundary  condition  applied  for  a  slotted  wall,  a  thorough  experimental  verification  was  felt 
to  be  indispensable  .  In  the  framework  of  existing  DFVLR-NLR  co-operation  programs  an  extensive 
experimental  investigation  was  set  up  to  be  executed  in  a  modeltunnel  (MUK)  of  the  formerly  planned 
German  large  low-speed  windtunnel  (GUK)  at  DFVLR-Braunschweig. 

Ski  tches  of  both  model  and  test  section  are  given  in  figure  11  snd  12.  The  model  geometry  was  chosen 
as  simp?  e  as  possible;  lift,  drag  and  pitching  moment  are  measured  by  an  internal  strain  gauge  balance. 
Initially  it  was  intended  to  investigate  a  large  variety  of  slotted-wall  configurations  (several  slot 
widths,  slot  lengths,  etc.),  but  up  to  now  only  the  configuration  shown  in  figure  12  has  been  measured. 

The  tent  section  has  a  square  cross  section,  its  open-area  ratio  is  12  %  and  the  slot  length  is  only 

1.3  times  the  test-section  width.  As  is  shown  in  the  sketch  in  figure  12,  the  test  section  was  surrounded 
by  a  plenum  chamber  of  rather  small  depth.  The  plenum  vails  can  also  be  removed,  however. 

T:p  to  now  a  complete  set  of  three-component  measurements,  i.e.  polar  runs  at  seven  different  tai. 
setii.  gs  and  without  tail  surfaces,  has  been  carried  out  three  times,  viz.  in  the  closed  test  section,  in 
the  .‘l.'tted  test  section  with  and  in  the  slotted  test  section  without  plenum  chamber  *) .  In  preceding 
calibration  tests  it  was  checked  that  the  empty  test-s<  ction  flow  did  not  change  significantly  with  wall 
configmation,  so  that  differences  in  results  indeed  could  be  ascribed  to  differences  in  wall  interference. 
Much  attention  was  paid  to  differences  in  pitching  moment,  which  are  directly  connected  with  streamline 
curvatuis,  thus  with  differences  in  wall-induced  upwash  between  wing  and  tail. 

An  example  of  the  lift  measurements  without  fail  is  given  in  figure  13.  The  most  remarkable  item  is  a 
rather  small  difference  between  the  uncorrected  results  for  the  closed  test  section  and  the  slotted  test 
section  with  plenum  chamber.  This  means  that  the  effectiveness  of  the  slots  is  very  small,  and  that  r 
positive  angle-of-attack  correction  will  exist  of  more  than  one  half  of  the  correction  of  the  closed  test 
section  (the  corrected  lift  curve  for  the  closed  test  section  is  indicated  by  flagged  circles  in 
figure  13). 

Presumably  this  feature  is  caused  by  the  small  volume  of  the  plenum  chamber.  When  the  slotted  test 
section  with  plenum  is  considered  as  a  closed  test  section  with  inner  dimension* equal  to  those 
of  the  plenum  chamber,  and  the  measured  data  are  corrected  in  tbe  usual  way  for  lift  interference,  a 
corrected  lift  curve  is  found  (indicated  by  flagged  squares  in  figure  13)  which  coincides  almost  perfectly 
with  the  former.  Although  the  results  for  the  slotted  test  section  without  plenum  are  not  yet  available, 
it  may  be  concluded  that  the  small  plenum  chamber  will  be  the  principal  cause  of  the  poor  slot  effectiv- 


For  the  pitching-moment  results  a  similar  correction  procedure  has  been  attempted  (Fig.l4),  but  in 
this  case  the  simple  interpretation  of  the  slotted  test  section  with  plenum  as  a  closed  test  section  with 
the  dimensions  of  the  plenum  is  not  as  convincing  as  in  the  former  case.  In  fact  the  uncorrected  results 
are  so  nearly  coinciding  that  any  correction  which  differs  between  the  closed  and  the  slotted  test  section  will 
disperse  the  curves.  On  the  other  hand,  if  from  the  difference  in  angle-of-attack  correction  between  the 
slotted  and  the  closed  test  section  a  value  of  Q  is  determined  (ataknovn  value  of  K)  so  that  the 
pitching-moment  correction  can  be  calculated,  a  very  good  agreement  is  found  between  the  measured  and 
calculated  pitching-moment  correction.  , 

As  the  results  of  the  tests  without  plenum  are  not  yet  available,  the  issue  of  this  investigation 
remains  undecided.  As  a  tentative  conclusion  it  may  be  stated,  that  the  slotted  test  section  surrounded 
by  this  small  plenum  chamber  may  have  some  characteristics  in  common  with  both  a  slotted  test  section 
of  small  effectiveness  (small  Q)  and  a  closed  test  section  of  enlarged  dimensions.  The  small  slot  length 
may  be  a  secondary  cause  of  the  small  effectiveness  of  the  slots ,  especially  concerning  the  pitching- 
moment  correction. 

3  SOME  SOLID-WALL  INTERFERENCE  STUDIES 

3.1  Two-dimensional  high  angle  of  attack  models  (attached  flow) 

Although  it  is  not  exactly  of  recent  date,  for  the  sake  of  completeness  the  work  of  de  Vries  and 
Schipholt  may  be  recalled  briefly  here  (Ref.  12).  It  was  argued  several  years  ago,  that  both  improvements 

*)  The  latter  measurements  have  not  yet  been  completed. 
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of  two-dimensional  testing  methods  and  developments  in  the  direction  of  more  advanced  high-lift  devices 
with  high  incidences  and  flap  deflections  necessitated  a  check  on  the  validity  of  the  commonly  used 
linearized  corrections  due  to  Glauert.  This  was  motivated  also  by  the  results  of  earlier  calcinations 
due  to  De  Jager  and  Van  de  Vooren  (Ref.  13),  from  which  strongly  nonlinear  effects  had  been  found  for  an 
infinitely  thin  hinged  flat  plate  at  zero  incidence. 

The  calculations  of  De  Vries  and  Schipholt  (Ref.  12)  were  restricted  to  the  inviscid  wall-interfer¬ 
ence  effects,  thus  excluding  wake  blockage  and  boundary  layer  effects,  and  were  carried  out  by  means 
of  an  adapted  version  of  a  computer  program  for  multi-element  aerofoils.  This  program  in  turn  was  based 
on  the  NLR  panel  method.  It  proceeds  from  source  panels  on  both  the  aerofoil  contour  and  the  tunnel  walls 
and  a  prescribed  vortex  distribution  on  the  mean  line  of  unknown  overall  strength. 

The  results  substantiated  the  suppositions  in  that: 

-  the  correction  theory  of  Glauert  breaks  down  for  large  flap  deflections  and  may  be  used  for  single 
aerofoils  only  at  not  too  large  angles  of  attack. 

-  the  results  of  De  Jager  and  Van  de  Vooren  cannot  be  used  to  estimate  the  pitching-moment  correction 
for  a  flapped  aerofoil  with  slat. 

These  results  are  very  similar  to  those  obtained  independently  by  Eriksson  (Ref.  1^),  who  states  generally 
that  the  very  displacement  of  the  aerofoil  (elements)  in  vertical  direction  causes  the  wall  corrections 
to  be  nonlinear  in  CT . 

h 

The  final  conclusion  of  reference  12  has  been  summarized  schematically  in  figure  15.  It  implies  that 
wall  corrections  for  multielement  aerofoils  can  be  estimated  from  a  combination  of  two  types  of  calculat¬ 
ions,  i.e.  lift  interference  from  a  new  non-planar  vortex  lattice  method  (thus  accounting  for  camber  but 
neglecting  thickness  effects)  and  blockage  from  conventional  techniques  (thus  including  thickness  and  wake 
effects).  This  approach  has  in  common  with  classical  approximations  that  lift  interference  and  blockage 
are  considered  to  be  independent  of  each  other. 

Intentions  to  derive  a  practical  wall-correction  procedure  alorg' these  lines  have  not  been  carried 
out,  however,  and  that  for  several  reasons.  The  main  reason  is  that  such  an  isolated  perfection  of  the 
two-dimensional  wall  corrections  received  less  attention  in  course  of  time.  In  the  mean  time  these  plans 
have  been  superseded  by  expectations  which  maybe  taken  from  a  more  integrated  approach:  a  two-dimensional  version 
of  a  new  method  based  on  measured  tunnel-wall  pressure  distributions  to  be  discussed  in  section 

Wall-interference  effects  for  two-dimensional  tests  in  presence  of  strong  viscous  effects,  either 
on  the  model  itself  or  at  the  tunnel  walls,  have  not  yet  received  special  attention  at  NLR.  However, 
when  flow  separation  occurs  only  on  the  model,  it  can  be  considered  as  a  special  case  of  the  next  item. 

3.2  Models  with  large  regions  of  separated  flow 

Especially  with  a  view  to  very  high  angle  of  attack  conditions  (deepstall  measurements)  of  complete 
models  of  transport  aircraft,  as  well  as  for  high-lift  conditions  in  general,  a  need  was  felt  at  NLR  of 
a  method  by  which  additional  blockage  effects  due  to  separation  bubbles  could  be  estimated. 

In  a  preliminary  investigation  (Ref.  15)  the  method  due  to  MasKell  (Ref.  1 6),  in  a  form  as  recommen¬ 
ded  by  Vayssaire  (Ref.  IT),  was  judged  to  be  a  tool  which  could  meet  a  direct  need.  An  important  feature 
of  this  method  is  its  suitability  to  be  used  in  a  fully  computerized  on-line  data  processing  system. 

Last  year  this  method  was  incorporated  into  the  general  low-speed  windtunnel  data  reduction  program 

(APVW)  of  NLR  by  way  of  trial.  This  new  part  of  the  program  is  being  tested  now  in  the  usual  standard 
reduction  process  of  representative  wind  tunnel  measurements.  A  schematic  picture  of  the  main  features 
of  this  procedure  is  given  in  figure  1 6.  Its  on-line  character  is  due  to  the  fact  that  the  procedure 
already  can  be  started  before  the  polar  run  is  finished.  For,  as  soon  as  the  ideal  polar  has  been  derived 

from  measuring  points  in  a  pre-selected  low  angle  of  attack  range,  the  correction  procedure  can  be 

applied  immediately  to  any  following  measuring  point  beyond  this  range. 

The  experience  gained  thus  far  at  NLR  with  this  method  is  not  unsatisfactory,  although  some  short¬ 
comings  declare  themselves  clearly.  A  practical  imperfection  is  found  in  the  case  of  model  configurations 
with  deflected  slats,  where  flow  separations  occur  even  in  the  low-angle-of-attack  portion  of  the  polar. 

In  these  eases  the  shape  of  the  ideal  polar  can  be  affected  to  such  an  extent  that  very  irrealistic  values 
of  A',  and/or  Ct  '  are  found.  But  even  then  the  dynamic -pres sure  correction  often  turns  out  to  be  of  a 
credible  magnitude.  A  more  fundamental  shortcoming  is  the  fact,  that  variations  over  the  model  of  this 
kind  of  wall  interference  can  not  be  accounted  for,  and  thus  additional  buoyancy  corrections  arc  neglected. 
This  means  that  only  the  most  striking  effects  are  considered  and  that  in  particular  the  drag  correction 
may  be  not  very  accurate.  Also  in  this  context  expectations  are  set  on  a  method  based  on  measured  wall 
pressure  distributions  as  has  been  developed  at  Lockheed-Georgia  (Ref.  18). 

k  NEW  METHODS  USING  MEASURED  FLOW  QUANTITIES  NEAR  THE  WALLS 

From  the  previous  sections  it  becomes  evident ,  that  in  the  theory  of  low-speed  wall  corrections 
in  slotted  and  closed  test  sections,  principally  two  types  of  problems  exist  which  can  hardly  be  solved 
satisfactorily  by  conventional  or  classical  means.  Firstly,  for  slotted  test  sections  a  key  problem  is 
the  lack  of  an  adequate  boundary  condition  for  a  slotted  wall.  An  analytical  approach  of  this  subject 
would  supply  a  need  certainly,  particularly  in  behalf  of  test-section  design,  but  a  solution  which  would 
be  satisfactory  also  from  a  practical  point  of  view  can  not  be  expected  to  be  reached  at  short  notice. 
Secondly,  in  closed  as  well  as  in  ventilated  test  sections  problems  ‘arise  as  to  the  theoretical  model 
representation.  These  problems  are  aggravated  by  several  typical  features  of  modern  aircraft  designs  and 
test  requirements  (e.g.  large  models,  separated  flow,  large  wakes  and  jets,  etc.). 
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To  some  extent  both  types  of  problems  can  be  discerned  also  in  transsonic  wind  tunnel  testing, 
although  the  origin  differs  of  course,  and  it  is  not  fortuitous  then  that  the  solution  is  sought  in  a 
similar  way.  At  NLR,  as  well  as  in  many  other  places  (e.g.  Ref.  *0,  it  is  anticipated  that  the  problems 
mentioned  ash  for  an  approach  which  uses  measured  flow  quantities  in  the  vicinity  of  the  tunnel  walls 
obtained  during  the  test  in  addition  to  measured  model  forces  and/or  pressure  distributions.  The  use 
of  measured  wall  flow  quantities  (e.g.  static  pressures)  overcomes  the  need  for  an  analytical  boundary 
condition,  as  has  been  pointed  out  in  reference  4,  but  this  experimental  information  may  be  used  also  to 
detect  possible  uncorrigible  conditions.  The  same  information,  possibly  in  combination  with  (extra) 
measured  pressures  on  the  model,  may  be  used  to  overcome  problems  with  respect  to  the  model  representa¬ 
tion.  In  this  way  more  certainty  will  exist  about  the  actual  flow  phenomena  during  the  test,  and 
appropriate  and  effective  means  may  be  taken  at  once  (e.g.  by  limited  wall  control)  to  such  an  extent 
that  corrigible  conditions  arise  again. 

Returning  to  figure  4  it  is  worth  mentioning,  that  the  vast  part  of  the  work  summarized  under  the 
heading  "Measured  boundary  conditions"  should  still  be  d^ne,  and  that  it  can  be  considered  to  be  of  a 
long-term  type.  Since  the  activities  in  this  field  have  been  started  only  recently,  results  can  not 
yet  be  shown. 

Firstly,  possibilities  are  being  considered  presently  of  developing  a  new  universal  method  using 
both  measured  wall  pressures  and  model  forces,  which  can  dispense  with  a  theoretical  model  representa¬ 
tion.  It  is  intended  to  be  applied  to  both  closed  and  slotted  walls  and  includes  both  lift  interference 
and  blockage.  The  basic  idea  is  due  to  mr.  J.W.  Slooff  of  NLR  and  has  some  aspects  in  common  with  the 
method  ol  Murman  described  in  reference  4.  Just  like  in  the  "self-correcting"  or  "adaptive-wall"  concept 
a  relatively  simple  fictitious  potential  flow  exterior  to  the  tunnel  is  calculated  from  a  measured 
distribution  of  the  horizontal  velocity  component  along  the  walls.  From  this  calculation  a  distribution 
of  the  vertical  velocity  component  is  obtained  which  generally  differs  from  the  measured  one.  This 
difference  is  a  measure  of  the  wall  interference.  Otherwise  than  in  a  self-correcting  tunnel,  however, 
it  is  no.  the  actually  measured  normal-velocity  distribution  which  is  changed  by  some  form  of  active 
wall  control  (by  blowing,  suction  or  wall  deformation),  but  the  calculated  distribution  is  influenced  by 
modifyirg  the  boundary  conditions  at  infinity,  i.e.  by  varying  the  horizontal  velocity  component  Uw  and 
the  vertical  component  V  .  These  components  are  determined  such  that  both  the  tangential  and  normal 
velocity  at  the  walls  are  matched  in  some  mean  sense  (i  e.  similar  to  Murman  in  matching  normal  and 
tangential  velocities  at  the  airfoil  surface,  see  reference  4).  The  components  and  Vm  thus  found  are 
the  sought  corrections  on  magnitude  and  direction  of  tVe  tunnel  flow  respectively,  corresponding  with 
blockage  and  downwash  (lift  interference). 

Thus,  as  described  so  far,  this  method  would  require  both  pressure  measurements  (to  determine  the 
tangential-velocity  distribution)  and  flow  direction  measurements  (to  determine  the  normal -velocity 
distribution)  in  the  vicinity  of  the  walls. .The  measurem -nt  of  the  latter  does  not  seem  to  be  attractive 
in  practice.  However,  this  difficulty  can  be  avoided  in  principle  by  utilizing  relations  between  the  wall 
velocity  distributions  and  the  measured  model  forces  which  can  be  derived  from  momentum  considerations 
(in  two-d.mensional  flow  equivalent  to  Blasius’  theorems).  In  this  way,  the  matching  procedure  could  be 
carried  out  through  longitudinal  wall  velocities  and  model  forces,  rather  than  normal  wall  velocities*). 
Presently  the  feasibility  study  concentrates  on  practical  and  fundamental  aspects  of  this  matching 
procedure,  as  it  is  felt  that  the  viability  of  this  method  will  depend  crucially  on  such  points. 

An  alternative  approach  is  being  taken  into  consideration  next,  intended  especially  for  complicated 
models  in  closed  test  sections.  It  proceeds  from  the  method  developed  at  Lockheed-Georgia  by  Hackett  et 
al.  (Ref.  18)  for  estimating  all  kinds  of  blockage  corrections,  viz.  model  volume,  wake  and  jet  effects 
and  separation  bubbles.  Typical  for  the  Lockheed-met hod  is  that  wall  pressures  are  used  to  define  the 
theoretical  model  representation.  The  calculation  of  the  wall  interference  is  similar  to  classical 
methods.  It  is  intended  at  NLR,  however,  to  include  lift  interference  right  from  the  onset  and  that  for 
high-lift  models  with  or  without  jet  simulation.  As  a  first  step  a  combined  numerical  and  experimental 
investigation  has  been  started,  preliminarily  restricted  to  two-dimensional  cases.  It  may  be  expected 
that  this  method  will  become  operational  at  relatively  short  notice  and  that  it  will  be  able  to  replace 
several  methods  for  special  purposes,  like  those  for  multi-element  aerofoils,  separated  flows,  etc. 

(see  section  *3,).  Dependent  on  the  results  of  the  present  preliminary  studies  it  will  possibly  be  decided 
to  stop  the  further  development  of  either  the  first  or  the  second  method,  so  that  only  a  single  method 
might  be  considered  for  extension  to  three-dimensional  flow. 
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5  CONCLUDING  REMARKS 

For  several  years  it  is  realized  at  NLR  that,  some  of  the  commonly  used  wall  correct  ion  methods  for  closed 
test  sections  should  be  re-examined  thoroughly  before  they  could  be  applied  satisfactorily  to  modern 
low-speed  wind  tunnel  testing.  Although  for  low  speeds  the  development  of  methods  based  bn  purely  computa¬ 
tional,  means  seems  to  be  conceivable,  it  may  not  be  expected  to  result  in  an  optimum  solution. 

For  the  new  German-Dutch  Wind  Tunnel  (DNW)  slotted  walls  are  chosen  as  a  means  to  minimize  the  wall 
interference  gradients  over  the  model.  In  order  to  derive  reliable  corrections  for  the  remaining  homoge¬ 
neous  wall  interference,  methods  are  considered  which  make  use  of  measured  flow  quantities  near  the  walls, 
although  also  work  on  analytical  boundary  conditions  is  continuing.  Besides  this  long-term  research,  also 
attention  is  payed  to  short-term  solutions  for  closed  test  sections.  In  this  context  expectations  are  set 
on  a  Lockheed-type  of  approach  as  a  means  for  defining  the  necessary  theoretical  model  representation. 


*)  It  should  be  noticed  that  the  procedure  described  here  applies  to  ventilated  vails  in  general.  In  the 
special  case  of  a  closed  wall  the  situation  is  more  simple  because  the  normal  velocity  component  equals 
zero  at  the  wall. 
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PREVAILING  DESIGN  PRINCIPLES  INVOLVED 

•  HIGH  AERODYNAMIC  ANO  AERD-ACOUSTIC  QUALITIES 

•  COMPREHENSIVE  AND  AOVANCEO  TEST  EQUIPMENT 

•  LARGE  VARIETY  OF  MODELS  ANO  TEST  CONDITIONS 

•  HIGH  TESTING  PRODUCTIVITY 

•  FLEXIBLE  AND  ECONOMIC  OPERATION 

MAIN  ACTIVITIES  FOCUSED  ON  : 

-  IMPROVEMENT  OF  AIRCRAFT  LOW-SPEED  CHARACTERISTICS 
(SAFETY,  PERFORMANCE  ANO  ECONOMY ) 

--  HIGH- LIFT  OEVICES 

-  V/STOL  AERDOYNAMICS 

-  ENGINE/AIRFRAME  INTERFERENCE 

-  AIRFRAME  ANO  ENGINE  NOISE 

-  ROTOR  AEROOYNAMICS 

-  (HIGH-SPEEO )  NE LICDPTERS 

-  FLUtTER  TESTS 

-  JETTISON  TESTS 

-  OPTIMIZATION  OF  FULL-SCALE  AIRCRAFT  COMPONENTS 

-  REAL  ENGINES  (INTAKE,  JET  PLUME  EFFECTS) 

-  NON-AERDNAUTICAL  INVESTIGATIONS 
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Fig.  1  Design  principles  and  expected  activities 
of  the  DNW  facility 


Fig.  3  Upstream  view  of  the  convertible 

8m  x  6m  test  section  with  slotted  vails 


TYPE  OF  TUNNEL 

SIZE  OF  WORKING  SECTION 
TYPE  OF  SECTION 

CONTRACTION  RATIO 
MAX.  SPEED  (- A) 

STATIC  PRESSURE  IH  TEST  SECTION 
REYHOLOS  NUMBER  *  10“6*) 

MAIN  DRIVE 
AUXILIARY  ORIVES 
FAN 


CLDSEO  RETURN  CIRCUIT 
(OVERALL  LENGTN  OF  CENTERLINE: 
320  m) _ 

9.S  m  «  9.5  mf"  a..U  I  6  m.  6  n, 


9.S  aiiS* 
CLOSED*) 


CLOSED*)  CLOSED*) 
ANO  OPEN 


•  |  C»o)  |  j 

ATMOSPHERIC  (I  BAR  ) 

3.9  1  5-2  |  5.8 

THYR.  SYHCHR.  MOTOR;  NORMAL  RATING: 
MAINLY  FOR  COMPR,  AIR!  5:7  MW  12  7  MW 

SINGLE  STAGE;  8  BLADES;  DIRECT  DRIVE 
22S  PPM;  CONST.  PITCH,  WIHD  SPEEO 
CONTROL  BY  MOTOR 


)  BASED  OH  Y*,,  AHO  0.1  VA  (A:  T  EST  SECTION  AREA  ) 

*)  ALL  CLOSED  TEST  SECTIONS  ARE  PROVIDED  WITH  LONGITUDINAL  WALL  SLOTS 
OF  A  COHTIHUOUSLY  VARIABLE  WIDTH  RANGING  FROM  0  TO  0.12  m  AT  I  m  PITCH 


STANOARO  EQUIPMENT 

•  MOOEL  STING  SUPPORT  (VERTICAL  MOVEMENT  WITH  5mA  ) 

•  MECHAHICALSIX-COMPOHEHT  BALAHCE 

•  HEAT  EXCHAHGER  (IN  THE  SETTLIHG  CHAMBER) 

•  CIRCUIT  AIR  EXCHANGE  (OUTLET  &  IHLET  HATCHES.  THROTTLE  1 

•  CONTROL,  DATA  AOUISITION  AHD  REDUCTION  SYSTEM 

(DISTRIBUTED  COMPUTER  SYSTEM) _ 

SPECIAL  EQUIPMENT 


•  COMPRESSED  AIR  PLAHT 

CONT..  35 Vq/%  IHTERMITT 

WITH  AIR  STORAGE 

SO  m3.  280  bo. 

•  NYORAULIC  ORIVES 

•  NYOROCEH  PEROXIOE  PLANT 

FOR  HOT  GAS  SIMULATION 

•  MDVIHG  BELT  GROUND  PLANE 

57m  ,  6  m,  60mA 

( IHCL  BLC) 

•  Q-STOPPER 

RETRO-JET;  SAFETY  OEVICE 

FOR  FLUTTER  TESTS 

•  SCOOP 

REMOVAL  OF  HOT  GAS  FROM 

TEST  SECTIONS  BY  INJECTOR 

|  UPE  OF  RESEARCH  | 

OBJECTIVE 

PROCEOURE 

EFFECTIVENESS 

OF 

SLOTTEO 

WALLS 

NUMERICAL 

STUOIES 

AHO 

EXPERIMENTAL 

VERIFICATION 

OEVELOFMENT 

AND 

EVALUATION 
OF  A 

HEW  METHOD 

^FEASIBILITY 

STUOY 

2)  2-O.OEVELOP' 
MEHT 

3)  J-D. 

EXVEHSIOH 

POSSIBLV 

LATER  ON 

EVALUATIOH 
OF  A 

METHOD  FOR 
COMPLICATED 
MOO ELS 
(HIGH  UFT  ANO 

OR  LARGE  JETS, 
WAKES.  INCL.  SE 
PARATEO  FLOW) 

‘LOCKHEED-TYPE' 
OF  APPROACH. 
(REF.  It) 
EXTENDED 

TO  LIFT 
INTERFERENCE 

J-O. SLOTTED 
TEST  SECTIONS. 
ATTACNEO 
FLOW 


REMARKS 

(PROSPECTS) 

1 1)  STILL  IN 
PROGRESS 
(THEORETICAL 
STUDIES  AT 
LOW  PRIORITY) 
2JEXPERIMENTS 
FORM  PART  OF 
DFVLR-NLR 
CO-OPERATION 


AND  CLOSED 
TEST  SECTIONS 
(DNW.  HST  ) 


CLOSED  TEST  | 
SECTIONS,  I 


TO  SLOTTED  IF  METHOD (A) 
EST  SECTIONS  FAILS 
LATER  ON 


Fig,  2  Main  data  and  test  equipment  of  the  DNW 
facility 


U  Survey  of  principal  low-speed  wall 
interference  studies  at  NLR 
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*  or  * 


STARTING-POINTS;  POTENTIAL  PLOW  (IHCOMPR  }  A(p-Q  (LAPLACE  EOUATION ) 

9m9m  »NERE  DUE  TO  MODE L  ( REPRESENTED  BY 

A  SINGLE  SWEPT  HORSESHOE  VDRTEX  OF  KNOWN 
STRENGTNlv*  OUE  TO  WALLS  (REPRESENTED  BY 

source  Panels  of  unknown  strengths®) 


BOUHOARY  CONDITIONS; 


ON  CLOSED  parts  of  WALL 


0  ON  CLOSED  Pj 
O"  , 

dp  1  dyt  dp 

dx  K  dn  dxdn 


=  C  ON  VENTILATED  PARTS  OF  WALL 


2C 

WHERE  C  IS  RELATEOTOpe  :  — 

Ucc 


Pg.-Pc 

VipUan 


K  IS  A  FUNCTION  OF  SLOT  GEOMETRY  (SLOT  WIDTH, 
SLOT  S(ACIHG) 


1  R«(1~0)'0  IS  AN  EMPIRICAL  QUANTITY  (INTENDED 
TO  ACCOUNT  FOR  VISCOSITY  EFFECTS) 


NUMERICAL  SOLUTION;  BECAUSE  pm  IS  KNOWN.  TNE  UNKNOWN  SOURCE  STRENGTHS  ® 
AS  WELL  AS  TNE  UNKNOWN  OUANTITY  C  CAN  BE  SOLVEO 
FROM  I) AND  (2) 


*9w  1  k°9»  .  ,  *9m  1  dgm  l 

dx  H  dn  hr  in  dx  R  dn  dxdn  * 

-  -  -  (VENTILATED  PART) 


VENTIL  TED  PARTS 


0  (ZERO  NETMASS  FLOW 

TNROUCN  VENTILATED  WALL) 


Fig.  5  Outline  of  the  NLR  calculation  method  for  Fig.  7  Example  of  lift-interference  reduction  in 

lift  interftrence  in  ventilated  test  a  perforated  test  section, h/b  =  0.75* 

sections  *  2s, b  =  0.7  (representing  an  unswept  wing 

vita  span  s»  0.8b) 


TEST-SECTION  AND  WALL  CONFIGURATION 

CONSTANT 

VARIABLE 

HEICNT /WIOTN;  HA  *=0.7$ 

* 

SLOT  LENGTN/T.S.  WIOTH  LA  -2  TD  3 

★ 

SLOT  PARAMETER  KN  (HOR.  WALLS). 

Ky  (VERT.  WALLS) 

★ 

POROSITY  (OR  VISCOSITY  )  PARAMETER  QN,  0V 

* 

MODEL  CONFIGURATION  AND  SIZE 

UHSWEPT  WING  (AHO  HELICOPTER  ROTOR) 
REPRESENTED  BY  A  SINGLE  HORSESHOE 

VORTEX  OF  STRENGTH  T  ANO  SPAN  2.  t 

WIHC  MODEL  SPAN3;0,8!>  -2»/b=0,7 

★ 

ROTOR  OIAMETER  2R-S/R=:0.«7 

★ 

dx  JJ  dxd*  * 


Fig.  6  Survey  of  model  and  tunnel  data  used  in 
slotted  wall  study 


Fig,  8  Distribution  of  lift  interference  factor 
6q  ‘along  the  centre  line  of  a  closed, 

perforated  and  slotted  test  section, 
h/b  =  0.75;  2s/b  =0.7  (unswept  wing  with 
span  »  0.8b) 


*  A 

0.1  0.2  0.3  0.4  0.5  0.6  0.7 


Fig.  9  Distributions  of  60 along  test-section 
centre  line  and  ving  span  for  a 
rectangular  test  section. h/b  *  0.75; 

2s /b  ~  0.7  (unsvept  ving  vith  span  »  0.8b) 


(CALI  1  > 10 
01UEKSIONS  !H 


Fig.  11  Geometry  of  model  and  sting  used  in 
slotted-vall  interference  st  idr 
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Fig,  10  Distribution  of  angle  of  attack  correc-  Fig.  12  Test-section  vail  geometry  of  model 

tion  along  the  longitudinal  axis  of  a  tunnel  at  DFVLR-Braunschveig 

rotor  model  of  hm  diameter  in  a  slotted 
8m  x  6m  test  section, h/b  *  0.75; 

_  .  .  n  7,  . 
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VORTEX-LATTICE  ON  WALLS  ANO 
MEAN  LINES  OF  AEROFOIL  ELEMENTS 


+ 


SINGLE  SOURCE  (+SINK  AT  INFINITY) 
REFRESE  NT  MG  VISCOUS  WAKE 


Fig.  13  Influence  o:'  the  vail  configuration  and 
of  the  correction  procedure  on  (a). 

Model  without  horizontal  tail. 

Re  »  0.68  x  106 
c 


Fig.  15  S  heme  of  an  intended  (not  realized) 

ci  rrection  procedure  for  two-dimensional 
&.U- ti -element  aerofoils  in  closed  test 
section 


4q/q  =  e~  CDs,  WHERE  2.8-0.068  A* 

S  -REFERENCE  (WING) AREA 
C  =  CROSS -SECTIONAL  AREA  OF  T«T  SECTION 


c0m  CD 


(CL-C,  )j 

•IDEAL* POLAR  tC/).j~CD  +  - 

*  U‘d  m  7t  A 


a0-aIcL+a2Cl‘ , 


WHERE  Aq-Co^  +  - 


ARE  CALCULATEO  FROM  A  LEAST  SOUARES  METHOD  APPLIEO  FOR  (ABOUT  7) 
POINTS  ON  A  PRE-SELECTEO  LOW  ANCLE-OF-ATTACK  PORTION  OF  THE 
MEASUREO  POLAR  CURVE  (I.E.  A  PART  WITHOUT  FLOW  SEPARATION) 

TNEH  Cj^,  C^  ANO  A'  BEING  KNOWN,  Cp*  FOLLOWS  FROM 
„  *  (CL-C.  )3 

co»"*cd  -c*,d~cv-cOm - Z7, — 


Fig.  lU  Influence  of  wall  configuration  and 
correction  procedure  on  C  (CT). 

m  It 

Model  with  horizontal  tail  (i^  -  0°), 

Re  =  0.68  X  106 
c 


Fig.  16  Outline  of  the  Maskell/Vayssaire 

correction  procedure  for  separation 
blockage  on  three-dimensional  (wing) 
models 


6-1 


A  REVIEW  OF  SOME  INVESTIGATIONS  ON  WIND  .TUNNEL  WALL 
INTERFERENCE  CARRIED  OUT  IN  SWEDEN  IN  RECENT  YEARS 

by 

S.-E.  Nyberg 

FFA,  The  Aeronautical  Research  Institute  of  Sweden  • 
P.O.  Box  11021,  S-161  11  BROMMA,  Sweden 


SUMMARY 

A  survey  is  presented  of  wind  tunnel  wall  interference  studies  in  Sweden  relating  to 
testing  of  models  at  high  angles  of  attack,  published  after  1974  or  in  progress.  For 
subsonic  incompressible  flow  the  mutual  circulation-induced  model  -  wind  tunnel  interfer¬ 
ence  has  been  calculated  by  panel  methods  for  large  multicomponent  two-dimensional  air¬ 
foils,  for  three- limensional  swept  wings,  full  or  half  models,  and  for  wing-tail  con¬ 
figurations.  Wake  blockage  effects  from  a  swept  wing  with  and  without  high  lift  devices 
have  been  studied  experimentally.  The  effects  of  air  flow  leakage  between  half  model 
fuselage  and  reflection  wall  have  been  investigated.  For  transonic  flow  the  flow  proper¬ 
ties  of  slotted  vails  and  the  influence  of  wall  boundary  layer  have  been  studied.  Based 
on  these  results  a  numerical  method  has  been  developed,  and  so  far  axisymmetric  calcula¬ 
tions  have  beer,  carried  out.  The  results  have  been  compared  with  experimental  results 
for  large  bloc)  a  ie  models.  The  interference-free  flow  field  around  a  lifting  wing-body 
model  has  been  investigated  to  establish  cross-flow  characteristics  for  ventilated  walls. 
The  interference  on  a  half  model  from  a  full  model  support  strut  has  been  determined. 


Mach  number 

radius  vector  in  cross  flow  plane 
(r=1  at  tunnel  walls,  Figs.  28-30), 
also  radial  distance  model-probe 
axis.  Fig.  16 

radius  of  body 
wing  thickness 

Cartesian  coordinates  (Fig.  16) 
distance  along  tunnel  axis 

center  of  pressure  in  wing-fixed 
coordinate  system 

effect  of  wind  tunnel  walls  on  X 

cp 

coordinate  of  line  in  slot  center- 
plane,  where  plenum  pressure  is 
imposed 

angle  of  attack 

disturbance  in  a  due  to  wind  tunnel 
walls 

mean  value  of  in  over  wing  area 
flap  deflection  angle 
leading  edge  flap  deflection  angle 
trailing  edge  flap  deflection  angle 

perturbation  velocity  potential  to 
"exact"  problem,  also  meridional 
angle  model-probe.  Figs.  16-23 

perturbation  velocity  potential  to 
approximate  (filtered)  problem 

wing  sweep  angle  of  quarter  chord 
line 

meridional  angle.  Fig.  28,  also 
flow  deflection  angle,  positive 
out  from  centerline.  Figs.  16-23 
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INTRODUCTION 

During  the  last  decade  a  considerable  part  of  the  aerodynamic  research  and  development 
work  in  Sweden  has  been  directed  to  studies  of  swept  wings  suitable  for  trainer  and  light 
attack  aircraft  configurations.  The  development  and  optimization  of  high-lift  devices  for 
take-off  and  landing  as  well  as  manoeuvre-devices  for  high  subsonic  and  transonic  speeds 
have  played  an  important  role.  The  flow  patterns  encountered  in  these  studies  often  show 
strong  Reynolds'  number  effects  and  it  is  therefore  important  that  each  wind  tunnel  test  is 
carried  out  with  the  largest  model  possible. 

The  model  size  is  limited  by  the  wind  '.lihnel  wall  induced  disturbances  that  can  be 
accepted.  If  the  disturbances  can  be  calculated  properly,  corrections  can  be  applied  for 
the  mean  values  over  the  model  of  the  disturbances  in  angle  of  attack,  velocity  and  in 
some  cases  pitching  moment  (e.g.  for  wing-tail  configurations) .  The  accuracy  of  the  cor¬ 
rected  results  will  be  directly  dependent  on  the  accuracy  with  which  the  disturbance  effects 
can  be  calculated.  The  variation  of  the  disturbances  over  the  model  around  .the  mean  value 
causes  effects  equivalent  to  the  testing  of  a  distorted  model  (with  possible  changes  in 
separation  pattern,  etc.),  which  can  not  be  corrected  for.  If,  however,  these  disturbance 
variations  can  be  calculated  with  reasonable  accuracy,  the  results  can  form  a  basis  for 
the  determination  of  acceptable  model  size.  Relatively  large  uncorrectable  disturbance 
variations  can  in  many  cases  be  tolerated  (for  instance  in  development  and  optimization  of 
high  lift  devices,  studies  of  Reynolds' number  effects,  etc.)  compared  to  what  can  be  accept¬ 
ed  when  taking  stability  and  control  or  performance  data. 

An  exploratory  study  by  Hedman  [2]  with  an  approximate  method  of  swept  wings  in  n  sub¬ 
sonic  closed  test  section  showed  that  the  wall  interference  across  the  span  increased  sig¬ 
nificantly  as  the  sweep  angle  increased.  Hedman  suggested  more  accurate  computations  with 
a  wing  vortex  lattice  method  for  representation  of  the  model  and  the  tunnel  and  that  the 
vertical  position  of  the  model  parts  should  oe  taken  into  account.  Calculations  of  this 
kind  were  made  for  two-dimensional  multi-component  airfoils  by  Eriksson  [ 1 3  ,  for  three- 
dimensional  swept  wings  by  Lind  [3]  and  for  wing-tail  configurations  by  Forslund  [41.  Some 
comparisons  were  made  with  results  calculat  ,d  by  conventional  methods.  Experiment.!  stud¬ 
ies  at  subsonic  speeds  were  made  by  Ridder  [5]  of  wake  blockage  effects  from  a  swept'  wing 
with  and  without  high  lift  devices  and  by  v  Segebaden  [6]  of  the  effects  of  air  leakage 
between  half  model  fuselage  and  reflection  wall. 

As  far  as  calculation  of  wall  interference  in  ventilated  transonic  wind  tunnels  is  con¬ 
cerned..  it  was  considered  that  before  accurate  calculations  could  be  made,  deeper  knowledge 
of  the  boundary  condition  at  the  wind  tunnel  walls  had  to  be  established.  A  theoret.  cal 
and  experimental  study  by  Berndt  and  SbrensSi  [10]  led  to  an  inviscid  theory  of  wall  inter¬ 
ference  in  slotted  test  sections  by  Berndt  [11].  A  method  to  calculate  the  influence  of 
the  wall  boundary  layer  in  slotted  wind  tunnels  has  been  studied  by  LO'fgren  [8,9].  Cross 
flow  characteristics  for  interference-free  ventilated  walls  for  testing  of  lifting  wing- 
bodies  at  low  supersonic  speeds  were  established  theoretically  by  Hedman  and  Rizzi  [14] 
and  experimentally  by  Nyberg  and  SorensSn  [14,15).  A  numerical  method  for  calculation  of 
wind  tunnel  interference  in  slotted  test  sec  .ions  based  on  the  theory  of  Berndt  [11]  has 
been  developed  by  Karlsson  and  Sedin  [12,13].  The  results  hitherto  published  concern  axi- 
symmetric  flows,  but  the  method  can  also  be  used  for  the  angle  of  attack  case,  and  some 
preliminary  results  have  been  obtained  from  the  authors  [16]. 

The  half-model  technique  has  been  developed  and  used  extensively  also  at  transonic 
speeds.  The  only  report  irom  this  work  published  so  far  is  a  theoretical  and  experimental 
investigation  of  the  interference  effects  from  a  support  strut  for  complete  models  on  the 
flow  around  a  reflection-plane  mounted  half  model  made  by  BrSnnstrom  and  Lindau  [7], 

It  is  the  purpose  of  this  paper  to  present  a  brief  review  of  those  of  the  above  mention¬ 
ed  investigations,  which  are  directly  related  to  testing  models  at  high  angle  of  attack. 

For  the  investigations  which  are  more  indirectly  related  to  this  topic,  only  a  brief  ab¬ 
stract  without  figures,  based  mainly  on  the  authorfe  summary,  is  presented. 


SUBSONIC  FLOW 
Two-dimensional  flow 

Calculation  of,  potential  £low  wall  intvi&uience  tfoA.  multi-component  al/i^oili  in  cloied  wind  tunneti  [1] 

The  numerical  method  by  Eriksson  [1],  assuming  incompressible  potential  flow  around  a 
mrlti-component  airfoil  can  be  summarized  in  the  following  points: 

•  Each  airfoil  component  is  approximated  by  a  closed  polygon  whose  corner  points  are 
situated  on  the  original  contour.  The  tunnel  walls  are  similarly  represented  by 
two  straight  "polygons"  (see  Figure  1). 

•  The  straight  polygon  elements  are  given  a  surface  vortex  distribution  that  varies 
linearly  over  each  element  and  is  continuous  at  all  corner  points. 

•  The  surface  vortex  strengths  at  the  corner  points  are  chosen  as  independent  vari¬ 
ables. 

•  A  control  point  is  placed  at  the  midpoint  of  each  polygon  element.  The  velocity 
component  normal  to  the  element  at  this  point  is  taken  as  a  measure  of  the  total 
flow  across  the  element. 
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•  The  boundary  condition  of  zero  normal  velocity  applied  at  every  control  point  results 
in  a  set  of  simultaneous  equations  for  the  vortex  variables.  The  solution  to  these 
equations  gives  the  desired  surface  vortex  distribution,  from  which  the  pressure  dis¬ 
tribution,  the  lift  and  pitching  moment  are  calculated. 

•  The  method  gives  the  inviscid  flow  coefficients  for  both  the  tunnel  case  and  the  free 
flight  case,  and  from  these  results  the  interference  corrections  are  determined. 

The  author  presents  results  for  two  test  cases.  The-  first  is  for  a  symmetric  airfoil 
RAE-102  (t/c  =  0.1  with  a  rather  large  chord  to  tunnel  height  ratio  c/h  =  0.339.  The  re¬ 
sulting  interference  corrections  are  shown  in  Figure  2  in  comparison  with  two  different 
analytic  methods,  a  flat  plate  linear  theory  method  and- a  finite  thickness  method  by 
Goldstein.  The  flat  plate  theory  is  correct  to  order  (c/h)2  but  takes  no  account  of  thick¬ 
ness  or  camber.  Goldstein's  method  is  correct  to  order  (c/h)'  and  includes  the  effect  of 
thickness  and  camber.  None  of  these  methods  is,  however,  in  contrast  to  the  method  by 
Eriksson,  valid  for  large  angles  of  attack,  whan  vertical  displacement  of  the  airfoil  can 
not  be  ignored.  The  results  show  that  Eriksson's  calculations  agree  well  with  Goldstein's 
method  for  small  iingles  of  attack,  but  the  discrepancies  between  the  two  methods  grow 
rapidly  with  increasing  angle  of  attack. 

The  second  test  case  concerns  a  slightly  modified  Wil.liams  airfoil  with  flap  as  depict¬ 
ed  in  Figure  3.  Calculations  were  made  for  c/h  =  0.184  and  0.339.  Some  of  the  results  for 
the  large  chord  are  shown  in  Figure  4 .  The  fact  that  the.  and  corrections  are  non¬ 
linear  in  Cl  ,whlch  is  also  the  case  for  the  small  chord,  is  interesting,  since  all  analyt¬ 
ic  methods  in  practical  use  give  a  linear  dependence.  The  nonlinearity  is  caused  by  the 
vertical  displacement  of  the  airfoil  and  flap  as  the  angle  of  attack  and  flap  deflection 
are  increased.  When  the  flap  approaches  the  lower  wind  tunnel  wall,  it  experiences  a  nega¬ 
tive  lift  increment  that  counteracts  the  normal  tunnel  induced  lift  increase. 


Three-dimensional  flow 

Calzutcution  viatt  intviffi)ivnx.z  on  ijcve  iivzpt  ioing-6  (3] 

The  panel  method  by  Lind  [3]  assuming  incompressible  potential  flow  around  a  thin  wing 
is  summarized  in  the  following  points: 

•  The  real  wing  in  the  wind  tunnel  is  replaced  by  the  mean  surface  of  the  wing  positioned 
in  a  constant  section  thin  walled  tube  modeled  from  the  wind  tunnel  test  section.  The 
combined  model  is  subject  to  a  uniform  wind  blowing  parallel  to  the  test  section  axis. 

•  The  surfaces  are  divided  into  constant  pressure  panels.  The  mutual  influences  from  the 
panels  are  calculated  using  vortex  singularities.  The  vorticity  on  all  panels  is  sim¬ 
ultaneously  determined  from  boundary  conditions  prescribing  attached  flow. 

•  The  calculations  are  fully  three-dimensional. 

•  The  vortices  and  control  points  move  with  the  wing  and  flap  as  the  flap  angle  and  the 
angle  of  attack  are  varied. 

•  In  the  cases  presented  here,  the  vortex  wake  of  the  wing  (and  test  section)  leaves 
the  trailing  edge  in  straight  lines  parallel  to  the  wind  tunnel  axis  and  free  stream 
vector.  More  detailed  models  of  rolled-up  wake  are,  however,  available  in  the  program. 

In  the  present  investigation  a  half-model  installation  in  the  FFA  low  speed  wind  tunnel 
LT  was  studied.  The  circular  test  section  was  therefore  doubled  by  a  reflection  plane,  see 
Figure  5.  Five  wings  were  studied,  see  Figure  6.  The  wings  were  divided  into  24  panels, 

except  for  the  flapped  wing  W45F,  where  48  panels  were  used.  The  circular  walls  of  the 

test  section  were  approximated  by  a  polygon.  The  mathematical  model,  which  was  mostly  used, 
is  depicted  in  Figure  7.  Three  other  models  with  increased  panel  numbers  were  used  to 
check  computational  model  influences. 

The  calculations  were  done  for  the  wing  alone  and  for  the  wing  together  with  the  test 
section.  Some  of  the  results  obtained  with  the  standard  test  section  model  M88  are  shown 
in  Figures  8-11.  The  effect  of  varying  the  test  section  panel  number  was  very  small. 

As  can  be  seen  in  Figure  8  the  span  to  width  ratio  b/B  is  most  important.  The  wall- 

induced  errors  vary  approximately  as  the  square  of  this  variable  within  the  range  invest¬ 
igated.  The  large  effect  on  span  efficiency  factor  ie/e  should  be  taken  into  account  when 
model  size  for  drag  measurements  is  considered. 

The  variation  of  the  wall  induced  angle  of  attack  disturbance  is  presented  spanwise 
along  the  quarter  chord  line  and  chordwise  along  the  root,  midspan  and  tip  chords  in 
Figures  9-11.  The  effect  of  angle  of  attack  is  seen  to  be  small  but  angle  of  sweep  has 
considerable  effect.  The  span  to  tunnel  width  ratio  b/B  has  as  mentioned  above  a  very 
large  effect  for  large  values  of  b/B  .  The  mean  value  over  the  wing  area  Ao/Cl  of  the 
angle  of  attack  disturbance  is  shown  with  full  lines  in  Figures  9,  11a  and  11b.  Corre¬ 
sponding  corrections  calculated  with  conventional  handbook  methods  are  shown  in  the  same 
figures  with  dotted  lines.  The  agreement  is  fairly  good  for  the  plane  wings,  but  not  so 
good  for  the  flapped  wing.  This  indicates  that  for  the  low-lift  case  the  conventional  mean 
value  corrections  are  fairly  good,  but  that  for  high-lift  configurations  more  sophisti¬ 
cated  methods  could  be  necessary  for  accurate  results.  Further  studies  of  the  high  lift 
case  are  suggested. 


Calculation  of  wall  inteAfenen ce  ion  wing-tail  coniigunations  [4] 

The  computer  program  by  Forslund  [ 4 ]  was  developed  for  calculation  and  correction  of 
measured  pitching  moments  for  wing-tail  configurations  in  low  speed  wind  tunnels.  The 
program,  which  is  valid  for  incompressible  potential  flow,  is  summarized  as  follows: 

•  The  wing  and  the  stabilizer  are  replaced  by  their  mean  surfaces  positioned  in  a 
tube  with  thin  walls  modeled  from  the  wind  tunnel  test  section. 

•  The  surfaces  of  the  wing,  which  can  have  a  flap,  the  stabilizer  and  the  wind  tunnel 
are  divided  into  quadrilateral  panels.  The  influence  of  each  panel  on  the  flow  field 
is  simulated  by  a  horse-shoe  vortex.  Each  panel  has  a  control  point,  and  the  Vortic- 
ity  on  all  panels  is  simultaneously  determined  from  the  boundary  condition  prescrib¬ 
ing  tangential  flow  in  all  control  points. 

•  The  vortices  and  control  points  on  the  vjing  and  the  stabilizer  move  with  the  model  as 
the  angle  of  attack  is  varied. 

•  The  trailing  vortex  sheet  is  made  up  from  a  number  of  horse-shoe  vortices.  There  is 
at  present  no  wake  relaxation  in  the  computer  program.  To  position  the  wake  as 
correctly  as  possible,  it  is  recommended  to  use  results  from  visual  observations 

in  the  wind  tunnel. 

The  computer  program,  which  is  written  in  FORTRAN  IV  and  can  deal  with  250  panels  dis¬ 
tributed  on  the  wing, the  stabilizer  and  the  wind  tunnel  surface, is  used  to  compute  the  dif¬ 
ference  in  angle  of  downwash  at  the  position  of  the  stabilizer  with  and  without  the  wind 
tunnel  walls.  It  is  assumed  that  a  change  of  downwash  at  the  position  of  the  stabilizer 
is  ejuivalent  to  a  change  in  the  stabilizer  angle  in  so  far  as  the  same  moment  contribu¬ 
tion  results.  Corrected  C_,-a  curves  are  obtained  by  linear  interpolation  between  Cn~a 
curves  for  different  stabilizer  angles. 

An  example  of  where  this  correction  method  has  been  applied  is  shown  in  Figure  '■■2.  A 
model  of  a  swept  wing-tail  combination  in  a  high-lift  configuration  was  tested  as  a  full 
modei  at  F+W  in  Emmen  and  as  a  half  model  '  che  same  model)  at  FFA.  The  agreement  between 
the  v+w  results  and  the  corrected  FFA  res  ilts  is  very  good.  A  corresponding  correction 
using  conventional  handbook  methods  gives  3t%  of  the  panel  method  correction  at  a  =  0°  and 
67%  at  a  »  16°. 


Expen -'mental  itudy  of,  utake.  blockage,  effecti  fnom  a  iuicpt  wing  with  and  Without  high  lift  devices  [5] 

Ai  a  background  for  his  experimental  study  Ridder  [5]  presents  a  review  of  the  current 
methods  of  treating  the  wake  blockage  problem  and  discusses  possible'  wake  flow  patterns  for 
lifting  wings.  It  is  concluded  that  the  real  wake  flow  behind  three-dimensional  wings  can 
in  most  cases  not  be  characterized  with  the  current  oversimplified  wake  models. 

The  wind  tunnel  model  used  for  the  measurements  was  a  35°  swept  (25%  chord)  ,  aspect  • 
ratio  4  wing  with  10%  thick  symmetrical  section  and  tested  in  three  configurations:  Plane, 
with  leading  edge  flap  only  and  with  leading  edge  slat  and  double  slotted  trailing  edge 
flap.  The  model  was  tested  as  a  half  model  on  a  reflection  plate  and  the  test  section 
wall  pressures  were  measured  along  a  streamwise  line  on  the  tunnel  ceiling  close  to  the 
chord  plane  of  the  wing  at  zero  angle  of  attack,  see  Figure  13. 

The  measured  drag  polars  for  the  three  wings  are  shown  in  Figure  14  and  the  wake  block¬ 
age  induced  pressure  distributions  along  the  test  section  wall  for  the  configurations  with 
high  lift  devices  in  Figure  15.  The  wake  blockage  induced  pressure  coefficients  ACp  were 
obtained  from  the  measured  pressures  by  subtraction  of  the  pressures  measured  with  the 
clean  configuration  set  at  zero  angle  of  attack  and  by  subtraction  of  the  lift  induced 
pressures  calculated  from  the  results  obtained  by  Lind  13], 

* 

The  pressure  distributions  in  Figure  15  support  the  concept  of  an  upstream  and  a  down¬ 
stream  velocity  level  with  a  rather  pronounced  acceleration  in  the  vicinity  of  the  model. 
The  effective  location  of  a  source  to  represent  the  wake  flow  seems  to  be  quite  a  distance 
downstream  of  the  model  reference  point.  Accordingly,  the  assumption  of  a  mean  value  of 
the  two  velocity  levels  should  not  be  valid  when  estimating  the  effective  free  stream 
velocity  at  the  model  position.  It  is  shown  in  the  report,  however,  that  the  velocity 
gradient  at  the  test  section  wall  is  considerably  higher  than  on  the  centerline.  The  con¬ 
clusion  should  be  that  the  mean  value  concept  will  cause  an  overestimation  of  the  wake 
blockage  effect. 

The  difference  in  upstream  and  downstream  static  pressure  coefficient  levels  is  called 
ACpjnaj;-  This  pressure  coefficient  is  correlated  with  the  increase  in  flow  separation  drag 
ACp.  The  ratios  ACpmax/ACp  are  listed  in  Figure  15.  This  ratio  is  found  to  vary  from 
0.12  to  0.24.  The  corresponding  value  calculated  by  the  classical  Maskell/Vayssaire  meth¬ 
ods.  for  this  test  configuration  is  ACp^/ACp = 0.585.  It  seems  likely  that  these  methods 
could  provide  the  upper  limit  for  wake  blockage  effects  particularly  at  post  stall  condi¬ 
tions.  However,  for  wings  of  the  type  investigated,  the  above  methods  would  seem  to  over¬ 
estimate  the  actual  wake  blockage  effect  by  at  least  a  factor  3  up  to  and  above  stall. 

Ridder  proposes  direct  measurement  of  the  wake  blockage  induced  pressure  difference 
upstream  and  downstream  of  the  model  as  an  alternative  to  the  methods  of  Maskell  and 
Vayssaire.  The  method  has  been  used  by  Ridder  for  several  years  in  wind  tunnel  tests 
with  thrust  reversers  where  blockage  effects  are  highly  pronounced. 
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Ex.pt/Ujnejntat  study  o(,  the  ejects  o J  cua  llow  leakage  between  half,  model  iuselagz  and  /levitation  wall 

When  testing  half  models  it  is  very  important  to  bring  friction  between  the  model  and 
its  reflection  wall  to  a  minimum.  At  the  same  time  it  in'  also  important  to  keep  the  leak¬ 
age  of  air  between  model  and  wall  as  small  as  possible  so  that  the  flow  field  remains  un¬ 
distorted. 

In  connection  with  a  high  lift  investigation  of  a  25°  swept  wing  half  model  v.  Sege- 
baden  [6]  has  studied  the  effects  of  air  flow  leakage  between  the  fuselage  and  the  reflec¬ 
tion  wall. 

The  effects  of  four  different  sealing  conditions  or  gaps  between  model  and  wall  have 
been  investigated.  The  different  sealing  conditions  show  marked  differences  that  point  to 
the  importance  of  using  carefully  designed  sealings  between  the  half  model  and  the  wall  if 
the  results  are  going  to  be  used  as  absolute  values  and  comparable  with  results  from  other 
tunnels  and  models.  It  is  also  important  to  design  the  sealing  inside  the  fuselage  so  that 
pressure  loads  on  the  sealing  do  not  de-load  the  true  fuselage. 

It  is,  however,  so  difficult  to  estimate  the  effects  of  leakage  that  it  is  probably 
better  to  accept  a  certain  amount  of  friction  from  an  efficient  seal. 

The  results  do  also  show  that  the  effects  on  lift,  drag  and  pitching  moment  characterist¬ 
ics  do  not  necessarily  appear  at  the  same  time  due  to  the  fact  that  large  forces  may  not 
give  any  moments  at  all  if  they  balance  around  the  reference  axis. 


TRANSONIC  FLOW 

Wind  tunnel  wall  conditions  ; 

flow  paopeAties  o{  slotted  mile  [10,11] 

A  theoretical  and  experimental  study  of  the  flow  properties  of  slotted  walls  for  tran¬ 
sonic  test  sections  made  by  Berndt  and  SSrensdn  [10]  led  to  a  new  inviscid  theory  of  wall 
interference  it.  .lotted  test  sections  by  Berndt  [11].  The  classical  theory  of  longitudinal¬ 
ly  slotted  walls,  which  substitutes  an  approximate  homogeneous  wall  boundary  condition  for 
the  true  mixed  conditions,  is  extended  in  several  respects.  Based  on  the  experimental 
findings  at  FFA  ar  inviscid  flow  model  is  adopted,  in  which  the  flow  entering  the  slot 
might  split  up  into  two  different  streams,  one  leaving  the  slot  penetrating  into  the  plenum 
chamber  as  a  thin  jet  and  one  re-entering  into  the  test  section,  admitting  quiescent  air 
from  the  plenum  chamber  into  the  test  section  and .inducing  a  longitudinal  separation  bubble 
at  plenum  pressure  along  the  slot  and  adjacent  parts  of  the  test  section  wall.  The  three- 
dimensional  analysis,  based  on  the  assumption  that  the  slots  are  narrow,  retains  quadratic 
cross-flow  terms  in  the  pressure  equation  and  allows  the  slots  to  be  few  in  number  and 
have  non-uniform  distribution  and  geometry.  A  family  of  homogeneous  boundary  conditions 
is  obtained,  each  of  successively  higher  accuracy.  Application  to  the  design  of  interfer¬ 
ence-free  transonic  test  sections  is  discussed. 

Unsteady  effects  also  are  considered. 


Slotted  wall  boundasiy  layea  edicts  [8,9] 

Starting  from  an  idealized  physical  model  of  subsonic  flow  past  a  slender  axisymmetric 
body  in  a  longitudinally  slotted  wind  tunnel  of  circular  cross  section  and  negligible  wall 
thickness,  Lofgren  [8,9]  has  derived  a  homogeneous  wall  boundary  condition.  This  is  to 
be  used  at  the  wall  in  connection  with  the  usual  mathematical  model  of  non-viscous  flow  in 
order  to  produce  approximately  the  same  flow  field  in  the  central  part  of  the  tunnel  as 
the  field  obtained  from  the  direct  mathematical  representation  of  the  idealized  physical 
model.  The  boundary  layer  flow  is  assumed  to  be  laminar,  non-viscous  and  in  Reference  [8] 
to  have  an  one-step  velocity  distribution,  while  in  the  method  of  Reference  [9]  there  is  a 
free  choice  of  the  velocity  profile. 


JnteA(esience-(Aee  wall  boundaAy  condition  {oa  tests  oi  lifting  wing-body  models,  at  low  supensonic  speed 

[14,15] 

The  interference-free  flow  field  at  low  supersonic  speed  around  a  lifting  delta-wing- 
body  configuration  (see  Figure  17)  was  studied  theoretically  by  Hedman  and  Rizzi  [14]  and 
experimentally  by  Nyberg  and  SorensSn  [14,15],  The  undisturbed  pressures  and  flow  angles 
were  obtained  at  radial  locations  in  relation  to  the  model,  where  in  wind  tunnel  tests  the 
walls  are  normally  situated.  The  flow  properties  were  used  to  establish  criteria  for  cross 
flow  characteristics  of  interference-free  wave-attenuating  ventilated  wind  tunnel  walls. 

Theoretical  calculations  were  made  for  the  wing-body  configuration  at. small  angles  of 
attack  with  a  relaxation  method  based  on  the  transonic  small  perturbation  (TSP)  equation 
by  Hedman  and  in  one  axisymmetric  sample  case  for  'the  body  alone  with  a  method  based  on 
the  solution  of  the  full  gas  dynamics  equations  (the  Euler  equations]  in  a  time -dependent 
finite  volume  formulation  by  Rizzi.  In  the  theoretical  results  obtained  using  the  TSP 
method  the  shocks  are -smeared  out,  the  pressure  and  flow  angle  peak  values  are  strongly 
underestimated,  and  the  agreement  with  experimental  data  is  generally  only  qualitative. 

In  the  results  obtained  by  the  finite  volume  method  for  the  body  alone  at  zero  angle  of 
attack  the  pressure  rise  through  the  shock  agrees  well  with  experiments  but  the  shock  posi- 
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tion  is  not  predicted  correctly.  Details  about  the  theoretical  results  will  not  be  present¬ 
ed  in  this  review  as  the  calculations  only  cover  zero  and  small  angles  of  attack.  The  pro¬ 
blems  experienced  in  these  calculations  indicate,  however,  that  the  development  of  accurate 
theoretical  correction  methods  for  high  angle  of  attack  testing  with  this  kind  of  flow 
pattern  (supersonic  flow  between  the  model  and  the  wind  tunnel  wall  with  rather  strong 
oblique  shocks)  will  be  very  difficult. 

The  experimental  test  set-up  is  shown  in  Figure  16.  The  model  can  be  translated  in  the 
x-direction  along  the  tunnel  centerline  and  it  can  also  be  rolled  an  angle  if  around  the 
x-ayi‘s.  The  radial  position  of  the  probe  used  to  determine  pressure  coefficient  Cp  and 
flow  inclination  angle  8  is  variable.  The  model  (Figure  17)  has  a  4%  thick  delta  wing. 

Examples  of  the  experimentally  determined  flow  field  parameters  Cp  and  6  are  depicted 
in  Figure  18  together  with  Schlieren  pictures  in  the  same  linear  scale  for  a  high  angle  of 
attack  case  at  M  =  1.20.  Data  points  obtained  downstream  of  the  reflected  bow  wa\e  are 
disturbed  by  wall  interference  and  are  marked  with  filled  symbols.  This  kind  of  interfer¬ 
ence  is  greatly  reduced  at  B  =1.30  and  therefore  results  for  this  Mach  number  were  used 
in  Figures  19-22  to  demonstrate  the  effect  of  angle  of  attack.  It  can  be  seen  in  Figures 
19  and  21  that  the  influence  of  angle  of  attack  is  small  on  the  strength  of  the  body  bow 
shock  but  noticeable  on  its  position  above  and  below  the  model.  The  a-influence  on  the 
wing-induced  flow  field  is  very  large. 

Tne  highest  pressures  immediately  downstream  of  the  body  bow  shock  and  of  the  wing 
leading  edge  shock  denoted  Cpmax  and  the  lowest  pressure  immediately  upstream  of  the 
wing  trailing  edge  shock  denoted  Cpm^n  (see  Figure  19)  are  plotted  versus  angle  of  attack 
in  Figure  20.  The  corresponding  flow  angles  are  depicted  in  Figure  22  in  the  same  way.  The 
a-infiuence  is  clearly  seen.  The  hatched  curves  in  the  wing  diagrams  representing  the  mean 
values  of  Cp  and  9  above  (cp  =  0)  and  below  ( g>  =  180°)  the  model  demonstrate  the  lift  in¬ 
duced  "volume"-ef feet  which  for  high  angles  of  attack  is  much  larger  than  the  geometrical 
volume  effect  at  a  =  0.  It  is  also  noted  that  this  pressure  disturbance  is  nearly  linear 
with  a  and  very  close  to  the  pressure  at  the  side  of  the  model  ( cp  =  9 0°)  .  The  large  in¬ 
crease  in  Cp  and  0  in  the  body  diagrams  at  a  ”28:3  below  the  model  (<P=  180°)  is  insignif¬ 
icant  and  indicates  only  that  the  wing  leading  edge  shock  has  moved  upstream  and  joined  the 
body  bow  shock,  which  also  can  be  seen  in  the  Schlieren  picture  of  Figure  18. 

The  characteristics  of  a  ventilated  test  section  wall  are  frequently  described  b>  the 
relationship  between  the  pressure  drop  of  the  flow  when  it  passes  through  the  wall  a.id  the 
flow  angle  of  this  flow  at  the  wall.  Each  diagram  of  Cp  versus  8  in  Figure  23  can  there¬ 
fore  be  interpreted  as  the  required  character istics  of  an  interference-free  wall  for  the 
location  and  flow  condition  investigated.  The  influence  of  angle  of  attack  and  of  distance 
from  the  model  at  =  1.20  are  shown.  It  can  be  seen  that  the  disturbance  due  to  lift  de¬ 
cays  faster  than  the  disturbance  due  to  geometrical  volume.  Other  results,  not  reproduced 
here,  show  that  the  slope  of  the  curve  varies  with  Mach  number.  The’ character  of  the  curves 
are  generally  linear  or  slightly  curved  with  the  largest  slope  in  the  outflow  region.  This 
is  in  contrast  to  the  cone-cylinder  flow  field  hitherto  used  as  requirement  for  tunnel  wall 
characteristics,  where  the  largest  slope  is  in  the  inflow  region. 

In  some  areas  of  the  flow  field  outflow. is  combined  with  underpressure  and  inflow  with 
overpressure.  Perforated  walls  with  constant  porosity  and  plenum  pressure  can  not  match 
this  requirement.  It  is  suggested  by  the  authors  that  combined  slotted-perforated  should 
be  investigated.  Another  possibility  is  of  course  to  use  perforated  adaptive  walls  with 
locally  variable  porosity  and/or  pressure.  The  use  of  solid  flexible  walls  will  be  diffi¬ 
cult  as  the  abrupt  local  slope  changes  at  the  shocks  (see  Figure  21)  have  to  be  at  least  5° 
for  high  angle  of  attack  testing. 


Wall  interference 

Calculation  oij  uxM  inteA&ejicncc.  in  itotted  tut  actio  ni  -  AxiiijmmetAic  $louu  (12,13] 

f 

Using  the  above  mentioned  inviscid  theory  by  Berndt  [11],  Karlsson  &  Sedin  [12,13]  have 
developed  a  method  for  calculation  of  transonic  wind  tunnel  interference  in  slotted  test 
sections.  The  numerical  program  can  also  be  used  in  an  inverse  mode  for  design  and  analysis 
of  optimal  slot  shapes.  The  results  published  so  far  concern  axisymmetric  flows,  but  the 
program  can  also  be  used  for  the  angle  of  attack  case,  which  will  be  commented  upon  in  the 
next  section. 

The  theory  by  Berndt  [11]  is  built  on  the  calculation  of  a  filtered  perturbation  veloc¬ 
ity  potential,  neglecting  higher  order  variation  caused  by  the  slots  and  the  walls.  The 
theory  results  in  a  homogeneous  wall  boundary  condition  including  the  dependence  on  slot 
geometry.  The  small  perturbation  potential  equation  is  iteratively  integrated,  repeatedly, 
with  the  rapid  finite-difference  method  of  Karlsson  &  Sedin  [13],  using  this  wall  bound¬ 
ary  condition. 

Some  calculated  results  from  analysis  and  design  of  slot  shapes  in  flows  with  axisym¬ 
metric  bodies  at  subsonic  free  stream  Mach  numbers  published  in  Reference  [13]  will  be 
presented  here  and  also  some  unpublished  data  comparing  calculated  and  experimental 
results  [16]. 

Three  different  bodies,  generated  by  a  generalized  parabolic  arc,  were  considered  in 
the  calculations.  Only  results  for  Body  I  will  be  presented  in  this  review.  This  body 
had  the  maximum  diameter  at  50%  of  the  length,  a  slenderness  ratio  of  1/(6/2j  and  was 
truncated  at  5/6  x  length  with  the  cylindrical  sting  following  from  that  station,  see 
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Figure  24.  The  body  was  studied  in  one  small  and  one  large  version  with  blockage  area 
ratios  of  0.42  and  2.23%  respectively.  The  tunnel  was  assumed  to  be  circular  in  cross 
section,  with  the  wall  defined  by  r  =  1  and  provided  with  8  uniformly  distributed  slots. 

As  a  start  a  slot  width  distribution  of  a  slot  called  the  "FFA"  slot  as  depicted  in  Figure 
25  was  studied.  This  slot  shape  is  close  to  the  slot  of  the  existing  transonic  wind  tunnel 
HT  at  the  Aeronautical  Research  Institute  of  Sweden^  (FFA) . 

The  small  model  of  Body  I  with  0.4%  blockage  had,  vhen  calculated  in  a  test  section 
with  the  non-optimal  "FFA"-slot,  nearly  negligible  interference,  while  the  large  model 
with  2.23%  blockage  had  dramatically  increased  interference  as  shown  in  Figure  24.  The 
level  of  interference  had  to  be  diminished  by  at  least  an  order  of  magnitude  to  be  accept¬ 
able.  An  optimal  slot  shape  was  calculated  by  running  the  inverse  mode  design  process. 

This  slot  was  named  M98(I),  which  means  that  the  optimal  design  reference  Mach  number  is 
0.98  and  that  the  slot  is  designed  for  Body  I.  The  shape  of  the  optimal  slot  M98(I)  is 
found  in  Figure  ’.5.  For  this  slot,  as  shown  in  Figure  26,  the  agreement  between  the  cal¬ 
culated  free  air  reference  pressure  distribution  and  that  of  the  tunnel  run  is  excellent 
on  the  body  surface  and  also  in  good  agreement  on  the  wall. 

It  was  observed  that  for  the  large  model  of  Figures  24  and  26  a  slight  discontinuity  in 
the  pressure  gradients  could  be  recognized  at  the  inflow  section,  which  is  a  symptom  of 
shortcomings  of  the  assumed  inflow  boundary  condition.  The  large  Body  I  is  currently 
being  experimentally  tested  at  FFA  by  H.  SdrensSn  to  get  results  for  further  work  on  the 
subject.  Some  preliminary  results  comparing  calculated  and  experimental  wall  pressure 
distributions  [16]  are  shown  in  Figure  27.  The  agreement  in  shape  of  the  curves  is  good, 
but  there  is  a  shift  in  level,  which  could  be  a  result  3f  the  simplified  inflow  conditions 
used  in  the  calc  ilations.  Another  reason  for  the  disagreement  could  be  that  no  viscous 
effects  have  bee i  considered  in  the  slot  width  distribution.  The  experiments  will  continue 
with  surface  pressure  measurements  on  the  model,  and  with  detailed  flow  surveys  in  the 
slots.  Later  studies  of  the  effect  of  slot  depth  are  planned  and  also  tests  with  optimal 
slots. 

Calculation  o&  mJl.1  inteai^eJicncc  in  ilotted  tut  le.ctioni  -  Delta  ting  at  angle,  o ^  attack  [16] 

The  axisjm.ietric  results  by  Karlsson  &  Sedin  [13]  presented  in  the  preceding  section 
look  very  promising  and  the  theoretical  and  experimental  work  will  continue  to  clarify  re¬ 
maining  problems.  However,  a  more  interesting  and  natural  continuation  will  be  to  go  on 
with  three-dimensional  flow  fields,  especially  when  lif :ing  configurations  are  considered. 
Berndt's  theory  [11]  and  the  numerical  method  of  Karlsson  &  Sedin  [13]  are  capable  of 
handling  three-dimensional  flow  fields.  Exploratory  calculations  with  lifting  flow  fields 
have  already  been  performed  and  some  of  the  results  have  kindly  been  made  available  by 
Karlsson  &  Sedin  [16].  _ 

Calculations  have  been  made  for  a  delta-wing  as  shown  in  Figure  28.  The  cross-flow 
condition  at  the  inner  boundary,  which  is  a  cylinder  with  r  =  0.5,  is  for  the  free-air  case 
obtained  by  slender-body  theory.  The  test  section  is  circular  and  has  8  similar  slots  of 
constant  width  and  with  9.2%  open  area.  Calculated  pressure  distributions  at  the  inner 
boundary  and  at  the  tunnel  wall  for  a  =  8.75°  and  M  =  0.95  are  depicted  in  Figures  29.  and 
30.  The  pressures  at  the  inner  boundary  are  higher  for  the  tunnel  case  than  for  the  free- 
air  case.  This  could  indicate  that  the  open  area  of  the  slots  is  too  large.  Also  from 
the  wall  pressures  appreciable  wall  interference  can  be  observed.  These  first  results  are 
certainly  very  encouraging  and  it  is  to  be  hoped  that  this  important  work  can  be  continued 
as  soon,  as  possible  and  that  suitable  supporting  experiments  will  be  carried  out. 


sort  interference 


InteAfieAencc  on  a  kali  model  ( Korn  a  hull  model  tuppoAt  itAut  [7] 

A  theoretical  and  experimental  investigation  of  the  interference  in  a  wind  tunnel  from 
the  support  strut  for  complete  models  on  the  flow  around  a  reflection  plane  mounted  half 
model  has  been  made  by  Brannstrom  &  Lindau  [7], 

In  the  theoretical  part  a  computer  program,  which  could  handle  up  to  three  components 
(wing,  tunnel  and  support  strut),  was  used  to  calculate  the  interference  effects.  The 
components  are  represented  by  a  mean  surface  distribution  of  quadrilateral  vortex-ring 
elements,  the  strength  of  which  is  determined  from  the  boundary  condition  of  zero  normal 
component  of  velocity  at  a  control  point  in  each  element.  The  trailing  vortices  are  re¬ 
laxed  into  the  local  flow  direction  by  an  iterative  procedure. 

The  calculations  were  carried  out  for  an  angle  of  attack  of  5  degrees  with  four  dif¬ 
ferent  configurations:  without  tunnel,  tunnel  without  support  strut,  tunnel  with  vertical 
and  tunnel  with  horizontal  support  strut.  Calculations  were  made  at  Mach  number  zero  and 
by  Prandtl-Glauert  transformation  at  Mach  numbers  0.6  and  0.8.  The  results  show  notice¬ 
able  interference  from  the  tunnel  walls,  but  very  small  interference  from  the  support  strut. 
The  calculated  results  for  the  strut  cases  are  plotted  together  with  the  experimental  re¬ 
sults  for  comparison  in  Figures  32,  34  and  35.  . 

The  experimental  investigation  was  carried  out  with  a  half  model  of  a  wing-body  combin¬ 
ation.  The  wing  had  an  aspect  ratio  of  4,  a  taper  ratio  of  0.4  and  a  sweep  angle  at'  the 
quarter-chord  line  of  35°.  Forces  and  pressure  distributions  at  three  wing  sections  were 
measured.  Four  model  test  cart  configurations  were  tested:  without  strut,  with  normal 
strut  in  vertical  or  horizontal  position  and  a  strut  simulating  the  strut  in  the  proposed 
FFA  T-2  tunnel.  The  test  arrangement  is  depicted  in  Figure  31. 


6-8 


Some  results  from  the  force  measurements  are  depicted  in  Figures  32-35.  The  strut- 
effecus  on  C^,  C„,  and  Cr  are  small  at  low  angles  of  attack  but  are  quite  large  at  high 
angle  of  attack.  For  Cp  the  interference  effects  are  noticeable  over  the  whole  range  of 
angle  of  attack  and  Mach  number.  This  result  in  combination  with  the  theoretical  result, 
that  there  is  practically  no  downwash  or  vort.ex  relocation  interference,  indicates  that 
it  is  primarily  a  blockage  effect  of  the  'support  strut.  This  is  supported  by  the  fact, 
that  the  T-2  strut,  which  has  the  largest  volume  with  part  of  it  close  to  the  model,  pro¬ 
duces  the  largest  interference  effects  and  a.l'30  that  the  normal  strut  interference  is  in¬ 
dependent  of  whether  the  strut  is  vertical  or  horizontal. 

This  result  is  interesting  and  indicates  that  very  likely  there  could  be  interfetence 
effects  of  this  kind  in  many  tunnels  also  in  full  model  tests. 

CONCLUDING  REMARKS 

It  has  been  demonstrated  that  for  subsonic  incompressible  flow  the  use  of  panel  methods 
for  the  calculation  of  wind  tunnel  wall  interference  will  in  many  cases  provide  more  accu¬ 
rate  corrections  than  the  conventional  methods.  This  is  certainly  the  case  for  instance 
for  wings  with  high-lift  devices  and  wing-tail  configurations,  where  various  model  parts 
are  vertically  displaced  from  the  horizontal  symmetry  plane  of  the  wind  tunnel.  Calcula¬ 
tions  of  this  kind  can  also  display  the  variation  of  the  tunnel-induced  disturbances  over 
the  model- surface ,  the  knowledge  of  which  is  of  importance  for  the  decision  on  proper  model 
size  with  respect  to  the  purpose  of  the  test. 

Tt e  classical  methods  by  Maskell/Vayssair 3  seem  to  provide  the  upper  limit  for  wake 
blockage  effects,  particularly  at  post  stall  conditions.  For  an  investigated  conventional 
swept -wing  with  high- lift  devices  these  methods  overestimate  the  wake  blockage  effect  by 
at  least  a  factor  of  3  up  to  and  above  stall.  As  an  alternative  to  these  methods,  it  is 
proposed  to  measure  the  wake  blockage  induced  pressure  difference  upstream  and  downstream 
of  the  model. 

An  improved  theory  for  the  wall  boundary  condition  at  transonic  slotted  test  section 
wall j  has  been  presented. 

Criteria  for  cross-flow  characteristics  of  interference-free  wave-attenuating  ventil¬ 
ated  wind  tunnel  walls  have  been  established.  The  results  might  also  have  application  in 
the  design  of  adaptive  walls  for  cases  where  strong  shocks  extend  to  the  walls. 

Eased  on  the  new  wall  boundary  condition  a  numerical  method  for  calculation  of  ransonic 
wind  tunnel  interference  in  slotted  test  sections  has  been  developed.  The  program  has  been 
used  mainly  for  axisymmetric  calculations,  but  exploratory  calculations  with  lifting  flow 
fields  have  been  performed  with  encouraging  results. 

If  during  a  half  model  test  in  a  wind  tunnel  the  full  model  support  strut  remains  in 
its  normal  position,  noticeable  strut-interference  effects  on  the  half  model  flow  could  be 
encountered. 
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Figure  5.  Test  section  of  low  speed 
wind  tunnel  LT  at  FFA  [ 3 ] . 
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Figure  6.  Wing  plan  forms  investigated  [3] 
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Figure  7.  Wing  W45  at  a  =  25  in 
test  section  M88  [3]. 
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Figure  8.  Effect  of  wind  tunnel  walls  on 
lift,  drag  and  center  of  pres¬ 
sure,  wing  W45  [3], 
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Figure  12.  Pitching  moment  curve, 
£AAB-wing  7803  with 
stabilizer  [4]. 
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Figure  15.  Wake  blockage  Induced  pres¬ 
sure  distribution  along  the 
test  section  wall  (ceiling)  [5] 
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Figure  17.  Wing-body  model  [15] 
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Figure  18.  Flow  field  and  Schlieren  picture  M 
a  =  27?9,  r  =  103.2  [ 15] . 
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Figure  24.  7 ressure  distributions  of 

body  I  (large)  in  tunnel  with 
"FFA"  slot  at  M  =  0.98  [13]. 
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Figure  26.  Pressure  distributions  of 
body  I  (large)  in  tunnel 
with  optimal  slot  M98  (I) 
at  =  0.98  [13]. 
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Figuri-  25.  Comparisons  between  the  optimal 
slot  M98  (I)  and  the  "FFA"  slot 
1 13]. 


Figure  27.  Comparison  of  calculated 
and  experimental  wall 
pressure  distributions  [16]. 
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Figure  28,  Computational  field  model  for 
the  angle  of  attack  case  [16], 


PRESSURE  (CpI  AT  INNER  BOUNDARY, 
r  =  0.5 


/  \  0  =  90°  (UPPER) 

\V__ _ _ 


I'  0 =-90°  (LOWER) 


M  =  0. 95 

- 

0=6.75° 


- TUNNEL  CALCULATIONS 

-  FREE-AIR  — 


Figure  29.  Calculated  pressure  dist.cibu 
tions  at  inner  boundary  for 
an  angle  of  attack  case  116] 
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Figure  30.  Calculated  pressure  distribu¬ 
tions  at  tunnel  wall  for  an 
angle  of  attack  case  [16], 


Figure  31.  Test  arrangement  [7]. 
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Figure  32.  Strut  interference  effects 
on  CL  [7], 


Figure  33.  Strut  interference  effects 
on  CD  [7]. 
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SUMMARY 

Only  relevant  work  completed  or  in  progress  since  the  FOP  meeting  on  Wind  Tunnel  Design  and  Testing 
Techniques  (lg75)  is  covered.  It  includes  the  use  of  adaptive  walls  (Southampton  University),  a  panel 
method  of  model  and  wake  representation  for  a  2-D  model  in  a  wind  tunnel  with  solid  walls  (British 
Aerospace),  the  use  of  measured  pressure  distributions  on  tunnel  floor  and  roof  also  for  a  2-D  model  and 
solid  walls  (RAE),  a  vortex-lattice  representation  of  the  tunnel  walls  to  take  account  of  wake  curvature 
(British  Aerospace),  interference  limitations  on  tests  on  V/STOL  models  with  lifting  jets  (RAE),  and  work 
on  blockage  corrections  on  models  with  reverse  thrust  (British  Aerospace). 

Some  discussion  is  offered  on  the  limitations  on  the  validity  of  current  methods  for  determining 
wind  tunnel  corrections  and  it  is  argued  that  these  limitations  are  least  severe  with  the  use  of  adaptive 

walls. 

1 .  INTRODUCTION 

Since  the  AGARO  FOP  meeting  on  Wind  Tunnel  Oesign  and  Testing  Techniques  of  1975  the  work  in  the 
U.K.  relevant  to  the  subject  of  wind  tunnel  corrections  for  high  angles  of  attack  has  been  somewhat 
limited. 

This  work  can  be  sumaarised  as  follows: 

a)  Work  at  Southampton  University  by  Dr.  Goodyer  on  the  use  of  adaptive  walls. 

b)  Work  at  Biitish  Aerospace,  Kingston-Brough  Oivision,  by  Holt,  Storey  and  Cross  on  the 
corrections  for  a  two-dimensional  model  in  a  wind  tunnel  with  solid  walls.  This  involves 
a  panel  representation  of  the  model  and  wake  coupled  with  a  transformation  which  dispenses 
with  the  need  to  consider  more  than  one  image. 

c)  Work  by  Ashill  and  Weeks,  R.A.E.,  using  measured  pressure  distributions  on  the  tunnel  roof 
and  flooi  to  determine  the  corrections  for  a  two-dimensional  model  in  a  tunnel  with  solid 
walls  up  to  high  subsonic  Mach  numbers. 

d)  Work  at'British  Aerospace,  Weybridge,  by  Storey  involving  a  vortex-lattice  distribution  to 
represent  the  tunnel  walls  but  with  the  limited  aim  of  determining  the  interference  effects 
arising  from  differences  of  wake  curvature  between  the  results  in  ‘free  air1  and  tunnel. 

e)  Work  by  Trebble,  R.A.E.,  using  wind  tunnel  measurements  for  models  with  lifting  jets  to 
determine  the  conditions  limiting  the  validity  of  wind  tunnel  testing  of  V/STOL  configurations. 

f)  Work  at  British  Aerospace,  Hatfield,  by  Farley,  McRae  and  Bryant  on  blockage  corrections  on 
models  with  reverse  thrust. 

The  salient  features  of*  each  of  the  above  will  be  discussed  in  turn  in  what  follows. 

2.  work  of  dr.  goodyer1,-2 ,3,1  Southampton  university. 

THE  SELF-STREAMLINING  WIND  TUNNEL. 

The  essential  principle  of  Goodyer's  tunnel  is  that  the  pressure  distributions  and  shapes  of  the  top 
and  bottom  flexible  walls  are  measured  with  a  model  in  the  tunnel  and  the  pressure  distributions  are 
compared  with  those  calculated  for  an  unlimited  flow  external  to  the  tunnel  walls.  If  the  two  pressure 
distributions  do  not  agree  then  the  wall  shapes  are  changed  so  as  to  bring  the  distributions  into  closer 
agreement.  The  measurements  are  repeated  until  an  acceptable  level  of  agreement  between  the  pressure 
distributions  is  attained.  Perfect  agreement  implies  that  the  flow  in  the  tunnel  is  then  the  same  as 
for  unlimited  flow  and  the  tunnel  walls  are  streamlines  for  unlimited  flow.  The  method  includes  an 
allowance  for  the  effective  displacement  of  the  tunnel  walls  due  to  the  boundary  layers  on  them.  The 
process  of  data  acquisition,  computation  and  operation  of  the  jacks  determini, ig  the  wall  shapes  can  be 
made  automatic  by  means  of  a  suitably  programmed  on-line  computer.  The  number  of  iterations  required 
has  been  reduced  to  two  or  three  by  the  application  of  a  process  devised  by  Judd. 

Up  to  the  present  Goodyer's  work  has  been  essentially  on  two-dimensional  models  in  small  tunnels 
6  in.  between  top  and  bottom  walls.  An  early  but  striking  result  that  he  obtained  at  low  speeds  is 
illustrated  in  Fig.  1  which  shows  results  he  obtained  for  the  flow  past  a  circular  cylinder.  It  will 
be  seen  that  the  self-streamlining  tunnel  gave  satisfactory  results  even  with  a  cylinder  of  29.3% 
blockage.  Some  results  obtained  on  an  aerofoil  with  tunnel  height  (h)/ chord  (c)  =1.1  at  incidences 
up  to  and  beyond  the  stall  are  illustrated  in  Fig.  2.  Here  we  see  encouraging  agreement  with  measurements 
made  on  the  same  model  in  a  much  larger  Langley  tunnel  (L.T.P.T.)  with  h/c  =  16,  except  near  and  beyond 
the  stall . 
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In  his  low  speed  tunnel  some  18  jacks  are  used  on  the  flexible  walls  extending  three  or  four  chord 
lengths  up  and  downstream.  The  method  allows  for  truncation  errors  due  to  the  finite  length  of  the 
flexible  walls  and  can  also  incorporate  corrections  for  any  small  differences  that  may  remain  between 
the  measured  and  desired  wall  pressure  distributions  at  a  suitable  stage  in  the  Iterative  process. 

There  appears  to  be  no  problem  in  extending  the  method  to  transonic  speeds^,  apart  from  the 
limitation  that  it  cannot  be  pursued  much  beyond  the  point  where  a  strong  shock  wave  would  extend 
outwards  beyond  the  wall  position  In  unlimited  flow. 

Goooyer  plans  to  extend  his  work  to  three  dimensional  models,  but  will  continue  to  keep  his  sidewalls 
rigid.  The  hope  is  that  the  major  part  of  the  corrections  needed  for  normal  wlng-bndy  arrangements  can 
be  dealt  with  by  the  shaping  of  the  upper  and  lower  walls  and  any  residual  corrections,  hopefully  small, 
can  be  computed  with  adequate  accuracy  from  the  measured  pressure  distributions.  He  has  constructed  a 
6  in.  square  working  section  for  transonic  tests  and  two  wing/body  models  of  4  in.  span  have  been  made. 
These  have  been  tested  at  Langley  in  the  7  ft.  x  10  ft.  transonic  tunnel  to  provide  comparative  free  air 
results.  The  tests  may  range  up  to  stalling  incidences.  In  his  work  to  date  he  has  reached  a  Mach 
number  of  0.94. 

It  will  be  clear  that  self-streamlining  walls  offer  the  possibility  of  attaining  higher  Reynolds 
numbers  in  a  wind  tunnel  of  given  size  than  would  otherwise  be  feasible,  and  their  use  in  a  cryogenic 
tunnel  is  particularly  attractive. 

3.  WORK  AT  BRITISH  AEROSPACE,  KINGSTON-BROUGH  DIVISION4 ,5,0N 

TWO-DIMENSIONAL  MODELS  IN  A  WIND  TUNNEL  WITH  SOLID  WALLS. 

The  work  is  based  on  the  use  of  a  conformal  transformation  (c  =  ?  +  in  «  ie”2^)  which  transforms 
the  region  between  the  tunnel  floor  and  roof  in  the  z  plane  on  to  the  upper  half  of  the  %  plane.  Here  n 
is  the  tunnel  height  and  the  real  axis  in  the  z  plane  i;  along  the  tunnel  centre  line.  In  the  ;  plane 
the  wall  effect  is  therefore  reduced  to  that  of  a  single  image  in  the  f  axis  of  the  transformed  model. 

The  uniform  flow  in  the  z  plane  transforms  to  the  flow  ■frOm  a  source  at  the  origin  in  the  c  plane. 

The  model  is  represented  by  a  system  of  linearly  distributed  vortices  and  uniformly  distributed 
sources  over  the  surface  of  the  model,  the  latter  are  s^metrically  disposed  as  described  in  the  metioi 
of  Maskew  A  Woodward6.  The  boundary  layers  are  represented  by  a  normal  transpiration  velocity  at  tht 
surface  related  to  the  rate  of  change  with  distance  along  the  surface  of  the  displacement  thickness. 

The  wake  ir.  similarly  represented  by  linear  distributions  of  vortices  and  sources  over  the  interface 
of  the  two  wake  layers  extending  downstream  from  the  trailing  edge  and  lying  in  the  local  streamline 
direction.  An  allowance  for  the  boundary  layers  on  the  walls  has  yet  to  be  included.  An  iterative 
process  of  solution  is  adopted  to  bring  Into  conformity  a.he  inviscid  and  viscous  parts  of  the  flow, 
for  the  latter  Thwaites'  method  is  used  for  the  laminar  boundary  layers  and  the  lag-entrainment  method 
of  Green  et  al7  is  used  for  the  turbulent  boundary  layers  and  wakes. 

The  method  is  still  under  development  and  the  program  so  far  completed  has  yet  to  include  the  wake 
contribution;  it  has  been  applied  to  a  single  aerofoi1  op  to  stalling  incidence  for  which  experimental 
results  at  two  Reynolds  numbers  are  available  for  comparsbn,  with  the  ratio  of  h/c  about  6.9.  The 
calculated  curves  of  Cl  and  Cm  LE  as  functions  of  incideice  are  about  6 %  and  83!,  respectively,  above 
the  experimental  data  and  much  closer  to  the  latter  than  the  results  for  inviscid  flow.  The  inclusion 
of  the  wake  effects  is  expected  to  improve  the  agreement  with  the  experimental  data.  The  calculations 
yield  values  of  Co  that  are  also  in  reasonable  accord  with  the  experimental  data. 

It  is  hoped  to  extend  this  work  to  multi-component,  high  lift  arrangements  for  which  an  existing 
wall  jet  program  could  be  adapted  to  deal  with  the  use  of  blowing  as  well  as  the  interaction  of  a  wake 
from  one  component  and  the  boundary  layer  of  the  following  component. 

The  inviscid  program  has,  however,  also  been  applied  to  calculate  the  lift  and  pitching  moment 
characteristics  of  a  wing  in  free  air  and  in  a  tunnel  with  a  chord/tunnel  height  ■  0.3,  for  the 
following  configurations: 

11  unflapped 

ii)  303!  trailing  edge  flap  at  30° 

Hi)  303!  trailing  edge  flap  at  45° 

iv)  303!  trailing  edge  flap  at  60°  plus  a  103!  leading  edge  flap  at  45°. 

In  each  case  the  corrections  (Aa,  aCl,  aCmj)  to  a,  Cl,  Cmj*  were  determined  for  the  calculated  tunnel 
results  using  the  classical  standard  linearised  theory  formulae  derived  from  AGARDograph  1098,  and  the 
results  for  Cl  and  Cmj  calculated  for  the  wing  in  free  air  at  the  'corrected'  incidence  were  compared  with 
the  corresponding  'corrected'  tunnel  values  (i.e.  Cl  +  aCl,  Cm  +  aCm).  Some  of  the  results  are  set  out  in 

Table  1  where  suffix  'error'  denotes  the  difference  between  the  calculated  free  air  values  of  Cl  and  Cm 

and  the  corresponding  'corrected'  tunnel  values.  The  ratio  of  this  difference  to  the  corresponding 
standard  theory  correction  is  a  measure  of  the  validity  of  the  latter.  It  will  be  seen  that  for  the 
unflapped  section  at  low  incidences  the  standard  theory  corrections  are  of  acceptable  accuracy,  but  for 
all  other  cases  they  are  not  adequately  accurate.  Fig.  3  shows  that  the  uncorrected  results  can  be 
closer  to^free  air  results  than  after  correction  by  the  standard  method!  Examination  of  the  pressure 
distributions  shows  that  significant  changes  in  the  loading  distribution  can  be  induced  by  the  tunnel 
walls  and  these  do  not  lend  themselves  to  any  simple  correlation  which  could  hopefully  lead  to  seme 
generally  applicable  improvement  of  the  method. 


*  denotes  the  pitching  moment  coefficient  about  the  quarter-chord  point. 


7-3 


4.  THE  WORK  OF  ASHILL  &  WEEKS9  FOR 

TWO-DIMENSIONAL  MODELS  IN  THE  R.A.E.  8'  x  8'  TUNNEL. 

This  work  is  directed  at  the  corrections  needed  for  aerofoils  at  subsonic  speeds  in  a  tunnel  with 
solid  walls.  It  starts  with  the  argument  that  within  the  tunnel  the  Prandtl-Glauert  equation  applies 
in  a  region  S  outside  a  boundary  a  enclosing  the  aerofoil,  its  wake  and  areas  of  supersonic  flow 
(Fig.  4).  In  S  the  usual  Prandtl-Glauert  transformation  will  lead  to  Laplace's  equation  and  the 
complex  perturbation  velocity  can  be  related  by  Cauchy's  theorem  to  an  integral  over  the  boundaries  of  S. 

By  considering  an  aerofoil  in  free  air  with  the  same  flow  over  j  as  that  of  the  aerofoil  in  the  tunnel, 

the  difference  between  them  can  be  ascribed  to  the  integral  over  the  tunnel' roof  and  floor.  In  this  way 

the  tunnel  interference  can  be  determined  from  the  measured  velocity  components  over  the  roof  and  floor. 

With  solid  walls  tht  measured  longitudinal  component,  or  the  pressure  distribution,  at  the  walls  provides 
the  basis  for  determining  the  tunnel  interference.  The  analysis  can  be  applied  to  a  porous  wall  tunnel 

but  then  the  normal  velocity  distribution  is  also  required.  It  will  be  evident  that  the  method  will 

cease  to  apply  if  a  shock  extends  from  the  model  to  a  tunnel  wall.  Simple  extrapolation  formulae  have 

been  derived  to  allow  for  the  finite  length  of  working  section  over  which  it  is  possible  to  make 

measurements. 

The  method  has  been  applied  to  some  measurements  on  aerofoil  RAE  5225  with  c/h  =  0.26  at  various 
Mach  numbers  and  Cl  =  0.5  approximately  and  for  various  lift  coefficients  at  a  Mach  number  of  0.73. 

Fig.  5  shows  the  distributionof  wall-induced  upwash  over  the  aerofoil  chord  determined  from  the  measured 
wall  pressures  by  the  above  method  (vortex  distribution  method)  compared  with  the  distributions  given  by 
two  versions  of  the  conventional  image  system  method.  In  the  first,  of  these,  a  distribution  of  vorticity 
is  placed  along  the  chord  line  (linearised  theory,  dashed  line)  while  the  second  is  an  approximation  to 
the  first,  being  accurate  to  second  order  in  c/Beh  (linearised  theory,  full  line),  where  Be  =  (1  -  Me^)2 
and  Me  is  effective  free  stream  Mach  number.  The  agreement  betweei  the  three  methods  is  reasonably  good 
except  towards  the  trailing-edge  where  the  prediction  of  the  approximate  version  of  the  image  method 
deviates  from  those  of  the  other  two  methods.  The  difference  in  degrees  between  the  induced  upwash 
angles  at  the  leading  and  trailing  edges  =  0.5  Cl/B6. 

These  results  clearly  provide  a  warning  that  the  conventional  approach  of  treating  tunnel  Induced 
upwash  as  equivalent  to  an  incidence  change  will  be  increasingly  in  error  as  either  Cl  or  Mach  number 
increases.  Furthermore,  it  may  be  expected  that,  In  the  transoni  1  speed  range  where  there  is  increased 
sensitivity  to  local  effective  changes  of  shape,  the  effect  on  tne  pressure  distribution  of  the  chordwise 
variation  in  upwash  is  serious. 

Some  comparisons  have  been  made  between  the  measured  pressure  distributions  at  various  Mach  numbers 
and  corresponding  distributions  calculated  using  a  program  developed  by  Collyer  and  Lock10  (the  so-called 
VGK  method  -  a  version  of  the  Garabedian  &  Korn  method11  coupled  With  an  allowance  for  viscous  effects 
using  the  lag-entrainment  method  of  Green  et  al7  for  calculating  the  development  of  the  boundary  layer). 
For  the  latter  distributions  the  calculations  were  made  for  the  wing  section  with  and  without  corrections 
for  effective  camber  changes  associated  with  the  upwash  distribution  given  by  the  approximate  version  of 
the  image  method.  Some  results  are  illustrated  in  Fig.  6  for  M  =  1,504  and  in  Fig.  7  for  M  -  0.749. 

They  are  not  conclusive  in  indicating  any  clear  improvement  due  tr  allowing  for  these  camber  changes, 
but  for  these  cases  the  differences  are  small  and  not  easily  distinguished  from  possible  experimental 
errors  (including  side-wall  effects)  and  errors  in  the  theories  use1. 

The  method  for  determining  tunnel  wall  interference  from  measurements  of  wall  pressures  can  in 
principle  be  extended  to  high  lift  aerofoils  as  well  as  models  in  three  dimensions  (by  use  of  Green's 
theorem  provided  there  is  an  enveloping  region  in  which  the  Prandtl-Glauert  equation  holds).  Work  on 
these  lines  has  not  yet  gone  beyond  the  conceptual  stage. 

5.  WORK  AT  BRITISH  AEROSPACE,  WEYBRIDGE. 

VORTEX-LATTICE  REPRESENTATION.  WAKE  CURVATURE  EFFECT. 

12 

The  method  adopted  was  that  due  to  Joppa  '  in  which  the  model  wing  is  simply  represented  as  a 
lifting  line  with  a  pair  of  trailing  vortices,  but  these  are  not  assumed  rectilinear  but  are  allowed  to 
curve  in  response  to  the  flow  field.  The  tunnel  walls  are  modelled  hy  a  vortex-lattice  distribution  as 
illustrated  in  Fig.  8.  The  trailing  vortex  system  and  associated  flow  field  are  first  calculated  for 
unrestricted  flow  ('free  air')  by  an  iterative  process  and  then  the  calculation  is  repeated  for  the 
system  in  the  wind  tunnel.  The  differences  between  the  two  sets  of  calculations  for  the  same  wing 
span  and  wing  circulation  are  interpreted  as  the  interference  effects  and  the  corrections  for  wind 
tunnel  constraint  are  derived  from  them. 

The  method  was  applied  to  two  cases:  The  first  was  for  a  1/15  scale  VC10  (span  =  9.34  ft)  in  the 
Keybridge  13  ft  x  9  ft  low  speed  tunnel.  From  the  calculations  the  interference  factors  5  and  6'  for  wing 
and  tail  plane  incidences,  respectively,  were  determined.  These  are  related  to  the  corresponding 
incidence  and  pitching  moment  corrections  by:- 

Aa  =  4  £  CL  and  ACj^  =  -  4y  *  f 

C  h  C  5  n 

Here  S  =  reference  area  (i.e.  wing  area),  C  =  tunnel  cross-sectional  area,  t  =  tail  arm,  h  =  tunnel 
height,  n  =  tail  plane  incidence,  a  and  n  are  in  radians. 

For  this  model  the  following  results  were  obtained:- 

CL  4  4' 

1.1  0.105  0.185 

1.9  0.106  0.182 
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The  classical,  conventional  approach  of  treating  the  wake  vortices  as  rectilinear  and  summing  for 
the  images  In  the  tunnel  walls  gives  for  all  values  of  CL 

S  =  0.104  and  s'  =  0.178 

It  is  clear  that  for  this  case  the  differences- between  the  conventional  method  and  the  Joppa 
method  are  practically  negligible  for  the  values  of  Cl  considered.  However,  for  a  similar  sized  model 
generating  much  higher  values  of  Cl  the  difference  in  pitching  moment  correction  could  become  large 
enough  not  to  be  safely  ignored. 

The  second  case  considered  was  that  of  a  large  rectangular  wing  model  with  end  plates  used  for 
Investigations  of  high  lift  arrangements.  This  model  hid  a  span  of  7.5  ft  and  a  chord  of  2.5  ft  and 
was  also  tested  in  the  13  ft  x  9  ft  tunnel. 

For  .his  model  the  Joppa  programme  resulted  in  the  following  values  of  5 


Cl  5 


1.0 

C.  105 

4.0 

0.116 

5.4 

0.126 

The  conventional  method  gave  S  =  0.104,  for  all  values  of  Cl  i  here  there  Is  clearly  a  major 
difference  in  'the  predicted  interference  effects  for  the  higher  values  of  Ci  considered.  This  Is 
illustrated  in  Fig.  9  and  10  showing  the  resulting  curves  for  Cl  as  a  function  of  a  and  Cp  as  a  function 
of  Cl,  respectively. 

It  should  be  noted  that  for  this  model  4a,  the  correction  to  a,  is  of  the  order  of  Cl  in  degrees  a.id 
it  is  questionable  whether  with  such  large  corrections  the  chordwise  variations  of  4a  resulting  In 
effective  changes  of  model  shape  can  be  ignored. 

6.  WORK  BY  TREBBLE,  RAE,  ON  LIFTING  JETS13. 

LlitI  TATIONS  DUE  TO  TUNNEL  INTERFERENCE. 


T-eible  investigated  the  flow  field  associated  with  a  simple  fuselage  model  with  a  vertically 
emerging  jet  in  both  the  24  ft  tunnel  and  in  the  No.  1  11}  ft  x  8}  ft  tunnel  at  the  RAE  (see  Fig.  11). 

The  tesis  in  the  latter  tunnel  covered  a  range  of  free  stream  to  jet  dynamic  pressures,  qo/qj,  and 
various  ratios  of  the  height  of  the  nozzle  exit  above  the  tunnel  floor  (ii)to  the  nozzle  diameter  (d)  and 
they  induced  pressure  measurements  on  the  floor  of  the  tunnel.  The  model  size  was  such  that  the  24  ft 
tunnel  res  (Its  could  be  taken  as  free  of  significant  tunnel  interference,  and  they  were  therefore  used  is 
a  basis  for  assessing  the  interference  effects  in  the  smaller  tunnel.  Earlier  tests  as  well  as  the  flow 
visualisation  experiments  in  these  tests  have  shpwn.that  it  low  values  of  qj/q0  the  jet  is  diverted 
downstream  in  a  fairly  orderly  manner  by  the  main,  stream.  However,  as  qj/q0  is  Increased  the  jet  comes 
closer  to  the  tunnel  floor  and  eventually  at  some  critical  value  of  qj/q0  it  interacts  with  the  tunnel 
floor  boundary  layer  at  some  distance  downstream  of  the  jet  exit.  This  results  in  major  changes  to  the 
flow  in  the  boundary  layer  and  to  the  surface  pressure  distribution.  At  still  greater  values  of 
q j/q0  part  of  the  jet  flow  spreads  laterally  and  upstream  for  some  distance  along  the  tunnel  floor  and 
is  associated  with  the  development  of  a  horseshoe  vortex  which  moves  upstream  and  becomes  stronger  with 
increase  of  qj/q0. 

Clearly  this  critical  value  of  q j/q0  for  a  given  model  in  a  tunnel  defines  a  value  above  which  the 
flow  In  the  tunnel  is  sufficiently  different  from  that  of  flight  away  from  the  ground  that  no 
correlation  by  means  of  simple  wind  tunnel  corrections  is  possible.  Fig.  12  illustrates  the  nature  and 
degree  of  the  changes  that  can  occur  in  the  downwash  in  the  vicinity  of  the  model  fuselage  with 
variations  of  R/d  for  a  given  value  of  qo/qj-  In  this  case  critical  conditions  occur  for  h/d  about  4.8. 


Trebble's  measurements  show  that  the  onset  of  critical  conditions  can  be  fairly  readily  identified 
whether  one  is  looking  at  dynamic  pressure  or  downwash  (Fig. 13),  and  they  can  be  similarly  inferred  with 
reasonable  consistency  from  tunnel  floor  pressure  measurements.  From  an  analysis  of  his  own  and  other 


workers 1 
Fi/d  is 


results  he  found  that  a  good  empirical  relation  for  the  critical  value  of  qo/qj  as  a  function  of 

qQ/ q j  -  0.25  (d/R)  4/3 


for  the  jet  discharging  normally  to  the  main  stream  flow. 


Using  a  semi -empirical  formula  for  the  shape  of  an  emerging  jet  in  unrestrained  flow,  Trebble  argued 
tentatively  that  for  a  jet  emerging  at  an  angle  to  the  main  stream  flow  the  corresponding  relation 
would  be  ,,,  „  *  ■  _  ■? 

(qQ/ q j  3/4  =  0.145  (cot50°  -  cot  )  (d/R)^  sin  , 

but  there  is  as  yet  no  experimental  result  to  test  this. 


Trebble  also  did  some  tests  with  a  wing  attached  to  his  model  and  found  that  the  critical  conditions 
were  reflected  in  dramatic  changes  of  lift  associated  with  similar  changes  of  downwash,  and  these  critical 
conditions  were  reliably  predicted  by  the  above  relation  for  a  normal  jet  in  the  absence  of  a  wing.  He 
did  some  preliminary  work  on  an  open  jet  tunnel  and  found  similarly  critical  effects  when  the  jet  impinged 
on  the  edge  of  the  mainstream  flow. 


This  work  clearly  needs  pursuing  further  to  examine  experimentally  what  happens  with  a  jet  inclined 
to  the  main  stream  and  to  analyse  the  interference  effects  for  values  of  q0/q j  greater  than  the  critical. 
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7.  B.fle.  HATFIELD.  ..  ,,. 

WORK  ON  REVERSED  THRUST  BLOCKAGE  CORRECTIONS14’ l3‘ 

Measurements  r.ade  on  a  HS-125-700  model  in  a  7ft  x  9ft  low  speed  tunnel  showed  changes  in  the 
pressure  distributions  on  the  tunnel  walls  when  reverse  thrust  was  simulated.  These  changes  were 
presumed  to  be  associated  with  changes  in  blockage  effects  which  need  to  he  assessed. 


The  method  adopted  to  do  this  was  briefly  as  follows.  The  pressures  were  measured  on  the  tunnel 
roof  and  along  the  wall  and  groundboard  close  to  their  junction.  The  pressure  distributions  with  and 
without  the  model  (but  with  the  nacelles  and  rig)  were  determined  for  various  mass  flows,  and 
differences  were  derived  due  to  the  various  reverse  thrust  conditions  tested.  It  was  assumed  that  In 
each  case  the  pressure  difference  distribution  could  be  associated  with  an  equivalent  source-sink 
distribution  along  the  tunnel  axis  which  was  determined  from  the  measured  pressure  difference 
distribution  (the  mean  being  taken  between  roof  and  wall)  by  means  of  J.Y.G.  Evans'  method'6  and 
influence  coefficients.  The  method  then  yields  a  distribution  of  interference  velocity  increments  due 
to  blockage  along  the  tunnel  centre-line.  The  results  showed  marked  variation  of  this  increment  along 
the  centre-line  and  it  was  decided  to  fix  on  a  point  half  way  between  the  nacelle  exit  and  the  mid-chord 
position  for  determining  a  representative  value. 


Attempts  at  seeking  a  collapse  of  the  results  for  different  mass  flows  (or  reverse  thrusts)  showed 
a  near  linear  relation  could  be  established  between  e  (the  blockage  correction)  and  Cy  (a  non-dimensional 
thrust  coefficient)  but  with  some  scatter  (see  Fig. 14). 

8.  CONCLUDING  REMARKS 


It  seems  clear  that  the  use  of  suitable  programs,  based  on  some  form  of  panel  representation  of 
models  and  tunnel  walls  or  on  measurements  of  pressure  distrlbut  ons  on  the  walls,  can  up  to  a  point  lead 
to  significantly  more  accurate  estimates  of  tunnel  corrections  than  the  conventional  formulae  derived 
from  the  consideration  of  simplified  image  systems.  However,  su'h  methods  must  reach  their  limits  of 
application  at  the  point  where  the  corrections  are  such  that  they  can  no  longer  be  adequately  represented 
as  equivalent  to  the  effects  of  a  small  change  of  incidence  couM^d  with  a  simply  represented  small  change 
of  camber.  Beyond  such  limits  the  corrections  become  large  enough  and  vary  sufficiently  over  the  body 
surface  to  be  equivalent  to  a  change  of  shape  which  cannot  be  described  as  small  enough  for  the  effects 
to  be  regarded  as  additive  and  easily  calculated.  At  that  point  in  trying  to  determine  the  wind  tunnel 
corrections  we  are  faced  with  the  need  to  calculate  the  flow  over  a  body  of  prescribed  shape.  If  we 
could  do  that  accurately  we  would  not  have  needed  the  wind  tunnel  tests  in  the  first  place. 

The  only  solution  at  that  stage  appears  to  be  the  use  of  adaptive  walls.  It  remains  for  future 
research  to  show  how  closely  we  need  adjust  the  walls  to  the  fln-'l  streamline  shapes  of  unlimited  flow. 

It  may  well  be  that  judicious  use  of  fairly  simple  flexible  arrangements  coupled  with  measurements  of 
wall  pressure  distributions  to  provide  residual  corrections  will  prove  adequate  for  most  purposes. 

Even  the  use  of  adaptive  walls  Is  limited  to  cases  where  strong  shocks  do  not  extend  to  the  tunnel 
walls,  but  the  use  of  non-linear  transonic  flow  theory  for  the  asternal  flow  may  ease  this  limitation. 
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TUNNEL  FLOOR 

Fig.  4  SKETCH  ILLUSTRATING  BOUNDARY  AND 
REGION  S  IN  WHICH  P-G  EOU*  APPLIES 
(ASHILL  %  WEEKS  REF.  9) 


Fig. 6  RAE  5225  PRESSURE  DISTRIBUTIONS 
COMPARISON  BETWEEN  VGK  THEORY 
AND  MEASUREMENT  M.-0-504,  CL-0-551 
R  -20*106 


Fig.  5  WALL-INDUCED  UPWASH  AT  AEROFOIL 
CHORD  LINE  RAE  5225,  M.- 0-749 
C'0-560,  C.-0-238,  Rc-20x10t 


Flg.7  RAE  5225  PRESSURE  DISTRIBUTIONS 
COMPARISON  BETWEEN  VGK  THEORY 
AND  MEASUREMENT  M.-0-749  Cf0  557 
R.*20*106 


GIVEN  THE  LATTICE  ORDINATES, 
PROGRAM  'WINGT'  CALCULATES 
THE  CONTROL  POINT  ORDINATES 


Fig.  8  REPRESENTING  13*9'  LOW  SPEED  WIND  TUNNEL  BY  VORTEX 
LATTICE.'  (Square  vortex  rings  of  side  1-9 ft.) 


Fig  9  TYPICAL  HIGH  LIFT  MODEL  RESULTS.' 
THE  EFFECT  OF  APPLYING  JOPPA'S 
.INTERFERENCE  CORRECTION 


Fig.  10  TYPICAL  HIGH  LIFT  MODEL  RESULTS!  THE  EFFECT  OF  APPLYING 
JOPPA'S  INTERFERENCE  CORRECTION 
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TUNNEL  SPAN  «  3' 50  metres  i  d  .203  mm 


RAE  LIFTING  JET  MODEL  [see  Trebble,  Ref.  t3) 


Fig.  13  VARIATIONS  IN  TUNNEL  CONSTRAINT  WITH  VELOCITY  RATIO  FOR 
SEVERAL  STREAMWISE  LOCATIONS 
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Fig. 14  BLOCKAGE  CORRECTION  FACTOR  £  AS  FUNCTION 
OF  REVERSE  THRUST  COEFFICIENT 
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